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DESCRIPTION OF THE RELATED ART 



1 Communications Background 

1.1 Twisted pair telephone lines 

Telephone semce is provided to most businesses and homes via a pair of copper 
wires (a "twisted pair"). A telephone cable contains many twisted pairs: 25 twisted pairs 
are grouped in close proximity into 2 "binder groups," and several binder groups are 
packed together to form a cable. The two terminations of a telephone cable are at the user 
(subscriber) end and at the telephone company (central office, CO) end. We will use the 
terms "twisted pair," "line," and "subscriber loop" interchangeably herein as one example 
of a communications channel. 

Voice telephony uses only the first 4 kHz of bandwidth available on the lines. 
However, one can modulate data to over 1 MHz with significant bit rates. Only recently 
have schemes been developed to exploit the additional bandwidth of the telephone 
channel. A plot of the frequency response of a typical telephone channel is given in FIG. 
1. 

1.2 Overview of services 

In the past few years, a number of services have begun to crowd the bandwidth of 
the telephone channel. Some of the important services are: 

POTS - "Plain Old Telephone Service." This is the basic telephone service carrying 
voice traffic in the 0 - 4 kHz bandwidth. Conventional analog modems also use 
the same bandwidth. 

ISDN - Integrated Services Digital Network. This service allows end-to-end digital 
connectivity at bit rates of up to 128 kbps (kilo-bits-per-second). 
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Tl - Transmission !• This is a physical transmission standard for twisted pairs that 
uses 24 multiplexed channels (each at 64 kbps) to give a total bit rate of 1.544 
Mbps (Mega-bits-per-second). It uses costly repeaters. 

HDSL - High bit-rate Digital Subscriber Line. This is a full-duplex (two-way) Tl- 
like (1.544 Mbps) signal transmission service using only two twisted pairs and no 
repeaters. 

ADSL - Asymmetric Digital Subscriber Line. Over one twisted pair, this service 
provides a high-speed (on the order of 6 Mbps) downstream (from central office 
(CO) to subscriber) channel to each user and a low-speed (on the order of 640 
kbps) upstream (from subscriber to the central office) channel. This service 
preserves the POTS service over a single twisted pair. 

VDSL - Very high bit-rate DSL. This yet-to-be-standardized service will provide a 
very high speed (on the order of 25 Mbps) dovmstream channel to subscribers and 
a lower speed upstream channel to the central office over a single twisted pair less 
than 3 to 6 kft long. Further, it will preserve the POTS service. 

HDSL2 - High bit-rate Digital Subscriber Line 2. This soon-to-be-standardized 
service will provide ftiU-duplex 1.544 Mbps signal transmission service in both 
directions (full duplex) over a single twisted pair (<18 kft long) v^thout repeaters. 

"GDSL" - General Digital Subscriber Line. This hypothetical service would (for 
illustration purposes) carry 25 Mbps fiill-duplex data rate over a single twisted 
pair (see Sections 2.2.2 and 4.6.10). 

^*VDSL2" - Very high bit-rate DSL Line 2. This hypothetical service would (for 
illustration purposes) carry 12.4 Mbps fiill-duplex data rate over a single twisted 
pair less than 3 to 6 kft long (see Sections 2.2.3 and 4.6.10). 

Currently, all the above mentioned services have an ANSI standard except for 
VDSL, HDSL2, "GDSL" and "VDSL2". The various DSL standards (such as generic 
DSLs, ADSL, VDSL, HDSL2, and others) are collectively referred to as xDSL. 
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1.3 Crosstalk interference 
1.3.1 NEXT and FEXT 

Due to the close proximity of the lines within a binder, there is considerable 
amount of crosstalk interference between different neighboring telephone lines. 
Physically, there are two types of interference (see in FIG, 2): 

Near-end crosstalk (NEXT): Literference between neighboring lines that arises 
when signals are transmitted in opposite directions. If the neighboring lines carry 
the same type of service then the interference is called self-NEXT; otherwise, we 
will refer to it as different-service NEXT. 

Far-end crosstalk (FEXT): Interference between neighboring lines that arises when 
signals are transmitted in the same direction. If the neighboring lines carry the 
same type of service then the interference is called self-FEXT; otherwise, we will 
refer to it as different-service FEXT. 

FIG. 3 shows that crosstalk interference can be modeled as additive interference. 
Since neighboring lines may carry either the same or a different flavor of service, there 
are three categories of interference (see FIG. 3): 

1. Self-interference (self-NEXT and self-FEXT) between Unes carrying the same 
service. 

2. Interference into a channel carrying service A from other lines carrying services 
other than A (DSIN-NEXT and DSIN-FEXT). 

3. Interference from a channel carrying service A into other lines carrying services 
other than^ (DSOUT-NEXT and DSOUT-FEXT). 

Chaimel noise will be modeled as additive Gaussian noise (AGN). 
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1.3.2 Notation for self-NEXT and self-FEXT 

Here is some notation to keep things clear. Number the M twisted pairs (lines) in 
the cable with index /g { 1 , . , . , M}, and denote the direction of transmission with index o 
e{w5 c/}, with u = upstream (to the central office) and d = downstream (from the central 
office). All the twisted pairs in the cable bundle are assumed to carry the same service. 
Let o be the complement direction of o: u-d, d-u. Denote the transmitters and 
receivers on line / as: 

Z)^ : transmitter (Tx) on twisted pair / in direction o. 

i?f : receiver (Rx) on twisted pair / in direction o. 

Ideally, intends to transmit information only to R%]ndi real system, however, 
's signal leaks into the receivers i?J and i?J . Using our notation, this self-interference 
corresponds to: 

Self-NEXT: Crosstalk from into R] for all j ^ i, o e {u, d), and 
Self-FEXT: Crosstalk from into for all j ^ i, o e { w, . 

In a full-duplex xDSL service, each twisted pair i supports transmission and 
reception in both directions (using echo cancelers), so each line / has a foil set of 
transmitters and receivers: {7]" , R% Tf , Rf}. With perfect echo cancellation, there is 

no crosstalk from into R^^ . We will assume this for the balance of this document, 

although this crosstalk could be dealt with in a fashion similar to self-NEXT and self- 
FEXT. 

1.4 Capacity and performance margin 

The Channel capacity C is defined as the maximum number of bits per second 
that can be transmitted over a channel with an arbitrarily small bit error probability. The 



Conley, Rose & Tayon, P.C. 



9 



Atty. Docket No. 1789-01910 



achievable rate for a channel is any transmission rate below or equal to capacity, i.e., 
< C. Another channel performance metric is performance margin (or margin). It is 
defined (in dB) as 

margin = 101og„[^^ 

\ nun / 

5 where iSTVRrec is the received signal-to-noise ratio (SNR) and SNRmm is the minimum 
received SNR required to achieve a fixed bit error probability (BER) at a given 
transmission rate. The performance margin of a channel for a fixed bit error probability 
measures the maximum degradation (fi"om noise and interference) in achievable bit rate 
that a channel can sustain before being unable to transmit at that bit rate for a fixed BER 

10 (see [12]). The higher the performance margin of a channel at a given transmission rate 
and fixed BER, the more robust it is to noise and interference, i.e., the better is its 
performance. 
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2 Problem Statement 

2.1 General statement 

Given an arbitrary communications channel with: 

1. Self-interference (self-NEXT and self-FEXT) between users of serviced, 

5 2. Interference from users of different services v^th users of service A (DSIN-NEXT 

and DSIN-FEXT), 

3. Interference from users of service A into users of different services (DSOUT- 
NEXT and DSOUT-FEXT), and 

4. Other interference (including noise), 

10 maximize the capacity of each user of service A without significant performance (capacity 
or margin) degradation of the other services. 

Here services could refer to different possible signaling schemes. Users refer to 
the generic Tx-Rx pairs. 



15 2.2 Particular statement for DSLs 
2.2.1 HDSL2 service 

As a special case of the general problem, we will look into a particular problem of 
subscriber loops. In particular, we can phrase our statement in the language of HDSL2 
[2]. Here, the communications channel is the collection of twisted pairs in the telephone 
20 cable, interference is caused by: 

1. Self-NEXT and self-FEXT between neighboring HDSL2 lines (self-NEXT 
dominates over self-FEXT [8]), 

2, DSIN-NEXT and DSIN-FEXT from Tl, ISDN, HDSL and ADSL, 
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3. Interference from HDSL2 into other services, such as Tl, ISDN, HDSL and ADSL, 

and 

4. Channel noise, which we will model as AGN. 

We wish to maximize the capacity of the HDSL2 service in presence of other 
5 HDSL2, Tl, ISDN, HDSL, ADSL, VDSL lines and even services not yet imagined while 
maintaining spectral compatibility with them. We will consider HDSL2 service in 
Sections 4.4 to 4.7. 

The HDSL2 service is intended to fill a key need for fast (1.544 Mbps) yet 
affordable Ml duplex service over a single twisted pair. Efforts to define the standard are 
10 being mounted by several companies and the TlEl standards committee. The two key 

issues facing HDSL2 standards committee are: 

Spectral optimization. All previously proposed schemes for HDSL2 achieve the 
required data rates with satisfactory margins only in complete isolation. 

However, due to the proximity of the hues in a cable, there is considerable 
15 DSINNEXT, DSIN-FEXT, self-NEXT and self-FEXT interference from Tl, 

ISDN, HDSL, ADSL and HDSL2 into HDSL2 ~ this interference reduces the 
capacity of the HDSL2 service. 

Simultaneously, there is considerable DSOUT-NEXT and DSOUT-FEXT 
interference from HDSL2 into Tl, ISDN, HDSL and ADSL. This problem is 
20 known as spectral compatibility. The scheme ultimately adopted for HDSL2 must 

not interfere overly with other DSL services like Tl, ISDN, HDSL, and ADSL. 

Modulation scheme* No prior system has been developed that systematically 
optimizes the HDSL2 spectrum and reduces interference effects both from and 
into HDSL2. Further, a modulation scheme for HDSL2 has not been decided upon 
25 at this time. 
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2.2.2 "GDSL" service 

Consider the hypothetical DSL service "GDSL" described above. The "GDSL" 
service v^U enable very high bit-rate full-duplex, symmetric traffic over a single twisted 
pair. We assume that the lines carrying GDSL service have good shielding against self- 
5 NEXT. In this case, interference is caused by: 

1. Self-NEXT and self-FEXT between neighboring "GDSL" lines (self-FEXT 
dominates over self-NEXT), 

2. DSIN-NEXT and DSIN-FEXT from Tl, ISDN, HDSL, HDSL2 and ADSL, 

3. Literference from "GDSL" into other services, such as Tl, ISDN, HDSL, HDSL2 
10 and ADSL, and 

4. Channel noise, which we will model as AGN. 

We wish to maximize the capacity of the "GDSL" service in presence of other 
"GDSL", Tl, ISDN, HDSL, ADSL, HDSL2 lines and even services not yet imagined 
while maintaining spectral compatibility with them. The spectral optimization issue is 
15 similar to the one discussed for HDSL2 case, and we need to find an optimal transmit 

spectrum for "GDSL". Further, a good modulation scheme needs to be selected. 



2.2.3 "VDSL2" service 

Consider the hypothetical DSL service "VDSL2" described above. Optical fiber 
20 Unes having very high channel capacity and virtually no crosstalk will be installed in the 
fiiture up to the curb of each neighborhood (FTTC). The final few thousand feet up to the 
customer premises could be covered by twisted pairs. In such a scenario, high bit-rate 
asymmetric-traffic services (like VDSL) and symmetric-traffic services (like "VDSL2") 
over short length twisted pairs would become important. For illustration of such a 
25 potential future service we propose a hypothetical "VDSL2" service that would carry very 
high bit-rate symmetric traffic over short distance loops on a single twisted pair. In the 
"VDSL2" case, the interference will be caused by: 
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1. Self-NEXT and self-FEXT between neighboring "VDSL2" lines {both self NEXT 

and self-FEXT are dominant), 

2. DSIN-NEXT and DSIN-FEXT from Tl, ISDN, HDSL, HDSL2, VDSL and ADSL, 

3. Interference from "VDSL2" into other services, such as Tl, ISDN, HDSL, HDSL2, 
5 VDSL and ADSL, and 

4. Channel noise, which we will model as AGN. 

Again, we wish to maximize the capacity of "VDSL2" in presence of all the other 
interferers. To achieve this we need to find optimal transmit spectra and a good 
modulation scheme. 

10 
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3 Previous Work 

Here we discuss prior work pertaining to HDSL2 service. 

3.1 Static PSD masks and transmit spectra 

The distribution of signal energy over frequency is knovra as the power spectral 
density (PSD). A PSD mask defines the maximum allowable PSD for a service in 
presence of any interference combination. The transmit spectrum for a service refers to 
the PSD of the transmitted signal. Attempts have been made by several groups to come 
up with PSD masks for HDSL2 that are robust to both self-interference and interference 
from other lines. One way of evaluating channel performance is by fixing the bit rate and 
measuring the performance margins [12]: The higher the performance margin for a given 
disturber combination, the more robust the HDSL2 service to that interference. The term 
crosstalk here implies self-interference plus interference from other lines. 

To the best of our knowledge, no one has optimized the PSD of HDSL2 lines in 
presence of crosstalk and AGN. The significant contributions in this area, MONET-PAM 
and OPTIS, [1, 2, 4, 5] suggest a static asymmetrical (in input power) PSD mask in order 
to attempt to suppress different interferers. The PSD masks suggested in [1, 2, 4, 5] have 
a different mask for each direction of transmission. Furthermore, the techniques in [1, 4] 
use different upstream and downstream average powers for signal transmission. However, 
the mask is static, implying it does not change for differing combinations of interferers. 

Optis [5] is currently the performance standard for HDSL2 service. 

When a constraining PSD mask is imposed, the transmit spectrum lies below the 
constraining mask. Specifying a constraining PSD mask only limits the peak transmit 
spectrum. We do PSDs (transmit spectra) and not masks in this document unless stated 
otherwise. In Section 4.1 1 we indicate ideas to get PSD masks. 
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3.2 Joint signaling techniques 

Self-NEXT is the dominant self-interference component in symmetric-data-rate, 
full-duplex, long-length Hne xDSL service (e.g., HDSL2). One simple way of completely 
suppressing self-NEXT is to use orthogonal signaling (for example, time division 
signaling (TDS), frequency division signaling (FDS), or code division signaling (CDS)). 
In TDSj vs^e assign different services to different time slots. In FDS, v^e separate in 
frequency the services that could interfere with each other. In CDS, a unique code or 
signature is used in each direction of service. Further, in CDS each service occupies the 
entire available bandwidth for all of the time. CDS is similar to code division multiple 
access (CDMA), but here instead of providing multiple access, CDS separates the 
upstream and dovmstream transmit spectra using different codes. 

The choice of orthogonal signaUng scheme depends on the intent. We will see that 
FDS is in a sense optimal under an average power constraint (see Section 4.5.12). 

To eliminate self-NEXT using FDS, we would force the upstream transmitters 
{ 7^" , z = 1, . . . , M},and the downstream transmitters { 7]"^ , z = 1, . . . , A/} to use disjoint 
frequency bands. Thus, in FDS signaling, the upstream and downstream transmissions are 
orthogonal and hence can be easily separated by the corresponding receivers. Since in a 
typical system FDS cuts the bandwidth available to each transmitter to 14 the overall 
channel bandv^dth, we have an engineering tradeoff: FDS eliminates self-NEXT and 
therefore increases system capacity; however, FDS also reduces the bandwidth available 
to each transmitter/receiver pair and therefore decreases system capacity. When self- 
NEXT is not severe enough to warrant FDS, both upstream and downstream transmitters 
occupy the entire bandwidth. In this case, the upstream and downstream directions have 
the same transmit spectrum; we refer to this as equal PSD (EQPSD) signaling. 

On a typical telephone channel, the severity of self-NEXT varies with frequency. 
Therefore, to maximize capacity, we may v^sh to switch between FDS and EQPSD 
depending on the severity of self-NEXT. Such a joint signaling strategy for optimizing 
the performance in the presence of self-NEXT and white AGN was introduced in [3]. 
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The scheme in [3] is optimized, but only for an over simplified scenario (and 
therefore not useful in practice). In particular, [3] does not address self-FEXT and 
interference from other lines as considered in this work. Further, [3] does not address 
spectral compatibility issue. 

All other schemes for joint signaUng employ ad-hoc techniques for interference 
suppression [1, 2, 4, 5]. 

3.3 Multitone modulation 

Multicarrier or discrete multitone (DMT) modulation [6] can be readily used to 
implement a communication system using a wide variety of PSDs. Multitone modulation 
modulates data over multiple carriers and adjusts the bit rate carried over each carrier 
according to the signal to noise ratio (SNR) for that carrier so as to achieve equal bit error 
probability (BER) for each carrier (see in FIG. 4). 

Orthogonal FDS signaling is easily implemented using the DMT: we simply 
assign transmitter/receiver pairs to distinct sets of carriers. Note, however, that multitone 
modulation is definitely not the only modulation scheme that can be used to 
implement (optimal) transmit spectra. We can just as well use other techniques, such 
as CAP, QAM, multi-level PAM, etc. 

3.4 Summary of previous work 

The current state of the art of DSL technology in general and HDSL2 in particular 
can be described as follows: 

• Ad-hoc schemes (sometimes referred to as "optimized") have been developed that 
attempt to deal with self-interference and DSIN-NEXT and DSIN-FEXT as well 
as spectral compatibility of the designed service with other services. However, 
these schemes by no means optimize the capacity of the services considered. 
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• An optimal signaling scheme has been developed in [3] for the case of self-NEXT 
and white additive Gaussian noise only. The development of [3] does not address 
crosstalk from other sources, such as DSIN-NEXT and DSIN-FEXT, or self- 
FEXT, or other types of additive Gaussian noise. The development of [3] also 
does not address spectral compatibility of the designed service with respect to 
other services. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

Other objects and advantages of the invention will become apparent upon reading 
the following detailed description and upon reference to the accompanying drawings in 
which: 

5 FIG. 1 is an example of the frequency-response for a twisted pair telephone channel; 

FIG. 2 shows NEXT and FEXT between neighboring lines in a telephone cable, with 
"Tx" and "Rx" indicating transmitters and receivers, respectively; 

FIG. 3 shows how NEXT (DSIN-NEXT and self-NEXT) and FEXT (DSIN-FEXT 
and self-FEXT) are modeled as additive interference sources, with DSOUT- 
10 NEXT and DSOUT-FEXT representing the interference leaking out into other 

neighboring services; 

FIG. 4 illustrates how multicarrier, or discrete multitone (DMT) modulation 
multiplexes the data onto multiple orthogonal carrier waves; 

FIG. 5 and FIG. 5A are representative views of a subscriber-line communications 
15 system and a well-logging system that use the present invention; 

FIG. 6 is a representative view of a home system using the present invention for DSL 
commimications ; 

FIG. 7 is a block diagram of one embodiment of the computer from FIG. 6; 

FIG. 8 is a block diagram of one embodiment of the DSL card from FIG. 7; 

20 FIG. 9 is a flowchart for determining transmission characteristics for a 

communications system in one embodiment of the invention; 

FIG. 10 is a flowchart for determining a transmit spectrum with preKminary analyses 
of self interference and FEXT levels; 

FIG. 1 1 is a flowchart for determining a transmit spectrum with preliminary analyses 
25 of self interference levels; 
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FIG. 12 is a flowchart for determining a transmit spectrum; 

FIG. 13 is a flowchart for method for transmitting data on a communications channel; 

FIG. 14 is a flowchart for initiating a data transfer on the conmiunications channel; 

FIG. 15 is a flowchart for determining transmission characteristics for a 
communications system in one embodiment of the invention; 

FIG. 16 is a frequency-response graph showing the channel sub-division into K 
narrow bins (subchannels), each of vddth fF(Hz); 

FIG. 17 shows the magnitude squared transfer function of the channel (CSA loop 6), 
with 39 self-NEXT interferers, and 39 self-FEXT interferers (see (l)-(3)); 

FIG. 18 shows transmit spectra for EQPSD, FDS and multi-line FDS signaling 
schemes in a single frequency bin k for the case where the number of lines is 3 
(this also works for any number of lines); 

FIG. 19 is a model for combined additive interference from other services (DSIN- 
NEXT and DSIN-FEXT) plus channel noise (AGN); 

FIG. 20 is a flowchart of a method for determining an optimal transmit spectrum 
using only EQPSD signaling; 

FIG. 21 is a graph of an optimal transmit spectrum of HDSL2 (on CSA loop 6) vvdth 
49 HDSL DSIN-NEXT interferers and AGN of -140 dBm/Hz; 

FIG. 22 is a graph of an optimal transmit spectrum of HDSL2 (on CSA loop 6) with 
25 Tl DSIN-NEXT interferers and AGN of -140 dBm/Hz; 

FIG. 23 shows upstream and downstream transmit spectra in a single frequency bin (a 
= 0.5 => EQPSD signaling and a = 1 ^ FDS signaling); 

FIG. 24 is a graph demonstrating that is monotonic in the interval ae (0.5, 1]; 

FIG. 25 shows EQPSD and FDS signahng in a single frequency bin; 
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FIG. 26 shows upstream and downstream transmit spectra with regions employing 
EQPSD signaling (in bins [1, Me2f]) and FDS signaling (in bins [Me2f + 1, K]); 

FIG. 27 is a flowchart of the optimal and suboptimal schemes to determine the 
transmit spectrum using EQPSD and FDS signaling (and EQPSD/FDS transmit 
spectrum); 

FIG. 28 shows joint EQPSD/FDS signaling for a channel with "discrete" and 
"contiguous" transmit spectra for upstream (top graphs) and downstream (bottom 
graphs) signaling; 

FIG. 29 is a graph of an optimal upstream transmit spectrum for CSA Loop 6 using 
HDSL2 with 39 self-NEXT and 39 self-FEXT interferers, with EQPSD signaling 
taking place to the left of bin 9 (indicated by solid line) and FDS signalmg taking 
place to the right (indicated by dashed line); 

FIG. 30 shows graphs of optimal "contiguous" upstream and downstream transmit 
spectra for CSA Loop 6 using HDSL2 with 39 self-NEXT and 39 self-FEXT 
interferers (EQPSD signaling taking place to the left of bin 9); 

FIG. 31 shows graphs of another set of optimal "contiguous" upstream and 
downstream transmit spectra for CSA Loop 6 using HDSL2 with 39 self-NEXT 
and 39 self-FEXT interferers, with the property that these spectra yield equal 
performance margins (equal capacities) and equal average powers in both 
directions of transmission (EQPSD signaling taking place to the left of bin 9); 

FIG. 32 shows transmit spectra of signaUng line (S), interfering line (Y and Z), and 
lumped channel noise (N) for two cases: the FDS scheme (Case 2) for interfering 
line yields higher capacity for signaling line (S) than other schemes like CDS 
(Case 1); 

FIG. 33 shows EQPSD and multi-Une FDS signaling in a single frequency bin k for 
theM=3 line case; 
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FIG. 34 shows FDS and multi-line FDS signaling in a single frequency bin k for the 
M=3 line case; 

FIG. 35 is an example of an upstream transmit spectrum of line 1 iS^(f)) employing 
EQPSD, FDS and multi-line FDS signaling schemes for the M = 3 line case, in 
which bins [1, MEmFDs] employ EQPSD, bins [Me2mfds + U Mmfds2fds] employ multi- 
line FDS, bins [Mmfds2fds + 1, Mfi)S2MFDs] employ FDS, and bins [Mjr^s2MFDs + 1. 
employ multi-line FDS; The downstream spectrum of line 1 (Sf(f))is similar to 

S^(f) except for putting power in the complementary halves of FDS bins; The 
upstream spectra of Unes 2 and 3 are similar to S^(f) except for putting power in 
complementary thirds of multi-line FDS bins; The downstream spectra for lines 2 
and 3 are similar to S^(f) except for putting power in the complementary halves 
of the FDS bins and in the complementary thirds of multi-line FDS bins; 

FIG. 36 illustrates practical observation 1, a case of an EQPSD/FDS/MFDS transmit 
spectrum in which there is no FDS spectral region; bins [1, Me2mfds\ employ 
EQPSD, and bins [Me2mfds + 1, ^ employ multi-line FDS; 

FIG. 37 illustrates practical observation 2, a case of an EQPSD/FDS/MFDS transmit 
spectrum in which there is no multi-line FDS spectral portion within the EQPSD 
region; bins [1, M^eds2fds] employ EQPSD, bins [Mmfds2fds + U Meds2mfds] employ 
FDS, and bins [Mfds2mfds + 1? ^ employ multi-hne FDS; 

FIG. 38 shows upstream and downstream transmit spectra in a single frequency bin (a 
= 0.5 => EQPSD signaling and a=l^ multi-line FDS signaling); 

FIG. 39 shows EQPSD and multi-line FDS signaling in a single frequency bin; 

FIG. 40 is a flowchart of a scheme for determining an optimal transmit spectrum 
using EQPSD, FDS, and multi-line FDS signalmg (an EQPSD/FDS/MFDS 
transmit spectrum); 
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FIG, 41 shows, for the case where the lines have different line characteristics, 
upstream and downstream transmit spectra in a single frequency bin (a = 0.5 ^ 
EQPSD signaling and a multi-line FDS signaling); 

FIG. 42 shows, for the case where the lines have different line characteristics, 
upstream and downstream transmit spectra in a single frequency bin (a = 0.5 =:> 
EQPSD signaling and a = 1 multi-Une FDS signaling); 

FIG. 43 is a graph of an optimal downstream transmit spectrum for HDSL2 (on CSA 
loop 6) imder an OPTIS downstream constraining PSD mask with 49 HDSL 
DSIN-NEXT interferers and AGN of -140 dBm/Hz (the 'o-o' line shows the 
peak-constrained optimal transmit spectrum and the line shows the 
constraining OPTIS PSD mask); 

FIG. 44 is a graph of an optimal upstream transmit spectrum for HDSL2 (on CSA 
loop 6) under an OPTIS upstream constraining PSD mask with 25 Tl DSIN- 
NEXT interferers and AGN of -140 dBm/Hz (the 'o-o' Une shows the peak- 
constrained optimal transmit spectrum and the line shows the constraining 
OPTIS PSD mask); 

FIG. 45 shows graphs of optimal upstream and downstream transmit spectra for 
HDSL2 (on CSA loop 6) under the OPTIS upstream and downstream constraining 
PSD masks with 39 HDSL2 self-NEXT and self-FEXT interferers and AGN of - 
140 dBm/Hz (the 'o-o' lines show the peak-constrained optimal transmit spectra 
and the lines show the constraining OPTIS PSD masks); 

FIG. 46 shows graphs of optimal upstream and downstream transmit spectra for 
HDSL2 (on CSA loop 6) under the OPTIS upstream and downstream constraining 
PSD masks with 24 HDSL2 self-NEXT and self-FEXT interferers, 24 Tl 
interferers, and AGN of -140 dBm/Hz (the 'o-o' lines show the peak-constrained 
optimal transmit spectra and the lines show the constraining OPTIS PSD 
masks); 
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FIG. 47 shows graphs of optimal "contiguous" upstream and downstream transmit 
spectra for HDSL2 (on CSA loop 4, with a non-monotonic channel function due 
to bridged taps) with 39 HDSL2 self-NEXT and self-FEXT interferers; these 
transmit spectra yield equal performance margins (equal capacities) and equal 
average powers in both directions of transmission (note that there is only one 
transition region from EQPSD to FDS signaling); 

FIG. 48 shows (in the top graph) the channel transfer function, self-NEXT, and self- 
FEXT transfer functions for a short loop with bridged taps employing "GDSL" 
service (note that self-NEXT is very low for this hypothetical service), and shows 
(in the bottom graph) the distributed EQPSD and FDS spectral regions for the 
upstream and downstream transmit spectra, with a 0 indicating EQPSD signaling, 
a 1 indicating FDS, and a 0.5 indicating EQPSD or FDS signaling (note that in 
this case the non-monotonicity of the channel transfer function leads to several 
distributed signaling regions); 

FIG. 49 shows an alternative signaUng scheme: in the presence of high degrees of 
self-NEXT and self TEXT between group of Unes 1 and 2 and lines 3 and 4, we 
employ multi-line FDS; there is EQPSD signaling within each group of lines (1 
and 2 employ EQPSD as do 3 and 4) that have low self-interference; 

FIG. 50 shows optimal transmit spectra (upstream and downstream) of ADSL (on 
CSA loop 6) with 49 HDSL DSIN-NEXT interferers and an AGN of -140 
dBm/Hz; 

FIG. 51 shows optimal transmit spectra (upstream and downstream) of ADSL (on 
CSA loop 6) with 25 Tl DSIN-NEXT interferers and an AGN of -140 dBm/Hz; 
and 

FIG. 52 shows upstream and downstream transmit spectra in a single frequency bin (a 
= 0.5 ^ SPSD signaling and a = 1 =^ FDS signaling). 
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While the invention is susceptible to various modifications and alternative forms, 
specific embodiments thereof are shown by way of example in the drawings and will 
herein be described in detail. It should be understood, however, that the drawings and 
detailed description thereto are not intended to limit the invention to the particular form 
disclosed, but on the contrary, the intention is to cover all modifications, equivalents and 
alternatives falling within the spirit and scope of the present invention as defined by the 
appended claims. 
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DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 



The present invention comprises an improved system and method for 
communicating information such as voice, images, video, data, or other information on a 
transmission medium. The present invention provides improved communications on the 

5 transmission medium in the presence of interference. More specifically, the present 
invention operates to model and then minimize the effects of interference on the 
transmission medium. The interference may take the form of similar services being 
transmitted on neighboring transmission mediums and/or may take the form of 
uncorrected interference from different services on neighboring transmission mediums 

10 and or interference from various noise sources which affect the transmission medium. 

FIG. 5- Subscriber Line Embodiment 

FIG. 5 illustrates a preferred embodiment for use of the present invention. FIG. 5 
illustrates a first location, e.g., a home 10 that is coupled through a subscriber line 12 to a 
15 second location, e.g. a telephone company central office (CO) 14. It is noted that the first 

and second locations may be any of various types of sites, such as a home, an office, 
business, an office building, or a CO. 

In this embodiment of the invention, the communication system and method is 
comprised in a digital subscriber line (DSL) device that operates to perform xDSL 
20 communications on subscriber line 12. Thus, this figure shows a configuration that 

includes subscriber line 12, e.g., a twisted-pair copper line, that is coupled between home 
10 and CO 14. The present invention is comprised in each of home 10 and CO 14. 

As discussed in the background section, subscriber lines are generally included in 
a cable that has a plurality of closely positioned transmission mediums, including other 
25 subscriber lines. Due to the close proximity of the transmission mediimis comprised in a 
subscriber cable, a given subscriber line is subject to interference from neighboring 
transmission mediums, including self-NEXT and self-FEXT interference, and different 
service interference (DSIN). Some of the transmit spectra discussed herein are 
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substantially optimized to maximize perfomiance margins and avoid the effects of this 
interference, thereby providing improved communication. 

By design, an optimal transmit spectra give increased performance margins 
(increased immunity against noise) and spectral compatibility margins as compared to 
5 one fixed transmit spectrum. The optimal transmit spectra described herein are typically 
obtained by fixing the average input power and choosing the best signaling strategies and 
optimal power distribution to maximize the bit rate or performance margin. The transmit 
spectra may also be used to minimize the required average input power where the desired 
performance margins or bit rates are fixed. 

10 

FIG. 5A - Well Logging Embodiment 

FIG. 5 A illustrates an alternate scenario for use of the communication system and 
method of the present invention. FIG. 5A illustrates a drill hole and/or well-logging 
scenario which utilizes the communication system of the present invention. As an 
15 example, in FIG. 5 A commxmication equipment 16 on the surface communicates through 

a commimication medium 12A to instrumentation 18 comprised in the borehole 
underground. The communication system and method operates to reduce the effects of 
interference in the well hole and to provide improved communications. 

Although FIG. 5 and FIG, 5A illustrate two embodiments for use of the system 
20 and method of the present invention, it is noted that the present invention may be used in 

any of various types of systems or scenarios which involve communication of data on a 
transmission medium that is subject to noise or other interference. The present invention 
is particularly usefiil in scenarios where the transmission medium is in close proximity to 
various sources of interference that can be ascertained, identified, and modeled. In 
25 general, the present invention is applicable to reduce the effects of interference on 
transmission media that are subject to interference from known or unknown sources 
where the spectral characteristics of the interference can be modeled. 
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In the analyses presented herein, we use a generic DSL (xDSL) model. For 
concreteness, we present results optimizing DSL services such as HDSL2, "GDSL", and 
"VDSL2" in the face of noise and interference jfrom neighboring services. The invention 
is not, however, limited to these services, but can be applied to any communications 
chaimel that exhibits crosstalk interference. 

Although FIG. 5 illustrates a subscriber line embodiment, it is noted that the 
present invention may be used for any of various types of transmission media, e.g., 
copper wire, fiber optic, lines, co-axial cable, wave guides, etc. For example, the present 
invention is well suited for use in local and vdde-area networks to minimize noise 
interference on networks, e.g., Ethemet, token ring, and wide area networks such as frame 
relay, Svdtched 56, ATM (asynchronous transfer mode), etc. Also, although FIG. 5 
illustrates use of the present invention between a home 10 and a central office 14 over a 
subscriber loop or subscriber line 12, it is noted that the present invention may also be 
used for the various trunks comprised in the PSTN. The present invention is also useful in 
the various backbones or lines used for the Internet. 

The present invention is also useful for wireless transmission applications, e.g., 
cellular telephones, cordless telephones, short wave radio, etc. as well as the various 
broadcast media such as the various digital satellite services providing television, 
Internet, data or voice services. In short, the present invention is applicable to any of 
various types of systems which involve the transmission of data over a wired or wireless 
medium. The present invention is also appUcable to any of the various types of signaling 
protocols such as frequency division multiple access (FDMA), time division multiple 
access (TDMA) and code division multiple access (CDMA) as well as hybrids of these, 
among others. 

Therefore, the present invention is applicable to a variety of communications 
charmels in a nimiber of communication scenarios. In the description that follows, the 
present invention is described with respect to the preferred embodiment, the preferred 
embodiment being a digital subscriber line application between a first location, e.g., a 
home or business 10, and a telephone company central office 14. 



Conley, Rose & Tayon, P.C. 



28 



Atty. Docket No. 1789-01910 



FIG. 6 - Home DSL System 

FIG. 6 illustrates a system 100 comprised in location 10, i.e., in the home or 
business 10 which performs digital subscriber line (DSL) communication operations over 
subscriber line 12. In a preferred embodiment, the DSL circuitry of the present invention 
is comprised in a computer system 102 coupled to subscriber line 12 through an xDSL 
port 106. In one embodiment, computer system 102 is also coupled to a telephone system 
104. However, it is noted that the DSL system of the present invention may be comprised 
in any of various types of systems including computer systems, Internet appliances, 
televisions or dedicated boxes. In the preferred embodiment, the DSL system of the 
present invention is comprised in a DSL device on an add-in card to the general purpose 
computer system 102. In the preferred embodiment, the DSL card includes a port for 
coupling to a standard telephone jack, or "splitter," which in turn couples to the 
subscriber Une 12, In this embodiment, the computer system 102 may be utilized as a 
virtual telephone which operates through the DSL device for voice commimications over 
the subscriber line 12. In another embodiment, a separate telephone system 104 is 
coupled to a second port of the DSL card, as shown in FIG. 6. 

FIG. 7 - Computer System Block Diagram 

Turning now to FIG, 7, a block diagram of one embodiment of computer system 
102 is shown. Other embodiments are possible and contemplated. The depicted system 
includes a microprocessor or CPU 110 coupled to a variety of system components 
through a bus bridge 1 14 . In the depicted system, a main memory 1 12 is also coupled to 
bus bridge 1 14, Finally, a plurality of PCI devices are coupled to bus bridge 1 14 through 
a PCI bus 1 16. In the depicted embodiment, the PCI devices include a video card 118 and 
a add-in card for the DSL device 120. 

Bus bridge 114 provides an interface between microprocessor 110, main memory 
1 12, and the devices attached to PCI bus 1 16. When an operation is received from one of 
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the devices connected to bus bridge 114, bus bridge 114 identifies the target of the 
operation (e.g. a particular device or, in the case of PCI bus 116, that the target is on PCI 
bus 116). Bus bridge 114 routes the operation to the targeted device. Bus bridge 114 
generally translates an operation from the protocol used by the source device or bus to the 
5 protocol used by the target device or bus. 

Main memory 112 is a memory in which application programs are stored and 
from which microprocessor 110 primarily executes, A suitable main memory 112 
comprises DRAM (Dynamic Random Access Memory), and preferably a plurality of 
banks of SDRAM (Synchronous DRAM). 

10 

FIG. 8 -DSL Device 

FIG. 8 is a block diagram illustrating DSL device 120 comprised in the computer 
102 of FIG. 6. As noted above, although in the preferred embodiment the DSL device is 
comprised as a computer add-in card, DSL device 120 may take any of various types of 

15 forms including being comprised in a television system, Internet apphance, or dedicated 
device, among others systems. As shown, the DSL device or add-in card 120 comprises a 
first port 160 for coupling to an expansion bus of the computer system, preferably a PCI 
expansion bus port as shown. DSL device 120 also includes at least one subscriber line 
port 170 for coupling to the digital subscriber line 12. DSL device 120 may include any of 

20 various hardware elements for performing the communication operations of the present 
invention. For example, in one embodiment, the DSL communication device includes one 
or more programmable processing units which implement instructions from a memory. 
For example, the DSL communication device may include a programmable digital signal 
processor (DSP) 152, a general purpose processor, or other processors that execute 

25 instructions from a memory 156 to implement the commimication operations of the 

present invention. Alternatively, or in addition, the DSL communication device 120 
includes one or more application specific integrated circuits (ASICs) 154 and 158 or 
programmable logic devices such as FPGAs etc. that implement a portion or all of the 
present invention. In short, the communication system and method of the present 
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invention may be implemented in any of various types of ways including programmable 
devices such as processing units, CPUs, DSPs, microcontrollers, etc., dedicated hardware 
such as ASICs, programmable logic devices such as FPGAs, or combinations of the 
above. 

5 

FIG. 9 - FIG. 12 - Method for Determining Transmission Characteristics 

FIG. 9 is a flowchart for determining a transmit spectrum for use in 
communicating data on a communications channel according to one embodiment of the 
present invention. This method may be used in commmiicating data on the 

10 communications channel when the communications channel is subject to interference 
from one or more other communications channels. The communications channel of 
interest and one or more of the other communications channels carry a particular type of 
service, such as xDSL, ISDN, Tl, or spread-spectrum, for example. The first steps in this 
method comprise determining a channel transfer function of the communications channel 

15 210 and an amount of self interference 220 into the communications channel from the 
other communications channels that carry the same type of service. In step 230, the 
transfer function and the amount of self interference are examined, and in step 240 a 
transmit spectrum for the channel is determined based on the examining. 

Determining the channel transfer function in step 210 of FIG. 9 may be done by 
20 directly measuring it. For example, a transmitter on one end of the communications 
channel, such as at CO 12 (in FIG. 1) or in well-logging instrumentation 18 (in FIG. lA), 
may be directed to send a signal or a series of signals with predetermined intensities as a 
function of frequency, with which the a receiver at the other end of the channel may 
measure the attenuation, and perhaps also the phase shift, as a function of frequency. The 
25 measurement may be extended to determine nonlinear responses of the channel by 
repeating the measurement with varying source strengths. Alternately, the channel 
characteristics may be determined in advance of the communication and stored, for 
example in a database at the CO or in a memory on a DSL card. Determining the channel 
transfer function could then entail receiving it from the CO, the local memory storage, or 
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other storage locations. In the case where the invention is used in a subscriber line 
system, receiving the channel transfer function from the CO is particularly useful since 
the CO may rapidly look up pre-stored information on the particular physical line being 
used for the commxmications channel. 

5 Similarly, the amount of self interference may be determined in step 220 of FIG. 9 

by receiving it or, if transmitter/receiver pairs are accessible on the other same-service 
channels, by measurmg it. Determining the amount of self interference in step 220 may 
comprise determining a total self interference power, or a power distribution, or a 
coupling coefficient from the other same-service channels into the channel of interest, or 

10 a coupUng coefficient with frequency dependence (such as a self-interference transfer 
function) from the other same-service channels into the channel of interest, or a 
combination of these characteristics, among others. The amount of self interference may 
also include a characterization of the self interference in terms of self-NEXT and self- 
FEXT interference. In a preferred embodiment of the invention, step 220 includes 

15 determining a self-NEXT transfer function and a self-FEXT transfer function from the 

other same-service channels into the channel of interest. 

FIG. 10 - FIG. 12 show various embodiments 240a-c of step 240, in which the 
transmit spectrum is determined. 

In one embodiment of the invention, as shown in steps 241 and 245 of FIG. 10, 
20 determining the transmit spectrum in step 240 comprises determining an EQPSD transmit 

spectrum if the amount of self interference is substantially low or negligible. An EQPSD 
transmit spectrum is a transmit spectrum in which EQPSD signaling is used on at least 
one portion of the available spectrum of communication frequencies. 

The method may also include steps 242 and 247, in which an EQPSD/FDS 
25 transmit spectrum is found if the amount of self interference is substantially high or non- 
negligible. 

In general, an EQPSD/FDS transmit spectrum has a number of frequency regions 
in which EQPSD signaling is used, and a number of frequency regions in which FDS 
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signaling is used. The locations of these regions in the available spectrum of 
communication jfrequencies and the transmission power as a function of frequency are 
preferably determined so that the data transmission rate on the chaimel is substantially 
maximized. An EQPSD/FDS transmit spectrum preferably includes at least one portion 
5 using FDS signaling and one portion usmg FDS signaling, but in a degenerate case, the 
maximization may be achieved by using only EQPSD or only FDS signaling. An 
EQPSD/FDS transmit spectrum is thus a transmit spectrum in which the available 
spectrum of communication frequencies includes at least one portion using EQPSD 
signaling or FDS signaUng. 

10 Still further, the method may also include step 249, in which an 

EQPSD/FDS/MFDS transmit spectrum is found if the amount of self interference is 
substantially high or non-negligible and if the amount of self-FEXT interference is 
substantially high. In multi-line FDS (MFDS) signaling, different channels carrying the 
same or similar services are orthogonally separated to reduce crosstalk interference. 

15 hi general, an EQPSD/FDS/MFDS transmit spectrum has a number of frequency 

regions in which EQPSD signaling is used, a number of frequency regions in which FDS 
signaling is used, and a number of frequency regions in which MFDS signaling is used. 
Again, the locations of these regions in the available spectrum of communication 
frequencies and the transmission power as a function of frequency are preferably 

20 determined so that the data transmission rate on the channel is substantially maximized. 
An EQPSD/FDS/MFDS transmit spectrum preferably includes at least one portion using 
FDS signaling, one portion using FDS signaling, and one portion using MFDS signaling, 
but in a degenerate case, the maximization may be achieved by using only EQPSD or 
only FDS or only MFDS signaling. An EQPSD/FDS/MFDS transmit spectrum is thus a 

25 transmit spectrum in which the available spectrum of communication frequencies 

includes at least one portion using EQPSD signaling or FDS signaling or MFDS 
signaling. 

In another embodiment of the invention, as shown in FIG. 11, determining the 
transmit spectrum in step 240 comprises determining an EQPSD transmit spectrum (in 
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step 245) if the amount of self interference is substantially low or negligible (according to 
step 241), and determining an EQPSD/FDS/MFDS transmit spectrum (in step 249) if the 
amount of self interference is substantially high or non-negligible. 

In yet another embodiment of the invention, as shown in FIG. 12, determining the 
transmit spectrum in step 240 comprises determining an EQPSD/FDS/MFDS transmit 
spectrum (in step 249). For the degenerate cases in which only one or two of the EQPSD, 
FDS, and MFDS signaling techniques are needed for maximizing the data transmission 
rate, the EQPSD/FDS/MFDS transmit spectrum will reduce to the appropriate signaling 
techniques. 

In a preferred embodiment of the invention, the method further comprises a step 
of determining an amount of uncorrelated interference, such as additive Gaussian noise 
(AGN), into the communications channel. If one or more of the other communications 
channels carry a different type of service than the service on the communications channel, 
then the uncorrelated interference may include different-service interference (DSIN) from 
the other communications channels carrying the different service. Thus, the uncorrelated 
interference includes a total noise interference that preferably comprises AGN, DSIN, and 
other noise and interference whose spectral characteristics are not controlled by the user. 
Determining the transmit spectrum in step 240 is then performed in response to the 
amount of uncorrelated interference. 

FIG, 13 - FIG. 14 - Method for Transmitting Data 

FIG. 13 is a flowchart of a method for transmitting data according to one 
embodiment of the present invention. This method may be used in communicating data 
on a communications channel when the communications channel is subject to 
interference from one or more other communications channels. The other 
communications channels may be located proximate to the communications channel, for 
example, in the case of multiple subscriber hnes in a binder group of a telephone cable, or 
in the case of multiple radio transmission systems with closely located transmitters or 
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overlapping coverage regions. The commxmications channel of interest carries a particular 
type of service, such as xDSL, ISDN, Tl, or spread-spectrum, for example. The method 
comprises the steps of determining a channel transfer function of the communications 
channel in step 310, initiating a data transfer on the communications channel in step 320, 
5 and transferring the data on the communications channel using the transmit spectrum in 
step 330. As shown in FIG, 14, the step 320 of initiating the transfer comprises 
determining interference characteristics of the interfering communications chaimels in 
step 322, and determining a transmit spectrum in response to the channel transfer function 
and the interference characteristics in step 324. 

10 As discussed earlier, step 310 of determining the channel transfer function may 

comprise measuring the channel transfer function, receiving the channel transfer function, 
or determining the channel transfer function through other means. The channel transfer 
function may be determined at power-up of a transmission system, or at regular intervals 
in time, or in response to temperature changes, or at other appropriate times. 

15 The transmit spectrum is preferably determined in step 324 to substantially 

maximize the data transmission rate for the communications channel, so that the 
maximum information may be communicated per unit time on the communications 
chaimel in light of the various sources of noise and interference. The transmit spectrum is 
also preferably determined in such a manner that the communications channel has equal 

20 upstream and downstream capacities, and that the transmit spectrum is spectrally 
compatible (that is, determined with regard to spectral compatibility) with the one or 
more other communications channels. 

In one embodiment of the present invention, the transmit spectrum is preferably 
determined so that it satisfies a predetermined average power constraint for the 
25 commimications channel. Note that if the channel capacities depend on the transmit 
spectra of other lines carrying the same service, for example in the case of self 
interference, then the water filling technique may be carried out as described in reference 
[16]. If the channel capacity depends on channel noise and/or different-service 
interference, then the classical water-filling technique is used, as described in [14], The 
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transmit spectrum is preferably determined dynamically so that it may be optimized in 
response to changing interference conditions or a changing channel transfer function. 

In another embodiment of the invention, the transmit spectrum is determined so 
that it satisfies both a predetermined average powder constraint and a predetermined peak 
power constraint for the communications channel, and may be determined using a peak 
constrained water-filling technique. Note that if the channel capacities depend on the 
transmit spectra of other hnes carrying the same service, for example in the case of self 
interference, then the peak constrained water filling technique may be carried out as 
described in section 4.8.3 (which presents a modification of the technique discussed in 
[16]). If the channel capacity depends on channel noise and/or different-service 
interference, then the peak constrained water-filling technique is used, as described in 
section 4.8.2. The transmit spectrum is preferably determined dynamically so that it may 
be optimized in response to changing interference conditions or a changing channel 
transfer function. 

In another embodiment of the invention, the transmit spectrum is determined so 
that it satisfies only a predetermined peak power constraint for the communications 
channel, and may be determined using a peak constraiaed water-filling technique. 

Dynamical determination of transmit masks 

In a preferred embodiment of the invention, steps 322 and 324 of determining the 
interference characteristics and of determining the transmit spectrum are performed more 
than once so that the transmit spectrum is modified appropriately as the interference 
characteristics change in time. These steps 322 and 324 may be performed each time a 
data transfer is initiated. Or, if step 330 of transferring data occurs repeatedly at regular or 
irregular intervals in time, then steps 322 and 324 of determining the interference 
characteristics and of determining the transmit spectrum are preferably performed prior to 
each occurrence of transferring data in step 330. In one embodiment of the invention, a 
new transfer function or a new set of interference characteristics may be determined 
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during a data transfer and used to calculate a new transmit spectrum. The new transmit 
spectrum may then be used in a subsequent portion of the data transfer. These measures 
of dynamically determining the transfer function enhance the data transfer by allowing the 
transfer function to adapt as the characteristics of the communications channel change in 
5 time. 

Orthogonality for upstream/downstream separation and multi-line separation 

hi one embodiment of the present invention, the transmit spectrum is determined 
so that it specifies a pair of complementary spectra: one for transmission in each of the 

10 two directions on the communications channel. These two spectra may be called the 
"upstream transmit spectrum" and the "downstream transmit spectrum." For example, in 
the case where the channel provides communication between home 10 and CO 14, the 
transmit spectrum used in transmission from home 10 may be designated the upstream 
transmit spectrum, while the transmit spectrum used in transmission from CO 14 may be 

15 designated the downstream transmit spectrum. Similarly, in other cases, such as a well- 
logging or a multiple-radio-transmitter embodiment, "upstream" and "downstream" 
indicate opposite directions of transmission as desired. 

The upstream and downstream transmit spectra may include one or more regions 
of the spectrum that use FDS signaling. In these regions, the upstream and downstream 

20 transmit spectra are orthogonal with respect to each other, hi a preferred embodiment of 
the present invention, this duplexing orthogonality is achieved by choosing two non- 
overlapping frequency subregions in the FDS region, using one of the subregions for 
upstream signaling, and using the other subregion for downstream signaling. More 
generally, the FDS region may be constructed by choosing two non-overlapping sets of 

25 frequency subregions in the FDS region, using one of the sets for upstream signaling, and 
using the other set for downstream signaling. In another embodiment of the invention the 
duplexing orthogonality is achieved by using code division signaling (CDS) to separate 
the upstream and downstream signals in the "FDS" region. In this embodiment, one 
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access code is used in upstream signaling, and a second, orthogonal, access code is used 
in downstream signaling. 

It is noted that there is an additional benefit to these transmit spectra with one or 
more regions of FDS signaling: as would be appreciated by one skilled in the art of 
communications electronics, using regions of orthogonally separated upstream and 
downstream signaling may reduce the overhead of echo cancellation. 

The method indicated in FIG. 13 and FIG. 14 may be used in communicating data 
on a communications channel in a situation where one or more of the other 
communications channels carries the same type of service as the communications 
channel. Under such a condition, step 322 of determining interference characteristics 
preferably includes determining an amount of self interference into the communications 
from the other same-service communications channels. Step 324 of determining the 
transmit spectrum may then include examining the channel transfer function and the 
amount of self interference. The transmit spectrum is then preferably determined in step 
324 in response to the channel transfer function and the amomt of self interference. 

In another embodiment of the present invention, the transmit spectrum is 
determined so that it specifies a number M of complementary spectra: one for 
transmission on each of M channels in a subset of the one or more of the other 
communications channels that carry the same type of service. These M transmit spectra 
may include one or more regions of the spectrum that use MFDS signaUng. In these 
regions, the M transmit spectra are orthogonal with respect to each other. In one 
embodiment of the present invention, this multi-line orthogonality is achieved by 
choosing M non-overlapping frequency subregions in the MFDS region, and using one of 
the subregions for transmission on each of the MUnes. More generally, the MFDS region 
may be constructed by choosing M non-overlapping sets of frequency subregions in the 
MFDS region, and using one of the sets for transmission on each of the M channels. In 
another embodiment of the invention, the multi-line orthogonality is achieved by using 
multi-line code division signaling (multi-line CDS) in the "MFDS" region. In this 
embodiment, different orthogonal access codes are used on each of the M channels. 
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In another embodiment of the present invention, the transmit spectrum is 
determined so that it specifies a number M' (>M) of complementary spectra: one for 
transmission on each of M channels in the subset of same-service channels, and additional 
spectra to provide orthogonal duplex separation on one or more of the M channels. These 
5 M' transmit spectra may include one or more regions of the spectrum that use FDS 
signaling as well as one or more regions of the spectrum that use MFDS signaling. 

In a preferred embodiment of the invention, determining the amount of self 
interference comprises determining (1) a self-NEXT transfer function and (2) a self- 
FEXT transfer function that describe the coupling from near-end and far-end transmitters, 

10 respectively, on the other same-service communications channels. In this preferred 

embodiment determining the interference characteristics in step 322 further comprises 
determining an amount of uncorrected interference arising jfrom factors such as additive 
Gaussian noise (AGN) and crosstalk from one or more different-service channels, which 
carry a type of service different than the service on the channel of interest, among the one 

15 or more other channels. The transmit spectrum is then determined in response to the 

chaimel transfer function, an average power constraint or requirement for the channel, the 
self-NEXT and the self-FEXT transfer functions, and the amount of mcorrelated 
interference. In regions of the communications spectrum where the self interference is 
substantially low, the transmit spectrum is determined to be an EQPSD transmit 

20 spectrum. In regions where the self-NEXT interference is substantially high and the self- 
FEXT interference is not substantially high, the transmit spectrum is determined to be an 
FDS transmit spectrum. And in regions where the self-FEXT interference is substantially 
high, the transmit spectrum is determined to be an MFDS transmit spectrum. Some 
specific examples of techniques for determining regions of EQPSD, FDS, and MFDS 

25 signaling are presented below. 

In other embodiments of the invention, determining the interference 
characteristics in step 322 includes determining some but not all of the self-NEXT and 
the self-FEXT transfer functions, and the amount of uncorrelated interference and the 
average power constraint or requirement for the channel may or may not be determined. 
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The transmit spectrum is then determined in response to the channel transfer function and 
the determined quantities, and is preferably optimized in response to these quantities. 

In one embodiment of the present invention, the transmit spectrum is determined 
in response to one or more characteristics of the communications channel and the sources 
5 of noise or crosstalk. Determining these characteristics comprises steps such as 
determining the channel transfer function, determining the self-NEXT transfer function, 
and determining the self-FEXT transfer function, hi one preferred embodiment, the 
transmit spectrum is determined in response to the power-transfer characteristics of the 
communications channel, so determining these characteristics preferably comprises 
10 determining only the squared modulus of the mathematical transfer functions. Thus, in 

this preferred embodiment, determining the channel transfer function means determining 

the function H{f)^\H^{f^ ^ determining the self-NEXT transfer function means 
determining the function X{f)^\H^{f^ , and determining the self-FEXT transfer 
function means determining the function F{f) = \Hp{f^ An this preferred embodiment, 
15 the phases of the transfer functions H^{f), H^{f), and Hp{f) may or may not be 

determined in addition to their squared modulii. hi the case where a distinction is to be 
made between the various functions for different hnes, the subscript / is used to indicate 
the different lines (as in H^f), ^Xf), and with channel number /=1 being the 

channel for which the transmit spectrum is being determined. 

20 Another characteristic of the communications channel and the sources of noise or 

crosstalk is the signal to noise ratio (?/(/). In the case where a distinction is to be made 
between different lines, G^(/) mdicates, at a frequency / the ratio of the signal 
(specifically, the signal power spectral density at f) in channel number / to the noise 
(specifically, the noise spectral density at f) in channel number 1. Here, channel number 

25 i=\ is the channel for which the transmit spectrum is being determined, and channel 

number / (for / > 1) is another channel that carries the same type of service as channel 
number 1, and which may provide interference into channel number 1. 
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Method for Determining Transmission Characteristics with Frequency Binning 

Another embodiment of the present invention comprises a method for determining 
a transmit spectrum for use in communicating data on a communications chaimel, 
5 preferably by determining signaling techniques in one or more frequency bins in the 
available frequency band of the communications channel. This method is outlined in the 
flowchart of FIG. 15. This method may be used in communicating data on a 
communications channel w^hen the communications channel is subject to interference 
from one or more other communications channels, some of vs^hich carry the same type of 
10 service as the communications channel of interest. Additionally, some of the other 
communications chaimels may carry different types of service than the communications 
channel of interest. 

The first steps in this method comprise determining a channel transfer ftmction of 
the communications channel 410. An amount of self interference 420 into the 

15 communications chaimel from the other communications channels carrying the same type 
of service is determined in step 420. An additional amount of uncorrelated interference is 
preferably determined in step 425. In step 430, the transfer fimction and the amount of 
self interference are examined, preferably along v^th the amount of uncorrelated 
interference. In step 440 a transmit spectrum for the chaimel is determined based on the 

20 examining. 

In a preferred embodiment of the method, the transmit spectrum is determined in 
step 440 so that different signaling techniques may be used in different frequency ranges 
in the communications band. These frequency ranges, or frequency bins, are non- 
overlapping ranges of the frequency spectrum, preferably with uniform frequency widths, 
25 and preferably chosen so that they cover the communications band. In other embodiments 

of the present invention, the frequency bins have non-uniform widths or do not cover the 
entire communications band. 
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In this embodiment of the invention, the transmit spectrum operates to specify an 
amovmt of transmission power used in each frequency bin for at least one direction of 
communication on at least one communications channel. The amount of transmission 
power in each bin is preferably determined by a water-filling technique or a peak 
constrained water-filling technique. 

In one embodiment of the invention, in a given frequency bin the transmit 
spectrum specifies EQPSD signaling if the amoimt of self interference is substantially 
low in that bin, and FDS signaling if the amount of self interference is substantially high 
in that bin. 

In one embodiment of the invention, in a given frequency bin the transmit 
spectrum specifies MFDS signaling if the amount of self-FEXT interference is 
substantially high in that bin. Otherwise, the transmit spectrum specifies EQPSD 
signaling if the amoimt of self-NEXT interference is substantially low in that bin, and 
FDS signaling if the amount of self-NEXT interference is substantially high in that bin. 

Under certain conditions, this method of the present invention may determine a 
transmit spectrum that includes one or more regions of neighboring bins using FDS 
signaling. In one embodiment of the present invention, the step of determining a transmit 
spectrum comprises determining a discrete FDS transmit spectrum in such regions of 
neighboring FDS bins. In the discrete FDS transmit spectrum, each bin has two 
subregions; one is used for transmission in the upstream direction, and one for 
transmission in the downstream direction. In another embodiment of the present 
invention, the step of determining a transmit spectrum comprises determining a 
contiguous FDS transmit spectrum in such regions of neighboring FDS bins. In the 
contiguous FDS transmit spectrum, the neighboring frequency bins are grouped into two 
sets of neighboring bins, one of the sets is used for transmission in the upstream direction, 
and the other set is used for transmission in the downstream direction. In one embodiment 
of the invention, the two sets of neighboring bins are chosen so that the contiguous FDS 
transmit spectrum provides equal upstream and downstream signaling capacities. 
Alternatively, the two sets of neighboring bins may be chosen so that the contiguous FDS 
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transmit spectrum provides equal upstream and downstream average power. In a preferred 
embodiment, the two sets of neighboring bins are chosen so that the contiguous FDS 
transmit spectrum provides equal upstream and downstream signaling capacities and 
equal upstream and downstream average powers. 

5 Similarly, under certain conditions, this method of the present invention may 

determine a transmit spectrum that includes one or more regions of neighboring bins 
using MFDS signaling, hi one embodiment of the present invention, the step of 
determining a transmit spectrum comprises determhiing a discrete MFDS transmit 
spectrum in such regions of neighboring MFDS bins. In the discrete MFDS transmit 

10 spectrum, each bin has M subregions. Each of the M subregions is used for bi-directional 
transmission on one of the M same-service channels. In another embodiment of the 
present invention, the step of determining a transmit spectrum comprises determining a 
contiguous MFDS transmit spectrum in such regions of neighboring MFDS bins. In the 
contiguous MFDS transmit spectrum, the neighboring bins are grouped into M sets of 

15 neighboring bins. Each of the M sets of frequency bins is used for bi-directional 
transmission on one of the M same-service channels. The M sets of neighboring bins are 
preferably chosen so that the contiguous MFDS transmit spectrum provides equal 
signaling capacities on the M channels. 

If the charmel transfer function and the interference characteristics are 
20 substantially monotonic in frequency, the determination of which frequency bins use a 
particular type of signaling may be simpHfied by determining frequency values at which 
the different types of interference become substantially large or substantially small. Thus, 
in one embodiment of the present invention, determining the transmit spectrum in step 
440 includes one or more steps of identifying "transition bins" that mark the endpoints (in 
25 the frequency spectrum) of different types of signaling techniques. These transitions bins 
may be rapidly identified by searching for bins in which certain characteristic quantities 
meet particular predetermined criteria. These searches, which are preferably implemented 
as binary searches, may be carried out in the step 430 of examining the channel transfer 
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fiinction and interference. The following list is a sample of transition bins that may be 
identified. 

Me. for bins with center frequencies < or < the center firequency of Me, EQPSD 
signaling is used. 

5 Mf\ for bins with center frequencies > or > the center fi-equency of Mp, FDS signaling 

is used. 

Me2F' for bins with center frequencies < or < the center frequency of Me2f, EQPSD 
signaling is used, and FDS signaling is used in higher-frequency bins. In other 
words, Me2f indicates a transition from EQPSD signaling to FDS signaling. 

10 Me2mfds' indicates a transition from EQPSD signaling to MFDS signaling. 

Mmfds2fds' indicates a transition from MFDS signaling to FDS signaling. 

Me2mfds' indicates a transition from EQPSD signaling to FDS signaling. 

Similarly, transition frequencies may be defined for particular frequencies that 
15 mark transitions from one form of signaling to another. For example, represents a 
transition frequency where EQPSD is used in a region with frequency less ihsxifE2Fr> and, 
and FDS signaling is used in a region with frequency less than^^F- 
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4 New, Optimized Signaling Techniques 



The proposed techniques combine a number of ideas into one signaling system that optimizes its 
performance given many different possible combinations of interferers. These ideas include: 

1 . Given expressions for the crosstalk from other services (DSIN-NEXT and DSIN-FEXT) into 
an xDSL channel and channel noise (AGN), our scheme computes the optimal distribution of 
power across frequency that maximizes the capacity (see Section 4.4), This distribution uses 
the same transmit spectrum (EQPSD signaling) in both upstream and downstream directions, 

2. Given expressions for the self-NEXT and self-FEXT crosstalk in an xDSL channel along 
with interference from other services (DSIN-NEXT and DSIN-FEXT) and channel noise 
(AGN), our scheme computes the optimal distribution of power across frequency that max- 
imizes the capacity. This distribution involves equal PSD (EQPSD) signaling in frequency 
bands with low self-interference, orthogonal signaling (FDS) in frequency bands where self- 
NEXT dominates other interference sources (Section 4,5), and orthogonal signaling (multi- 
line FDS introduced in Section 43) in frequency bands where self-FEXT is high (Section 
4.6), 

3. Given different channel, noise, and interference characteristics between lines, our scheme 
chooses the optimal signaling strategy (EQPSD, FDS or multi-line FDS) in each frequency 
bin (see Section 4,7) to maximize the channel capacity. 

4. Given an additional peak-power constraint in frequency, our scheme computes the optimal 
transmit spectra that maximize the capacity and choose the optimal joint signaling strategy 
(EQPSD, FDS and multi-line FDS) for a given channel, noise and interference characteristics 
(see Sections 4,8 and 4.9)."^ 

5. We present optimal and near-optimal signaling strategies in case of non-monotonic channel, 
self-NEXT and self-FEXT transfer ftmctions (see Section 4.10 on bridged taps). 



We will present the above ideas in the following sections in the context of a generic xDSL line 
carrying symmetric-data rate services like HDSL2, "GDSL", and "VDSL2" services. Note that the 
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techniques developed here can be applied to a more general communications channel with inter- 
ference characteristics characterized by self-interference and different-service interference models. 
Further, v^e can extend this work to apply to channels that support asymmetric data rates (different 
in each direction) (see Section 5), for e.g„ ADSL, and VDSL. We can follow a similar approach 
of binning in frequency and then analyzing the signaling strategy in each bin. In the asynmietrical 
data-rate case, the ratio of the average power between upstream and downstream directions needs 
to be known. 

We will present background material and our assumptions in Section 4.L In Section 4.2 we 
give details about the interference models and the simulation conditions. Section 4.3 looks at the 
various signaling schemes we will employ. We will present the optimal transmit spectrum using 
EQPSD signaling in Section 4.4 in the presence of only different-service interference and AGN, 
Sections 4.5 and 4,6 detail the new signaling strategies to obtain an optimal and/or suboptimal 
transmit spectrum in the presence of self-interference, different-service interference and AGN. 
Section 4,7 derives some results applicable when neighboring lines vary in channel, noise and 
interference characteristics. Sections 4.8, and 4,9 present optimal transmit spectra under additional 
peak-power constraint in frequency. We present optimal and near-optimal signaling schemes for 
non-monotonic channel, self-NEXT, and self-FEXT transfer functions in Section 4.10. We discuss 
Optimal signaling for asymmetrical data-rate channels in Section 5. Finally, Section 4.12 presents 
several new ideas, extending the results presented here. 

Note: All the transmit spectra are optimal (le,, yield the maximum possible bit rates or per- 
formance margins) given the assumptions in Section 4 J (see Sections 4,5.3, and 4,6.3 for 
additional assumptions) and that one of the specific joint signaling strategies is employed over the 
channel (see Sections 4.4, 4.5, and4.6). 

4.1 Assumptions, Notation, and Background 

We present background material and some of the standard assumptions made for simulations. 
These assumptions apply throughout the document unless noted otherwise, 

1. Channel noise can be modeled as additive Gaussian noise (AGN) [13], 
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2, Interference from other services (DSIN-NEXT and DSIN-FEXT) can be modeled as additive 
colored Gaussian noise [13]. 

3, We assume the channel can be characterized as a LTI (linear time invariant) system. We 
divide the transmission bandwidth B of the channel into narrow frequency bins of width 
W (Hz) each and we assume that the channel, noise and the crosstalk characteristics vary 
slowly enough with frequency that they can be approximated to be constant over each bin 
(For a given degree of approximation, the faster these characteristics vary, the more narrow 
the bins must be. By letting the number of bins -K' -> oo, we can approximate any frequency 
characteristic with arbitrary precision).^ We use the following notation for line i on the 
channel transfer function [10] 



\Hcif)? = <! 
self-NEXT transfer function [8] 

\HN{f)f = { 

and self-FEXT transfer function [9] 

\HF{f)f = < 



Hi,k if|/-M<f, 

0 otherwise, 



Xi,k if|/-M<f, 

0 otherwise, 



Fij, if|/-/fc|<f, 
0 otherwise. 



(1) 



(2) 



(3) 



Here are the center frequencies (see Figures 16 and 17) of the K subchannels (bins) with 
index k e{l,..., K}. We will employ these assumptions in Sections 4.5.4, 4.6.6, 4.6,8 and 
4.7.1. The DSIN-NEXT and DSIN-FEXT transfer functions are also assumed to vary slowly 
enough that they can be similariy approximated by a constant value in each frequency bin. 

Note that the concept of dividing a transfer function in frequency bins is very general and 
can include nonuniform bins of varying widths or all bins of arbitrary width (i.e., the bins 
need not be necessarily narrow). 



^ We divide the channel into narrow frequency bins (or subchannels) for our analysis only. This does not necessarily 
mean thai we need to use DMT as the modulation scheme. 
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4. Echo cancellation is good enough that we can ignore crosstalk from into iZf . We can relax 
this assumption in some cases where spectral regions employ FDS signaling (see Sections 
45, 4.6,4.7,4.9, and 4.10). 

5. All sources of DSIN-NEXT can be lumped into one PSD DSn{f) and all sources of DSIN- 
FEXT can be lumped into one PSD DSpif). 

6. All sources of self-NEXT can be added to form one overall self-NEXT source. 

7. All sources of self-FEXT can be added to form one overall self-FEXT source. 

8. Spectral optimization is done under the average input power constraint, i.e., the average input 
power is limited to P^ax (Watts) in each direction of transmission in the symmetric data-rate 
case . In the asymmetric data-rate case the average input power in upstream and downstream 
directions is limited to Pup (Watts)and Pdn (Watts) (see Section 5). 

9. The PSDs of the upstream and downstream transmission directions can be written using the 
notation introduced in Section 1.3.2. There are M interfering lines carrying the same service 
with index i G {1, . , . , M}. Denote the direction of transmission with index o € {u, d}, with 
u = upstream (to CO) and d = downstream (from CO). Denote the upstream and downstream 
PSDs on line i as: ' 

Sf{f): PSD on twisted pair i in upstream direction u, 
Sf if): ?SD on twisted pair i in downstream direction d. 

Further, we denote the upstream and downstream PSD on line i in a generic frequency bin 
(or subchannel) k as: 

5 •*(/): PSD on twisted paii i in upstream dircction u. 
Si{f): PSD on twisted pair i in downstream direction ± 

Note: When we refer to (/) we mean PSD on twisted pair i in a generic bin, demodulated 
to baseband (/ E [-W, W]) for ease of notation. When we refer to we mean PSD on 
a generic twisted pair in a generic bin, demodulated to baseband (/ e W]) for ease of 
notation. 
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10, We assume a monotone decreasing channel transfer function. However, in case the chan- 
nel transfer function is non-monotonic (e.g., in the case of bridged taps on the line), our 
optimization techniques can be applied in each individual bin independently. This scenario 
makes the power distribution problem more difficult however (see Section 4.10). 

11. We assume we desire equal channel capacities in upstream and downstream directions (ex- 
cept when the channel, noise, and interference characteristics between lines vary as in Sec- 
tion 4.7 and when we desire asymmetric data rates from both directions of transmission as 
in Section 5), 

4.2 Interference models and simulation conditions 

The interference models for different services have been obtained from Annex B of Tl, 4 13- 1995 
([9], the ADSL standard), with exceptions as in TlEl.4/97-237 [7]. The NEXT coupling model 
is 2-piece Unger model as in TlEl.4/95-127 [8]. BER was fixed at IQ-*^. Our optimal case re- 
sults were simulated using Discrete Multitone Technology (DMT) and were compared with that of 
MONET-PAM [1]. MONET-PAM uses Decision Feedback Equalizers (DFE) [20] in the receivers 
along with multi-level pulse amplitude modulation (PAM) scheme. The margin calculations for 
DFE margins were done per T1EL4/97-180R1 [11], Section 5.4.2,2.1,1. AGN of power -140 
dBm/Hz was assumed in both cases. MONET-PAM uses PAM with 3 bits/symbol and a baud rate 
of fbaud = 517.33 ksymbols/s. The actual upstream and downstream power spectra can be ob- 
tained from [1]. MONET-PAM spectra is linearly interpolated from 2x1552/3 Hz sampled data. 
The PAM line-transformer hpf comer, that is, the start frequency is assumed to be at 1 kHz. A 500 
Hz rectangular-rule integration is carried out to compute margins. The required DFE SNR margin 
for 10-"^ BER is 27.7 dB. _ 

To implement our optimal signaling scheme, we used DMT with start frequency 1 kHz and 
sampling frequency of 1 MHz. This gives us a bandwidth of 500 kHz and 250 carriers with carrier 
spacing of 2 kHz. No cyclic prefix (used to combat intersymbol interference (ISI)) was assumed, 
so the DMT symbol rate is same as the carrier spacing equal to 2 kHz. However, the scheme 
can easily be implemented by accounting for an appropriate cyclic prefix. The addition of cyclic 
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prefix lowers the symbol rate and hence lowers the transmission rate. No limit was imposed on the 
maximum number of bits per carrier (this is often done for simulations). Even with a 15 bits/carrier 
limit, the results should not change very much, as some of the test runs show, 

4.3 Signaling schemes 

The joint signaling techniques used in the overall optimized signaling schemes use one of the basic 
signaling schemes (see Figure 18) in different frequency bins depending on the crosstalk and noise 
combination in those bins, 

- Figure 18 illustrates the three signaling schemes: EQPSD, FDS and multi-line FDS (in the case 
of three lines).^ The Figure shows in frequency bin k the PSDs for each case (recall the notation 
introduced in Section 4.1, Item 9): 

• When crosstalk and noise are not significant in a frequency bin, EQPSD signaling is pre- 
ferred as it achieves higher bit rate than the other two orthogonal signaling schemes (see 
Section 4,5,5). In EQPSD signaling, the upstream and downstream PSDs are the same 
istif) = sfif)). 

r 

• When self-NEXT is high and self-FEXT is low in a bin and there are a large number of 
neighboring lines carrying the same service together, FDS signaling yields the highest bit 
rates by eliminating self-NEXT (we prove this in Section 4.5.5). In FDS signaling, each 
frequency bin is further divided into two halves, with all the upstream PSDs being same for 
all the lines and all the downstream PSDs being same for all the lines (5^(/) ± (/)), 
This type of orthogonal signaling completely eliminates self-NEXT but does not combat 
self-FEXT. _ 

• In frequency bins where self-FEXT is high, using FDS is not sufficient since self-FEXT still 
exists. In this case, doing multi-line FDS eliminates self-FEXT as well as self-NEXT and 
this achieves the highest bit rates when there are ony a few Unes and self-FEXT is high and 
donoinant over self-NEXT (we prove this in Section 4.6). In multi-line FDS signaling each 

^The signaling schemes EQPSD, FDS, and multi-line FDS work in general for M lines. 
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line gets a separate frequency slot (W/M for M lines carrying the same service) in each bin 
and the upstream and downstream PSDs for each line are the same (s^{f) ± s^{f) Vj ^ 
i, 0 e {% d}x 

We will see in future sections the exact relationships that allow us to determine which scheme 
is optimal given an interference and noise combination, 

4.4 Optimization: Interference from other services (DSIN-NEXT and DSIN- 
FEXT) - Solution: EQPSD signaHng 

In this scenario, each xDSL line experiences no self-interference (Figure 19 with neither self- 
NEXT nor self-FEXT), Hiere is only DSIN-NEXT and DSIN-FEXT from other neighboring ser- 
vices such as Tl, ADSL, HDSL, etc, in addition to AGN. The solution is well known, but will be 
useful later in the development of the subsequent novel (Sections 4.5, 46, 4.7, and 4.12) signaling 
schemes, 

44.1 Problem statement 

r 

Maximize the capacity of an xDSL line in the presence of AGN and interference (DSIN-NEXT 
and DSIN-FEXT) from other services under two constraints: 

1. The average xDSL input power in one direction of transmission must be limited to F^ax 
(Watts). 

2. Equal capacity in both directions (upstream and downstream) for xDSL. 

Do this by designing the distribution of energy over frequency (the transmit spectrum) of the xDSL 
transmission. 

44.2 Additional assumption 

We add the following assumption to the ones in Section 4.1 for this case: 
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12. Both directions (upstream and downstream) of transmission experience the same channel 
noise (AGN) and different service interference (DSIN-NEXT and DSIN-FEXT). 



4.4.3 Solution 

Consider a line (line 1) carrying xDSL service. Line 1 experiences interference from other neigh- 
boring services (DSIN-NEXT and DSIN-FEXT) and channel noise No{f) (AGN) but no self- 
NEXT or self-FEXT (see Figure 19). 

The DSIN-NEXT and DSIN-FEXT interference can be modeled as colored Gaussian noise for 
calculating capacity [13]. Recall that DSNif) is the PSD of the combined DSIN-NEXT and let 
DSrif) is the PSD of the combined DSIN-FEXT. Let 5"(/) and 5^(/) denote the PSDs of line 
1 upstream (u) direction and downstream (d) direction transmitted signals respectively. Further, 
let and denote the upstream and downstream direction capacities of line 1 respectively. Let 
Hcif) denote the channel transfer function of line 1. The twisted pair channel is treated as a 
Gaussian channel with colored Gaussian noise. In this case the channel capacity (in bps) is given 
by [14] 



TOO 

C" = sup / logs 



!ffc(/)r^(/) 



and 



No{f)+DS^{f) + DSF{f) 

1 + 



df (4) 



\Hc{f)\'SV) 



N„{f)) + DSrf{f) + DSF{f)\ 



df. (5) 



r 

sup riog, 

S^{f) 

The supremum is taken over all possible 5"(/) and S'^{f) satisfying 

5"(/)>0V/, 5^(/)>0V/, 

and the average power constraints for the two directions 

2/ S«(/)d/ and 2 / Sy)df < Pniax. (6) 

^0 Jq 

It is sufficient to find the optimal S"" (/) which gives since setting 5^(7) = 5"* (/) V/, "gives the 
capacity — as seen from (4) and (5). Thus, the optimal upstream and downstream channel 
capacities are equal (C^ = C^). 



52 



The optimal power distribution in this case is obtained by the classical "water-filling" tech- 
nique [14]. The optimal is given by 



5?pt(/) = 



0 otherwise; 

with A a Lagrange multiplier and E the spectral region where > 0, We vary the value 

of A such that SQ^^{f) satisfies with equality the average power constraint in (6), The equality is 
satisfied for a single value of A giving us a unique optimal PSD S^^. (/). Plugging the optimal PSD 
^opt(/) (4) yields the capacity under the average power constraint This procedure yields a 
unique optimal transmit spectrum 5'opt(/) [14]. 

Keynote: S^{f) = V/ — EQPSD signaling. 

Figure 20 gives a flowchart to obtain the optimal transmit spectrum using only EQPSD sig- 
naling in the presence of DSIN-NEXT, DSIN-FEXT and AGN. It uses the classic water-filling 
solution to obtain the transmit spectrum. The novelty is in applying this to xDSL scenario to 
achieve a dynamic transmit spectrum (different for each interference type). 

The channel capacities can be calculated separately for each direction of transmission in case 
of nonuniform interference between the two directions, i.e., when the additional assumption in 
Section 4.4.2 does not hold. The transmit spectra in general will be different (S'^if) ^ S'^{f)) for 
this case, but will still occupy the same bandwidth. 



4AA Examples 

In this Section, we present some examples for the HDSL2 service. An average input power (Pmax) 
of 20 dBm and a fixed bit rate of 1.552 Mbps was used for all simulations. The performance margin 
was measured in each simulation and the comparison with other static transmit spectra (obtained 
from static PSD masks) proposed is presented in Section 4.5.11. Figure 21 shows the optimal 
upstream and downstream transmit spectrum for HDSL2 in the presence of DSIN-NEXT from 49 
HDSL interferers and AGN (-140 dBm/Hz). Note the deep null in the transmit spectrum fi*om 
approximately 80 to 255 kHz. This results from "water-filling" — the peak of the first main lobe 
of HDSL lies in the vicinity of 80 to 255 kHz. 
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Figure 22 shows the optimal upstream and downstream transmit spectrum for HDSL2 in the 
presence of DSIN-NEXT from 25 Tl interferers and AGN (-140 dBm/Hz). 

The optimal transmit spectra for the two cases are significantly different, evidence of the 
fact that the optimal transmit spectra will change depending on the nature of the interfer- 
ence. 

Summary: Recall the discussion on static PSD masks of Section 3.1. We have seen that the 
optimal transmit spectrum varies significantly with the interference combination. The water-filling 
solution yields a unique transmit spectrum for each interference combination [14]. The optimal 
transmit spectrum adapts to mininiize the effect of the interference combination. The optimal 
transmit spectra for upstream and downstream direction are the same (EQPSD signaling) and thus, 
employ the same average power in each direction. 

4.5 Optimization: Interference from other services (DSIN-NEXT and DSIN- 
FEXT) plus self-interference (self-NEXT and low self-FEXT) - Solution: 
EQPSD and FDS signaling 

In this scenario each xDSL line experiences self-interference (high self-NEXT and low self-FEXT) 
in addition to AGN and DSDs[-NEXT and DSIN-FEXT from other services (see Figure 3) in a 
generic xDSL service. This is the case of interest for HDSL2 service. 

4.5.1 Self-NEXT and self-FEXT rejection using orthogonal signaling 

As we saw in Section 3.2, orthogonal signaling can completely reject self-NEXT In addition, FDS 
gives better spectral compatibilit}^with other services than other orthogonal schemes like CDS or 
TDS (see Section 4.5.12 for a proof). Therefore, we choose to use the FDS scheme for orthogonal 
signaling. Recall the FDS signaling tradeoff: FDS eliminates self-NEXT and therefore increases 
system capacity; however, FDS also reduces the bandwidth available to each transmitter/receiver 
pair and therefore decreases system capacity. 

To eliminate self-FEXT using orthogonal signaling, we would force each upstream transmitter 
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If to be orthogonal to all other transmitters Tf, j i. Using multi-line FDS, we would separate 
each 37 into different frequency bands. Unfortunately, this would reduce the bandwidth available 
to each transmitter to 1/M the overall channel bandwidth. In a typical implementation of HDSL2, 
M will lie between 1 and 49; hence orthogonal signaling (multi-line FDS) for eliminating self- 
FEXT is worth the decrease in capacity only when self-FEXT is very high. We will show later 
in Section 4.6 that multi-line FDS gives gains in capacity when there are only a few number of 
interfering lines carrying the same service (M = 2 to 4) 

In this scenario, we assume self-NEXT dominates self-FEXT and self-FEXT is not very 
high (see Figure 17 and [8]), so we will design a system here with only self-NEXT suppression 
capability. However, self-FEXT still factors into our design in an important way. This is a new, 
non-trivial extension of the work of [3]. 

4^.2 Problem statement 

Maximize the capacity of an xDSL line in the presence of AGN, interference (DSIN-NEXT and 
DSIN-FEXT) from other services, and self-NEXT and self-FEXT under two constraints: 

1. The average xDSL input power in each direction of transmission must be limited to F^ax 
(Watts), and 

2. Equal capacity in both directions (upstream and downstream) for xDSL. 

Do this by designing the distribution of energy over frequency (the transmit spectrum) of the up- 
stream and downstream xDSL transmissions. 



4.53 Additional assumptions'^ 

We add the following assumptions to the ones in Section 4. 1 for this case: 

12. The level of self-FEXT is low enough in all bins that it is not necessary to use orthogonal 
signaling between different transmitter/receiver pairs operating in the same direction (see 
Section 4.5.1), 
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13. All the M lines considered are assumed to have the same channel and noise characteristics 
and face the same interference combination (interference combination refers to combmation 
of different interfering services) in both transmission directions (upstream and downstream). 
We will develop some results in Section 4 J for when this does not hold true. Thus, we 
assume that the upstream PSDs of all lines are the same {S'^if)) and the downstream PSDs 
of all lines are the same {S^{f)y That is, 

S\f) = 5f(/), iG{l,...,M}. (8) 

14. The coupling transfer fxmctions of NEXT and FEXT interference are synmietrical between 
neighboring services. For example, each line has the same self-NEXT transfer function 
i?jv (/) and self-FEXT transfer function Hp{f) for computing coupling of interference power 
with any other line. However, we develop some results in Section 4.7 when there are different 
NEXT and FEXT coupling transfer functions between lines. 

4.5.4 Signaling scheme 

Since the level of self-NEXT will vary with frequency (recall Figure 17), it is clear that in high 
self-NEXT regions of the spectrum, orthogonal signaling (FDS, for example) might be of use in 
order to reject self-NEXT. However, in low self-NEXT regions, the loss of transmission bandwidth 
of FDS may outweigh any gain in capacity due to self-NEXT rejection. Therefore, we would like 
our signaling scheme to be general enough to encompass both FDS signaling, EQPSD signaling, 
and the spectrum of choices in between. Our approach is related to that of [3], 

Key to our scheme is that the upstream and downstream transmissions use different transmit 
spectra. All upstream (to CO) transmitters transmit with the spectrum 5'*(/) All downstream 
(from CO) transmitters transmit with the spectrum S\f) Implicit in our scheme is the fact that 
in this case, self-NEXT dominates self-FEXT and self-FEXT is small If not, it would not be wise 
to constrain all 7^ to the same transmit PSD, 

Our goal is to maximize the upstream capacity (C^) and the downstream capacity (C^) 
given an average total power constraint of P^nax and the equal capacity constraint = C^. 
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Consider the case of two lines with the same service. Line 1 upstream capacity is and line 2 
downstream capacity is C"^. Under the Gaussian channel assumption, we can write these capacities 
(in bps) as 



sup / logs 

and 



1+ mf)\'s-if) 



Noif)+DSr,{f) + DSj.{f) + \H:,{f)\' S<^{f) + \Hp{f)\^ S-{f)\ 



df, 
(9) 



sup 



1^. \Hc{f)?SM) 



Noif) + DSM + DSpif) + \HM\' S-{f) + \Hpif)\' S'^{f)\ 



df. 
(10) 



The supremum is taken over all possible S^{f) and S^{f) satisfying 

5"(/)>0V/, 5V)>0V/, 
and the average power constraints for the two directions 

2 'r < P^«, and 2 T S''{f)df < P^^, (11) 

JQ JQ 

We can solve for the capacities (7" and using "water-filling" if we impose the restriction 
of EQPSD, that is = 5^{/) V/. However, this gives low capacities. Therefore, we employ 
FDS (5"(/) orthogonal to 5^(/)) in spectral regions where self-NEXT is laige enough to limit our 
capacity and EQPSD in the remaining spectrum. This gives much improved performance. 

To ease our analysis, we divide the channel into several equal bandwidth subchannels (bins) 
(see Figure 16) and continue our design and analysis on one frequency bin k assuming the subchan- 
nel frequency responses (l)-(3). Recall that Figure 17 shows that the channel and self-interference 
frequency responses are smooth and justifies our assuming them flat over narrow subchannels. For 
ease of notation, in this Section set 



H^Hi^k, X = Xi^k, F^Fi^^ in(l)-(3), 
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(12) 



and 



N = Noih) + DS^ih) + DSpih), (13) 

the noise PSD in bin k. Note that N consists of both AGNplus any interference (DSIN-NEXT and 
DSIN-FEXT) from other services. Let s^[f) denote the PSD in bin fc of line 1 upstream direction 
and s^{f) denote the PSD in bin A; of line 2 downstream direction (recall the notation introduced 
in Section 4.1, Item 9). The corresponding capacities of the subchannel k are denoted by c"* and 

We desire a signaling scheme that includes FDS, EQPSD and all combinations in between in 
each frequency bin. Therefore we divide each bin in hal? and define the upstream and downstream 
transmit spectra as follows (see Figure 23): 



a 



2Pm 

w 



and 



Af) = 



a 



2Rn 

w 



if|/|<f, 
iff <|/|<P^, 

otherwise 



if|/l<f, 

iff <i/r< w, 

otherwise. 



(14) 



(15) 



Here is the average power over frequency range [0,T1^] in bin fc and 0.5 < a < 1. When 
a = 0.5, 5"(/) = s^if) V/ e [0,W] (EQPSD signaling); when a = 1, 5«(/) and s^if) are 
disjoint (FDS signaling). These two extreme transmit spectra along with other possible spectra (for 
different values of a) are illustrated in Figure 23. The PSDs 5"(/) and s'^if) are "symmetrical" or 
power complementary to each other. This ensures that the upstream and downstream capacities are 
equal (c^* = c^). The factor a controls the power distribution in the bin, and W is the bandwidth of 
the bin. 

Next, we show that given this setup, the optimal signaling strategy uses only FDS or EQPSD 
in each subchannel. 



^The power split-up in a bin does not necessarily have to be 50% to the left side of the bin and 50% to the right 
side of the bin as shown in Figure 23. In general any 50% - 50% power-complementary split-up between opposite 
direction bins wiH work. 
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4^.5 Solution: One frequency bin 



If we define the achievable rate as 



RA{s^{f),s'if))= r log. 



1 + 



N + s^{f)X + s^{f)F 



df, 



(16) 



then 



c'' = ,f^^^,RA{s''if),sV)) and c'=m^^^RA{s'(f),s-{f)). (17) 

Due to the power complementarity of s"(/) and s'*(/), the channel capacities c" and are equal. 
Therefore, we will only consider the upstream capacity c" expression. Further, we will use Ra for 
s''(/)) in the remainder of this Section. Substituting for the PSDs from (14) and (15) 
into (16) and using (17) we get the following expression for the upstream capacity 



W 

— max < logo 
2 0.5<a<l I ^2 



1 + 



W 



AT , {l-a)2PmX , a2P^F 



1 + 



il-a)2P„H 
W 



■i* T ^ -r ^ 



(18) 



Let G = ^ denote the SNR in the bin. Then, we can rewrite (18) as 



W ( 
max — < 

0.5<a<l 2 I 



1 + 



1 + (1 - a)GX + aGF 



+ log2 



1 + 



(1 - a)GH 



1 + aGX + (1 - a)GF 



1) 



(19) 



Note from (17) and (19) that the expression after the max in (19) is the achievable rate Ra^ 
Differentiating the achievable rate,(i?>i) expression in (19) with respect to a gives us 

dRA W {\ 1 + (1 - a)GX + aGF 



da 



W fr l + (l-a) 
21n2 \[l + (l-a)GX 



X 



t)GX + aGF + aGH 
GH{1 + (1 - a)GX + aGF) - aGHj-GX + GF) 



{l + {l-a)GX + aGFy 

1 + aGX + (1 - a)GF 
1 + aGX + (1 - a)GF + (1 - a)GH^ 
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+ 



-GHjl + aGX + (1 - a)GF) - (1 - a)GH{GX - GF) ] 1 

(l + aGX + (l-Q;)GF)2 Jj ^^"^ 

= G{2a - 1) [2(X - F) + G(X2 - F^) - H(l + GF)] L, (21) 

with L > 0 Va 6 (0, 1]. Setting the derivative to zero gives us the single stationary point a = 0.5. 
The achievable rate Ra is monotonic in the interval a € (0.5, 1] (see Figure 24). If the value 
a = 0.5 corresponds to a maximum, then it is optimal to perform EQPSD signaling in this bin. If 
the value a = 0.5 corresponds to a minimum, then the maximum is achieved by the value a = 1, 
meaning it is optimal to perform FDS signaling in this bin. No other values of a are an optimal 
option. See Figure 25. 

The quantity a = 0.5 corresponds to a maximum of Ra (EQPSD) if and only if ^ < 0 
Va € (0.5, 1]. For all Q! € (0.5, 1], the quantity (2a - 1) is positive and ^ is negative if and only 
if (see (21)) 

2{X - F) + G(X2 _ _ jj^^ ^ Qp-^ ^ Q 

This implies that 

G(X2 -F^-HF)<H- 2{X - F). 
Thus, the achievable rate Ra is maximum at a = 0.5 (EQPSD) 



lfX^-F^-HF<0 andG> ^,-^jf_-J) 



(22) 



or 



ifX^-f--HF>0 ^iG<§f^^^^. (23) 

In a similar fashion a = 0.5 corresponds to a minimum of Ra if and only if ^ > 0 Va € 
(0.5, 1]. This impUes that a = 1 corresponds to a maximum ofRA (FDS) since there is only one 
stationary point in the interval oTe [0.5, 1] (see Figure 24). For all a G (0.5, 1], ^ is positive if 
and only if 

2{X - F) + G(X2 - f2) - H{1 + GF) > 0. 

This implies that 

G(X2 -F^-HF)>H- 2{X - F). 
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Thus, the achievable rate Ra is maximum at Of = 1 (FDS) 

if X=_f^_Hf<o and G<|^^f^ (24) 

or 

if X' - - gf > 0 and G > ^~J^:Pp . (25) 

Thus, we can deteraiine whether the value a = 0.5 maximizes or minimizes the achievable 
rate by evaluating the above inequalities. If a = 0.5 corresponds to a maximum of Ra, then we 
achieve capacity c« by doing EQPSD signaling. If a = 0.5 corresponds to a minimum of Ra, 
then we achieve capacity c" by doing FDS signaling. This can be summed in test conditions to 
determine the signahng nature (FDS or EQPSD) in a given bin. Using (22) and (24) we can write 

Jf X^-F^-HF <0 then 

^_2P. ^Q^^^ if-2(Z-F) 

NW pgg X^-F^-HF' ^^^^ 

Also, using (23) and (25) we can write 

If - - ifF > 0 then 

• 2P^ ^QP^P ff-2(X-F) 

NW p^g X'^-F'^-HF' ^"^'^ 

Thus, we can write the upstream capacity c" in a frequency bin A; as 



c« = < 



(28) 



Note: Its always optimal to do either FDS or EQPSD signaling; that is, a = 0.5 or 1 only, 
FDS signaling scheme is a subset of the more general orthogonal signaling concept. However, 
of all orthogonal signaling schemes, FDS signaling gives the best results in terms of spectral 
compatibility under an average power constraint and hence is used here (see proof in Section 
4.5.12). In the case of a peak power constraint in frequency, other orthogonal schemes, such as 
CDS, could be more appropriate (see Section 4.12.6). 
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4^.6 Solution: All frequency bins 



We saw in Section 4.5.5 how to determine the optimal signaling scheme (FDS or EQPSD) in one 
frequency bin for the upstream and downstream directions. In this Section we will apply the test 
conditions in (26) and (27) to all the frequency bins to determine the overall optimal signaling 
scheme. Further, using "water-filling" (tiiis comprises of the classical water-filling solution [14] 
and an optimization technique to compute capacity in the presence of self-interference [16]) opti- 
mize the power distribution over the bins given the average input power (Pmax). 

We divide the channel into K narrow subchannels of bandwidth W (Hz) each (see Figure 
16). For each subchannel k, we compute the respective channel transfer function (Hcifk), self- 
NEXT (Hr,{fk)X self-FEXT (HF{fk)X DSIN-NEXT (DS^ifk)). DSIN-FEXT (DSpih)) and 
AGN {No{fk)). Then, by applying (26) and (27) to each bin in the generic xDSL scenario (with 
the usual monotonicity assumptions as outiined in Section 4.1),^ we can divide the frequency axis 
(K bins) into 3 major regions: 

1. The right side of (26) < 0 for bins [1,Me]* These bins employ EQPSD signaling (since 
power in every bin is > 0). 

2. The right side of (27) < 0 for bins [Mp, K]. These bins employ FDS signaling (since power 
in every bin is > 0) and < Mp. 

3. The signaling scheme switches from EQPSD to FDS signaling at some bin Me2f, which lies 
in the range of bins (M^, Mp). 

Figure 26 illustrates the situation of the 3 bins Me, Mp and Me2F^ hi the next Section we 

develop an algorithm to find the^timal bin Me2f and the optimal power distribution. 

^When the channel transfer function is non-monotonic (as in the case of bridged taps) a bin-by-bin approach may 
be required to achieve the optimal power distribution (see Section 4. 10). 
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43 J Algorithm for optimizing the overall transmit spectrum 

To find the optimal EQPSD to FDS switch-over bin Me2f and the optimal power distribution over 
all bins: 

L Set up equispaced frequency bins of width W (Hz) over the transmission bandwidth B of 
the channel The bins should be narrow enough for the assumptions (1H3) of Section 4.1 
to hold. 

2. Estimate the interference (DSIN-NEXT, DSIN-FEXT, self-NEXT and self-FEXT) and noise 
(AGN) PSDs. Lump the corresponding interference PSDs together into one PSD. 

3. Compute the bins Me and Mp using (26) and (27) as outiined in Section 4.5.6, 

4. Choose an initial estimate of Me2f is a great start). 

5. Choose an initial distribution of how much proportion of the total power (Pmax) should go in 
the spectrum to the left of Me2f and how much should go to the right. Denote these powers 
by Pb and Pp ^ F^ax - Pe respectively. 

6. Use water-filling to distribute these powers {Pp and Pp) optimally over frequency [14, 16] 
with EQPSD signaling in bins [1, Me2f] and FDS signaling in bins [Me2f + ^,KY Compute 
the subchannel capacity c" in each bin using (28). Calculate the channel capacity by 
sununing all subchannel capacities. 

7. Re-estimate the powers Pe and Pp. 

8. Repeat steps 6 to 7 for a range of powers Pp and Pp in search of the maximum channel 
capacity C". This search is guaranteed to converge [3]. 

9. Re-estimate the optimal EQPSD to FDS switch-over bin Me2F. 

10. Repeat steps 5 to 9 for a range of bin values for Me2f. 

1 1. Choose die bin number which yields the highest channel capacity as the true optimal bin 
Me2f after which the signaling switches from EQPSD to FDS. 
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Notes: 



1 . Standard minimization/maximization routines (like jfem in the software package MATLAB) 
can be used to search for the optimal powers Pe and Pp. 

2. We can use fast algorithms like the Golden Section Search [19] to find the optimal bin Me2f- 
This routine tries to bracket the minimum/maximum of the objective function (in this case 
capacity) using four function-evaluation points. We start with a triplet (p, q, r) that brackets 
the minimum/maximum. We evaluate the function at a new point x e (q, r) and compare this 
value with that at the two extremities to form a new bracketing triplet (p, q, x) or {q, x, r) for 

- the minimum/maximum point. We repeat this bracketing procedure till the distance between 
the outer points is tolerably small 

4.5.8 Fast, suboptimal solution for the EQPSD to FDS switch-over bin 

In the estimation of the optimal bin Me2f we have observed in practice that Me2f ^ Me, typi- 
cally within 1 or 2 bins especially when self-interference dominates the total crosstalk (see Section 
4.5. 1 1). In the case of low AGN and different-service interference the suboptimal solution is a sub- 
stantially optimized solution. Thus, with significantly less computational effort than the algorithm 
described in Section 4.5.7, a near-optimal solution can be obtained. Even if a search is mounted 
for Me2f, we suggest that the search should start at Me (and move to the right). 

Algorithm to implement the suboptimal solution: 

1. Perform Steps 1 and 2 of the algorithm of Section 4.5.7. 

2. Compute the bin Me using-<26) as outlined in Section 4,5.6. 

3. Set the EQPSD to FDS switch-over bin Me2f equal to Me^ 

4. Obtain the optimal power distribution and the channel capacity C by performing Steps 5 
through 8 of the algorithm in Section 4.5.7. 
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4^.9 Flow of the scheme 

Consider a line carrying an xDSL service satisfying the assumptions of Sections 4.1 and 4.5.3. 
Lines carrying the same xDSL service and different xDSL services interfere with the line under 
consideration. We wish to find the optimal transmit spectrum for the xDSL line under consideration 
(see problem statement in Section 4.5.2). 

1. Determine the self-NEXT and self-FEXT levels due to other xDSL lines, bin by bin. These 
can be determined either through: 

(a) a worst-case bound of their levels determined by how many lines of that xDSL service 
could be at what proximity to the xDSL line of interest; or 

(b) an adaptive estimation (training) procedure run when the modem "turns on." In this 
process the CO will evaluate the actual number of active self-mterfering xDSL lines 
and the proximity of those lines with the line of interest. 

2. Determine DSIN-NEXT and DSIN-FEXT levels, bin by bin. These can be determined either 
through: 

(a) a worst-case' bound of their levels determined by how many lines of which kinds of 
service could be at what proximity to the xDSL line of interest; or 

(b) an adaptive estimation (training) procedure run when the modem "turns on". In this 
procedure no signal ti-ansmission is done but we only measure the interference level 
on the xDSL line at the receiver. Finally, the combined DSIN-NEXT and DSIN-FEXT 
can be estimated by subtracting the self-interference level from the level measured at 
the receiver. _^ 

3. an adaptive estimation (training) procedure run when the modem "turns on". 

4. Optimize the spectrum of transmission usmg the algorithms of Section 4.5.7 or 4.5,8. 

5. Transmit and receive data. 
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6. Optional: Periodically update noise and crosstalk estimates and ti^nsmit spectrum from 
Steps 1-3. 

Figure 27 illusti-ates a flowchart showing the steps for the optimal and the suboptimal solution. 
4^.10 Grouping of bins and wider subchannels 

The optimal and near-optimal solutions of Sections 4.5.7 and 4.5.8 divide die channel into narrow 
subchannels (bins) and employ die assumptions as discussed in Sections 4.1 and 4.5.3. In the case 
of self-interference, the resulting optimal tiransmit spectiimi uses FDS and is "discrete" (a "line 
spectirum"). Such a transmit spectrum is easily implemented via a DMT modulation scheme, but 
is not easy to implement with other modulation schemes like PAM, multi-level PAM, or QAM 
[20]. In addition, the DMT scheme can introduce high latency which may be a problem in some 
applications. Thus, one may want to use other low-latency modulation schemes. In such a scenario, 
we can combine or group FDS bins to form wider subchannels and then employ other broadband 
modulation schemes. This may result in different performance margins but we believe that the 
change in margins would not be significant. An alternative broadband modulation scheme like 
multi-level PAM or QAM would use a decision feedback equalizer (DFE) [20] at the receiver to 
compensate for the channel attenuation characteristic (see Section 4.12.4 for further discussion). 

Figure 28 shows one possible way of grouping die bins. The left-hand-side figures show the 
optimal upstream and downstream "discrete" transmit spectra S''(/) and 5'*(/) as obtained by 
the algorithm of Section 4.5.7. The right-hand-side figures show die same optimal transmit spectra 
after appropriate grouping of bins resulting in "contiguous" tiransmit spectra. While grouping, only 
die bins employing FDS signaling are grouped togetiier and flie leftmost bins employing EQPSD 
signaling are retained as they areT In diis particular case, we have grouped die bins such that die 
upstream and downstream capacities are equal (C" = C**). The upstream tiransmit spectiimi is 
completely "contiguous" while the downstream spectrum is "contiguous" except for one "hole" as 
shown in Figure 28. 

Note: This is not the only way that the bins can be grouped. The bins can be grouped in 
a variety of different ways giving many different optimal transmit spectra. Particular modula- 
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tion schemes and spectral compatibility with neighboring services may influence the way bins are 
grouped. Further, grouping of bins may lead to different input powers for opposite directions of 
transmission. 

We look at another possible way of grouping bins such that we achieve equal performance 
margins and equal upstream and downstream average powers. This could be a preferred grouping 
for symmetric data-rate services. 

Algorithm for "contiguous'^ optimal transmit spectra: Equal margins and equal average 
powers in both directions: 

1 . Solve for the optimal transmit spectrum iS" (/) according to the algorithms in Sections 4,5,7, 
4.5.8, or 4,6, where is the water-filling solution (refer to [14, 16] if the spectral region 
employs EQPSD or multi-line FDS signaling and to [16] if the spectral region employs FDS 
signaling) (see Sections 4.5 and 4.6). This gives a discrete transmit spectrum 5''(/). 

2. Denote the spectral region employing FDS signaling as E'fds and the spectral region em- 
ploying EQPSD signaling as -Eeqpsd* 

Obtain S^[f) from 5«(/) by synametry, i.e., S\f) = 5"(/) in EQPSD and multi-line FDS 
regions and S\f) 1 S\f) in FDS spectral regions. Merge S\f) and to form S{f) 
as 

S{f) = S^{f)^S\f) V/in^EQPSD, 
S{f) = 5«(/)u5^(/) V/inJE^FDS, (29) 
where U represents the union of the two transmit spectra. 

3. Estimate bins Mc € {M^2^,K], and Mq € {Mc,K]. Group the bins of 5(/) to obtain 
upstream and downstream transmit spectra as 

S{f) V/inEEQPSD, and 
-^optC/) = \ V/ in bins {Mc, Mc], (30) 

0 otherwise, 
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sUf) = { 



S{f) V/in^EQPSD, and 

V/ in bins (Me2f, Mc], and 
V/ in bins {Mg^K], 

0 otherwise. 



(31) 



4. Iterate previous step for various choices of Mc and Mq. The bin Mc is chosen such that we 
get equal performance margins in both directions of transmission and the bin Mq is chosen 
such that upstream and downstream directions have equal average powers. 



The resulting transmit spectra S^^^{f) and S^^^{f) are another manifestation of the grouping of 
bins and yield equal performance maigins (equal capacities) and equal average powers in both 
directions of transnndssion. 



4.5.11 Examples and results 

In this Section, we present some examples and results for the HDSL2 service. AGN of —140 
dBm/Hz was added to the interference combination in all simulations. Table 1 lists our simulation 
results performance margins and compares them with results from [1]. The simulations were 
done for the Carrier SeiVing Area (CSA) loop number 6, which is a 26 AWG, 9 kft line with no 
bridged taps. The columin "Our-PAM" refers to our implementation using T1EL4/97-180R1 [11] 
of the PAM scheme (MONET-PAM) suggested by the authors in [1] using their transmit spectra. 
We believe the shght differences in margins between MONET-PAM and "Our-PAM" exist due 
to slight differences in our channel, self-NEXT and self-FEXT models. The use of "Our-PAM" 
maxgins allows us a fair comparison of our optimal results with other proposed transmit spectra. 
The columns Up and Dn refer to the upstream and downstream performance maigins respectively. 
The column Optimal refers to the performance margins obtained using the optimal transmit spectra. 
The column Diff shows the difference between the performance margins for the optimal transmit 
spectrum and the MONET-PAM transmit spectrum (using "Our-PAM" margins). A full-duplex bit 
rate of 1.552 Mbps and a BER of lO"*^ was fixed in order to get the performance margins. The 
HDSL2 standards committee desires a high uncoded margin (preferably more than 6 dB), Table 1 
shows that we achieve very high uncoded margins far exceeding current schemes. 
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Table 1 : Uncoded performance margins (in dB) for CSA No, 6: MONET-PAM vs. Optimal. 



Crosstalk source 


xDSL service 


MONET-PAM 


"Our-PAM" 


Optimal 


Diff 


Up 


Dn 


Up 


Dn 


49 HDSL 


HDSL2 


9.38 


3.14 


10.05 


3.08 


18.75 


15.67 


39 self 


HDSL2 


10.3 


6.03 


11.18 


6.00 


18.39 


12.39 


25 Tl 


HDSL2 


19.8 


20.3 


14.23 


20.29 


21.54 


7.31 



Bit rate fixed at 1.552 Mbps. 

Diff = Difference between Optimal and worst-case "Our-PAM", 

Table 2 shows the difference between the optimal solution of the signaling scheme (using the 
optimal Me2f) and the fast approximate suboptimal solution (using Me2f = Me) for a variety of 
interfering lines. The column Diff (in dB) notes the difference in performance margins between the 
optimal scheme and the suboptimal scheme. Note that there is hardly any difference between the 
two when self-interference dominates the total crosstalk. This is a very significant result fi-om an 
implementation view point for it shows that near-optimal signaling can be obtained with very 
little computational effort. The optimal solution requires a somewhat complicated optimization 
over the bins starting firom Me and moving towards the right Our results clearly indicate that 
the near-optimal solution can give extremely attractive results with no search for the optimal bin. 
Further, this suggests that the optimal bin Me2f is closer to Me than and so one should search 
for it to the immediate right of Mg, 

An optimal upstream transmit spectrum in the case of self-interference is illustrated in Figure 
29, The Figure shows the optimal upstream transmit spectrum for HDSL2 service in the presence 
of self-NEXT and self-FEXT from 39 HDSL2 disturbers and AGN of -140 dBm/Hz, The down- 
stream transmit spectra for the HDSL2 service are synmietric with the upstream transmit spectra 
as discussed earlier. 

Figure 30 illustrates optimal "contiguous" transniit spectra for the same case of 39 self-NEXT 
and self-FEXT disturbers with AGN of -140 dBmuTHz, The "contiguous" transmit spectra were 
obtained by grouping the bins as outlined in Section 4.5.10 (C^ = C*^). The upstream and down- 
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Table 2: Uncoded performance margins (in dB) for CSA No. 6: Optimal vs. Suboptimal. 



Crosstalk source 


xDSL 
service 


Optimal 
scheme (dB) 


Me2F 


Fast, suboptimal 
scheme (dB) 


Mb 


Diff 


Iself 


HDSL2 


27.68 


11 


27.68 


10 


0 


10 self 


HDSL2 


21.94 


10 


21.94 


10 


0 


19 self 


HDSL2 


20.22 


8 


20.22 


8 


0 


29 self 


HDSL2 


19.13 


8 


19.13 


8 


0 


39 self 


HDSL2 


18.39 


9 


18.39 


9 


0 


lOself+lOHDSL 


HDSL2 


12.11 


60 


11.46 


19 


0.65 


lOself+lOTl 


HDSL2 


7.92 


27 


7.90 


23 


0.02 



Bit rate fixed at 1,552 Mbps. 

Diff = Difference between Optimal and suboptimal scheme, 

stream directions exhibit the same performance margins and use different powers. 

Figure 31 illustrates another set of optimal "contiguous" transmit spectra for the same case of 
39 self-NEXT and self-FEXT disturbers with AGN of -140 dBm/Hz. These "contiguous" transmit 
spectra were obtained by grouping the bins as outlined in the algorithm of Section 4,5. 10 such that 
now we have both equal performance margins (equal capacities) and equal average powers in both 
directions of transmission. 

4.5.12 Spectral compatibiUty 

When we optimize the capacity of an xDSL service in the presence of interferers, we must ensure 
that the optimized xDSL servicers not spectrally incompatible with other services. That is, the 
performance margins of other services must not significantly degrade due to the presence of that 
xDSL. Our optimal xDSL transmit spectra involve water-filling (after choosing the appropriate 
joint signaling strategy). To maximize xDSL capacity we distribute more power in regions of 
less interference and vice versa. This implies the services which interfere with xDSL see less 
interference in spectral regions where they have more power and vice versa* This suggests that the 
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spectral compatibility margins for other services in the presence of optimized xDSL PSD should 
be high. 

Table 3 lists our simulation results for HDSL2 service and compares them with results from [1]. 
The simulations were done for the CSA loop number 6 (26 AWG, 9 kft, no bridged taps) and CSA 
loop number 4 (26 AWG, bridged taps). The column "Our-PAM" refers to our implementation 
using T1E1.4/97-180R1 [1 1] of the PAM scheme (MONET-PAM) suggested by the authors in [1] 
using thek transmit spectra. We believe the sUght differences in margins between MONET-PAM 
and "Our-PAM" exist due to the differences in our channel, self-NEXT and self-FEXT models. 
The column Optimal lists the performance margins of tiie xDSL service under consideration using 
the -optimal transmit spectrum only when HDSL2 is a crosstalk source. The use of "Our-PAM" 
margins allows us a fair comparison of our optimal margins with the other proposed transmit 
spectira. From Table 3, we can clearly see that the optimal transmit spectrum has a high degree of 
spectral compatibility with the surrounding interfering lines. 

Our optimal results in case of self-NEXT and self-FEXT give rise to FDS signaUng, which 
has a peaky PSD in bins employing FDS. All ortiiogonal schemes tike FDS, TDS, and CDS give 
self-NEXT rejection and can transmit at the same bit rate. But, using FDS is better than CDS since 
tiiere is a gain in the perfprmance margin of tiie interfering line. We now prove that FDS signaling 
gives higher spectral compatibility margins than other orthogonal schemes like CDS. 

Theorem: Let the line under consideration be tiie signaling Une (with PSD 5 in a 
single bin) and tiie line that interferes with this line be tiie interfering line (witii PSD 
s'^if) and 5''(/) in a single bin). Then, using an FDS scheme instead of CDS scheme 
for tiie interfering line results in higher capacity for tiie signaling Une under an average 
power constraint and a Gaussian channel model. 

Proof: Consider, as usual tiie scenario of one single frequency bin of widtii W (Hz) as illustirated 
in Figure 32. In tiiis Figure, S is tiie transmit spectrum of tiie signaling Hne under consideration 
(for example Tl, HDSL, ADSL, etc.), y and Z represent tiie different service interference powers 
from a neighboring interfering Une (for example HDSL2) and N represents tiie lumped channel 
noise (AGN) and otiier different-service interference. There are two cases of interest: 
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Table 3: Spectral-compatibility margins: MONET-PAM vs. Optimal 







MONET-PAM Dn 


"Our-PAM" 


Optimal 


Crosstalk Src 


xDSL Srvc 


CSA6 


CSA4 


CSA6 


CSA4 


CSA6 


CSA4 


49 HDSL 


HDSL 


8.53 


8.09 


8.09 


7.78 






39 HDSL2 Up 


HDSL 


10.1 


10.9 


9.74 


10.53 


15.44 


15.60 


39 HDSL2 Dn 


HDSL 


8.28 


7.99 


7.74 


7.53 






39 HDSL 


EC ADSL 


8.43 


9.55 


7.84 


9.02 






39 HDSL2 


EC ADSL 


9.70 


11.7 


8.17 


10.00 


6.93 


9.10 


49 HDSL 


EC ADSL 


8.12 


9.24 


7.52 


8.7 






49HDSL2 


HDSL 






7.10 


6.91 


14.95 


15.12 



1: The interfering line uses a CDS signaling scheme. In this case the power in a single bin k 
(Pm) is uniformly distributed throughout the bin resulting in a flat PSD, i.e., 

= s'if) = a. 

We assume the subchannel frequency responses (1H3) and the notation introduced in (12) 
and (13). We assume here that the NEXT and FEXT coupling transfer functions between 
different service lines are the same as that for same-service lines. Thus, we can write the 
different service interference power in signaling line bin k as 

DSMiD + DSpif) = s'^iDX + s^DF 

= aX+oF. (32) 

We define Y and Z as 

Y = a{X-F) 

Z = 2aF (33) 

Using (33) we can write the interference power in (32) as 

DSNif) + DSF{f) = Y + Z. 
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Case 2: The interfering line uses an FDS signaling scheme. In this case the power in a single bin k 
{Pm) is distributed in only half the bin, resulting in a peaky PSD, i.e., 



and. 



2a, ifl/l<f, 

0, iff <\f\<W, 



0, if|/l<f, 



2o, iff <!/!<P^. 

We assume the subchannel frequency responses (l)-(3) and the notation introduced in (12) 
and (13). We assume here that the NEXT and FEXT coupling transfer functions between 
• different service lines are the same as that for same-service lines. Thus, we can write the 
different service interference power in signaling line bin k as 

DS^if) + DSfH) = s-{f)X + s'{f)F 



2aX, if|/|<f, 



2aF, iff <|/|<W. 
Using (33) we can write the interference power in (34) as 



(34) 



DSNif) + DSF{f)-^< 



2 ' 



2F + Z, if|/|<^ 
Z, iff<|/|<W. 



Getting back to the problem, we consider a single signaling line (line 1). We divide the signal- 
ing line channel into narrow subchannels (or bins) and we analyze a narrow subchannel k. We use 
the standard assumptions of Section 4.1. We can write the upstream subchannel capacity of bin k 
of the signaling line in Case 1 as 



1 + 



Y + Z + N 



and in Case 2 as 



1 + 



2Y + Z + N 

Compute the capacity differences in the two cases as 



1 + 



Z + N 



]}• 



(35) 



(36) 



D = c5'(Case 2) - c^(Case 1) = b 

73 



(37) 



Taking the partial derivative of D with respect to Y we get 



W~ln2 [{Y + Z + N){Y + Z + N + S) {2Y + Z + N){2Y + Z + N + 



Let 



U = Y+Z+N, 
V = Y + Z + N + S. 

Note that 17, V" > 0 and that we can rewrite the partial derivative of D with respect to Y as 



dD_W_ 
dY ~ ln2 



UV {U + Y){V + Y) 



W 
= —S 
In 2 



Y^ + {U + V)Y 



UV{U + Y){V + Y) 



>0. 



(38) 



Further, = 0. The slope of D with respect to Y is always positive and hence, c?(Case 2) - 

c5'(Case 1) is always increasing with Y, which implies that 

c5'(Case 2) - c^(Case 1) W > 0. 

When y < 0, i.e., when FEXT is higher than NEXT in a bin (F > X), we can redefine F and Z 

as 

Z = 2aX, and, Y = a{F- X). 
We can then follow the same analysis and show that the capacity c5^(Case 2) is greater than c?(Case 1). 

Thus, we have proven that FDS scheme rather than CDS scheme for interfering lines, results 
in higher capacities for signaling lines under an average power constraint. Q.E.D. 

Interestingly, the power-peaky FDS transmit spectra should be very compatible with the ADSL 
standard, since ADSL can balance how many bits it places in each of its DMT subchannels using 
a bit loading algorithm [17]. — 

Note: FDS beats CDS in terms of spectral compatibility margins only in the case of an average 
power constraint. In the case of a peak power constraint in frequency, we may not be able to use a 
power-peaky scheme like FDS in some spectral regions. Here, we may find that other orthogonal 
signaling schemes like CDS offer better spectral compatibility margins. 
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4.6 Optimization: Interference from other services (DSIN-NEXT and DSIN- 
FEXT) plus self-interference (self-NEXT and high self-FEXT) - Solu- 
tion: EQPSD, FDS and multi-line FDS signaling 

In this scenario we have self-interference (self-NEXT and high self-FEXT) in addition to AGN and 
DSIN-NEXT and DSIN-FEXT from other services (see Figure 3) in a generic xDSL service. This 
is the case of interest for "GDSL", '*VDSL2", and HDSL2 (with a small number of lines). 

4.6.1 Self-FEXT and self-NEXT rejection using multi-Une FDS 

To reject self-FEXT and self-NEXT, we use muW-Une FDS (see Section 4.3 and Figure 18). In 
multi-line FDS we separate each line by transmitting on each in different frequency bands. This 
reduces the transmission bandwidth to 1/M the total channel bandwidth, with M the number of 
lines carrying the service under consideration. Thus, multi-line FDS signaling can increase the 
capacity only when there are a few number of lines. 

We will design a system here that has both self-NEXT and self-FEXT rejection capability. 
Thus, this serves as the complete solution under the assumptions in Section 4.1 and the constraints 
of limited average input power (Pmax) and equal capacity in both directions. 

4.6.2 Problem statement 

Maximize the capacity of an xDSL line in the presence of AGN, interference (DSIN-NEXT and 
DSIN-FEXT) from other services, and self-NEXT and self-FEXT under two constraints: 

1. The average xDSL input power in each direction of transmission must be limited to Pmax 
(Watts), and 

2, Equal capacity in both directions (upstream and downstream) for xDSL. 

Do this by designing die distribution of energy over frequency (the transmit spectrum) of the up- 
stream and downstream xDSL transmissions. 
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4.6.3 Additional assumptions 



We add the following assumptions to the ones in Section 4.1: 

12. All the M lines carrying the xDSL service are assumed to have the same channel and noise 
characteristics and face the same interference combination in both transmission directions 
(upstream and downstream). Refer to Section 4.7 for results when this does not hold true. 

13, The coupling transfer functions of NEXT and FEXT interference are synmietrical between 
neighboring services. For example, each line has the same self-NEXT transfer function 
Hn if) and self-FEXT transfer function Hpif) for computing coupling of interference power 
with any other line. However, we develop some results in Section 4,7 when there are different 
NEXT and FEXT coupling transfer functions between lines. 

4.6.4 Signaling scheme 

The level of self-NEXT and self-FEXT varies over frequency (recall Figure 17), In regions of 
low self-NEXT and low self-FEXT, EQPSD signaling is the best choice. In spectral regions of 
high self-NEXT but low self-FEXT, orthogonal signaling scheme like FDS is preferred (due to its 
self-NEXT rejection, as we saw in Section 4.5). But, in regions of high self-FEXT, multi-line FDS 
signaling might be required for gaining capacity. 

Key to our scheme is that the upstream and downstream transmissions of each of the M lines 
use different transmit spectra, 

4.6.5 Solution using EQPSD and FDS signaling: All frequency bins 

First, we assume that self-FEXT is small and then, using EQPSD or FDS signaling in each bin, 
we find the solution for all frequency bins as outlined in Sections 4.5.4 — 4.5.8. Thus, we ob- 
tain the optimal (or suboptimal) EQPSD to FDS switch-over bin Me2F under the low self-FEXT 
assumption. 

Next, we relax the self-FEXT assumption and open the possibility of multi-line FDS. We 
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search each bin to see if we need to switch from EQPSD to multi-line FDS or FDS to multi-line 
FDS. This may not necessarily yield the optimal solution for the transmit spectrum given that we 
use a joint signaKng scheme comprising of the three signahng schemes (EQPSD, FDS and multi- 
line FDS). But, this analysis is tractable and gives significant gains in channel capacity and is 
presented next 

4.6.6 Switch to multi-line FDS: One frequency bin 

Consider the case of M lines with significant self-FEXT interference between them. We divide 
the channel into several equal bandwidth (W Hz) bins (see Figure 16) and perform our analysis on 
one frequency bin k assuming subchannel frequency responses (l)-(3). We employ the notation 
introduced in (12) and (13). Let 3?(/) denote the PSD in bin k of line 1 upstream direction and 
sfif) denote the PSD in bin k of line 1 downstream direction (recall the notation introduced in 
Section 4.1, Item 9). Let Pm be the average power over the frequency range [0, W]. 

Next, we determine when we need to switch to multi-line FDS in a given bin to completely 
reject self-FEXT: 

EQPSD to multi-line FDS: Figure 33 illustrates the two possible signaling schemes EQPSD 
and multi-line FDS in bin k of each line for the case of M ^ 3 lines . We will consider line 1 
for our capacity calculations. Line 1 upstream and downstream capacities for EQPSD signaling 
are denoted by c^' ^qpsd and cIeqpst> respectively. Similarly, line 1 upstream and downstream 
capacities for multi-line FDS signaling are denoted by c^^mfds and cf ^mfds respectively. Since the 
upstream and downstream transmit spectra of line 1 in bin A; for EQPSD and multi-line FDS are 
the same, we have: 

*^1,EQE§P ~ n,EQPSD3 ^IjMFDS ~ ^L,MFDS 

Thus, we will consider only the upstream capacities in our future discussion. 

Under the Gaussian channel assumption, we can define the EQPSD upstream capacity (in bps) 

as 
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C?,EQPSD=^l0g2 



1 + 



N + sfX + s^F 



(39) 



where 



= siif) = < 



^, if|/le[0,n . 
0 otherwise. 



Let G = ^ denote the SNR in the bin. Then we can rewrite cf ,eqpsd 



1 + 



GH 



2 + GX + GF\ 

Similarly, we can define the multi-line FDS upstream capacity (in bps) as 



(40) 



_W 



1 + 



N 



(41) 



where 



0, otherwise, 
and G — ^ is the SNR in the bin. Then we can rewrite Cj* mfds 



Define the difference between the two capacities as 



(42) 



n — /.» — 

^ — n,MFDS n,EQPSD- 



(43) 



We wish to determine when it is better to do multi-line FDS than EQPSD, i.e., when is the capacity 
Ci,MFDS greater than c" ^Qpgu. This means we need a condition for when D > 0. Substituting from 
(40) and (42) into (43) we get I? > 0 iff 



F> 



[2 +V{X + H)] + fGH) ^{2 + GX) 



Similarly, EQPSD is better (gives higher capacity) than multi-line FDS when D < 0, i.e., iff 

[2 -f G{X -h H)] - (l -}- f GF)^ (2 + GX) 
G {{l+^■GH)^ -1^ 
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(44) 



F< 



(45) 



We can combine (44) and (45) into one test condition that tells us the signaling scheme to use in a 
single frequency bin 



multi - Hne FDS [2 + G{X + H)] - (l + f Gff) ^ (2 + GX) 
EQPSD 



(46) 



FDS to multi-line FDS: Figure 34 illustrates the two possible signaling schemes FDS and multi- 
line FDS in bin k of each line for the case of M = 3 lines. We will consider line 1 for our capacity 
calculations. Line 1 upstream and downstream capacities for FDS signaling are denoted by p^^g 
and cf FDS respectively. Similarly, line 1 upstream and downstream capacities for multi-line FDS 
signaling are denoted by c^^mfds and cf^MFDs respectively. Since the upstream and downstream 
transmit spectra of line 1 in bin A: for EQPSD and multi-line FDS are the same, we have: 

<^1,FDS ~ '^IjFDS? *^1,MFDS " n,MFDS 

Thus, we will consider only the upstream capacities in our future discussion. Under the Gaussian 
channel assumption we can define the FDS upstream capacity (in bps) as 



W 



1 + 



Sim 



where 



slif) = { 



if|/l€[0,f], 



0, otherwise. 

Let G = ^ denote the SNR in the bin. Then we can rewrite c^,fds ^ 



_ W 



1 + 



GH 1 



1 + GF\ 

Similarly, we can define the multi-line FDS upstream capacity (in bps) as 



W 



1 + 



N 



where 



if|/|e[o,f], 

0, otherwise, 



(47) 



(48) 



(49) 
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and G = ^ is the SNR in the bin. Then we can rewrite c^,mfds as 



W 



Cl,MFDS — logs 

Define the difference between the two capacities as 



(50) 



D = cl 



MFDS ~ Ci^Dg. 



(51) 



We wish to find out when it is more appropriate to perform multi-line FDS than FDS, i,e., when the 
capacity dl^uFT>s is greater than c^^pog. For this, we need a condition for when D > 0. Substituting 
from (48) and (50) into (51) we get £>> 0 iff 

{l + GH)-(l + fGH)^ 
G (il + fGH)^ -Ij 
Similarly, FDS is better (gives higher capacity) than multi-line FDS when I> < 0, i.e., iff 

{l + GH)-(l + fGH) 



F> 



(52) 



_2_ 

I M 



F< ^ X N 

G Ul + f GH)" -Ij 



(53) 



We can combine (52) and (53) into one test condition which tells us the signaling scheme to use 

multi - Hne FDS + qjj^ ^GH)^ 
FDS g(^(i + Mgh)^-i) 

Thus, we can write the generic upstream capacity c" for bin k of line 1 as 
' Wlog, [l + ^;^^tfx+^) ] , if EQPSD, 



(54) 



flo&[l + 



if FDS, 



(55) 



, f log2 [l + ^] , if multi-line FDS. 



4.6.7 Switch to multi-line FDS: AU frequency bins 



We saw in the previous Section how to determine if we need to switch to multi-line FDS fi"om 
EQPSD or FDS in a given bin. We already have the optimal solution assuming EQPSD and FDS 
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signaling scheme (from Section 4.5), Now, we apply the conditions (46) and (54) to each bin 
A;. Interestingly, due to the assumed monotonicity of self-FEXT, self-NEXT and channel transfer 
function, we can divide the frequency axis (all K bins) into 4 major regions: 

1. Using test condition (46), we find that bins [1, Me2mfds] employ EQPSD signaling. 

2. Using test condition (46), we find that bins [Me2mfds + 1, Mmfds2Fds] employ multi-line 
FDS signaling. Note that Mmfds2Fds = ^E2F obtained from optimization procedure of 
Section 4.6.5. 

3. Using test condition (54), we find that bins [Mmfds2FDS + 1? -W'fi>S2Mfds] employ FDS sig- 
naling. 

4. Using test condition (54), we find that bins [Mfds2Mfds + employ multi-line FDS 
signaling. 

Figure 35 illustrates the 3 bins Me2mfds. ^mfds2Fds and Mfds2MFds the EQPSD, FDS 
and multi-line FDS regions. In practice we mainly see 2 scenarios: 

1. If Me2mfds < ^mfds2Fds thcu Mfds2Mfds = Mmfds2Fds. and we get only 2 distinct 
spectral regions as shown in Figure 36: 

(a) Bins [1, Me2mfds] employ EQPSD signaling. 

(b) Bins [Me2mfds + 1? -K'] employ multi-line FDS signaling. 

FDS signaling is not employed in this case, 

2. If Me2mfds = Mmfds2fi>s = Me2f then we get 3 distinct spectral regions as shown in 
Figure 37: 

(a) Bins [1, Mmfds2Fds] employ EQPSD signaling. 

(b) Bins [Mmfds2fds + h -^fds2Mfds] employ FDS signaling. 

(c) Bins [Mfds2Mfds + Ij ^] employ multi-line FDS signaling. 
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There is no switch to multi-line FDS signaling within the EQPSD signaling region (bins 
[1,Me2f]}- 



Note that the bin Mmfds2Fds = -^E2F is fixed from the optimization procedure from Section 4.6,5. 
4.6»8 Special case: Performance of 2 lines 

Often in practice we may have only two twisted pair lines carrying the same service and interfering 
with each other. It is important to derive the optimal transmit spectrum for such a scenario. In this 
Section we focus on this special case of only 2 lines. We will see that in this case it is optimal 
to perform either multi-line FDS or EQPSD signaling in each bin. In this scenario with arbitrary 
self-FEXT and self-NEXT we easily see that there is no need to perform FDS signaling (reject 
self-NEXT only) as multi-line FDS rejects both self-NEXT and self-FEXT while achieving the 
same capacity as FDS. Thus, we choose between EQPSD and multi-line FDS signaling schemes 
for each bin to achieve the optimal transmit spectrum. 

Let S'^ (/) and Sf{f) denote the upstream and downstream transmit spectra of line 1 and (/) 
and 51 (/) denote the upstream and downstream transmit spectra of line 2 respectively. Let the line 
1 upstream capacity hc^Ci and let the line 2 downstream capacity be Cg. Under the Gaussian 
channel assumption, we can write these capacities (in bps) as 




sup / 




l0g2 1 + 



\Hc{f)\'snf) 



df, (56) 



and 



= sup / 

S|(/).S?(/),5f (/) Jo 




N,{f) + DS^if) + DSp{f) + \HMf Snf) + \HFU)f Sfif). 




df. (57) 



The supremum is taken over all possible S^{f), 5J(/), Sf{f) and 5|(/) satisfying 



sm > 0, Sfif) > 0, S^if) > 0, S|(/) > 0,V/, 
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and the average power constraints for the two directions 

2 r Sfif)df < P^ax, and 2 T S^iDdf < P^^. (58) 

We employ multi-line FDS and 5f(/) orthogonal to 5|(/) and in spectral 

regions where the self-FEXT is lai^e enough and EQPSD in the remaining spectrum. This gives 
optimal performance. 

To ease our analysis, as usual, we divide the channel into several equal bandwidth subchannels 
(bins) (see Figure 16) and continue our design and analysis on one frequency bin k assuming 
subchannel frequency responses (l)-(3). We use notation introduced in (12) and (13). Let 3" (/) 
and (/) denote the PSDs in bin k of line 1 upstream and downstream directions and 53 (/) and 
s^{f) denote the PSDs in bin k of line 2 upstream and downstream directions. The corresponding 
capacities of the subchannel k are denoted by , cf , and c^. 

We desire a signaling scheme that can have multi-line FDS, EQPSD and all combinations in 
between in each frequency bin. Therefore we divide each bin in half^ and define the upstream and 
downstream transmit spectra as follows (see Figure 38): 



si{f) = siif) = { 

and 

4if) = siif) = \ 



«^ if|/l<f, 

iff <|/1<W, (59) 

0 otherwise, 



(l-c)^ ifl/|<f, 



iff<l/l<W, (60) 



0 otherwise. 

Here, Pm is the average power fiver frequency range [0, W] in bin k and 0.5 < a < 1. In this 
discussion we wiU only use the PSDs 5?(/) and When a = 0.5, = s^(/) V/ € [0, W] 
(EQPSD signaling); when = 1, 5^(/) and s^(/) are disjoint (multi-line FDS signaling). The 
PSDs s^{f) and s^{f) are "symmetrical" or power complementary to each other. This ensures the 



^The power split-up in a bin does not necessarily have to be 50% to the left side of the bin and 50% to the right side 
of the bin as shown in Pigure 38. In general any 50% - 50% power complementary split-up between different-line 
bins will work. 
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capacities of the two lines are equal {c^ = c^). The factor a controls the power distribution in the 
bin and W is the bandwidth of the bin. 

Next, we show that the optimal signaling strategy uses only multi4ine FDS or EQPSD in each 
subchannel 



The achievable rate for one frequency bin can be written as 



RA{snf).si{f),sUf))= / log: 

JO 



1 + 



N + si{f)X + 4{f)F\ 



df, 



(61) 



then 



c?=^max^i2^(5^:(/),4(/).s5(/)) and 4 = ,^^^,RA{si{f),sUf),s^f)). (62) 

Due to the power complementarity of Sj'(/) and ^{f), the channel capacities are equal (c" = c^). 
Therefore, we will only consider the upstream capacity c" expression. Further, we will use Ra for 
■^^(si (/), 52(/), §2 (/)) the remainder of this Section. Substituting for the PSDs from (59) and 
(60) into (61) and using (62) we get the following expression for the upstream capacity 



W 

— max { log, 

2 0.5<o<l ' 



1 + 



W 



iv -r ^y -r w 



1 + 2^ 



N + 



w ^ w 



(63) 



Let G = denote the SNR in the bin. Then, we can rewrite (63) as 



W I 

— max < 

2 0.5<a<l 1 



l0g2 



1 + 



aGH 



1 + (1 - a)GX + (1 - a)GF 



+ log2 



1 + 



(1 - a)GH 
1 + aGX + aGF 



(64) 



Using (62) and differentiating the achievable rate (Ra) expression in (64) with respect to a 
gives us 

dRA 



da 



= {2a - 1) [2{X + F) + G{X + Ff - H] L, 



(65) 



with L > 0 Va G (0, 1]. Setting the derivative to zero gives us the single stationary point a = 0.5. 
Thus, the achievable rate Ra is monotonic in the interval a G (0.5, 1] (see Figure 24). If the value 
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a = 0.5 corresponds to a maximum of Ra, then it is optimal to perform EQPSD signaling in this 
bin. If the value a = 0.5 corresponds to a rninimum of Ra^ then the maximum of Ra is acMeved 
by the value a = 1, meaning it is optimal to perform multi-line FDS signaling in this bin. No other 
values of q: are an optimal option (see Figure 39). 

The quantity a = 0.5 corresponds to a maximum o/Ra (EQPSD) if and only if ^ < 0 
Va 6 (0.5, 1]. For all a € (0.5, 1], the quantity (2a - 1) is positive and ^ is negative iff (see 
(65)) 

2{X + F) + G{X + Ff - H <Q. 

This implies that 

In a similar fashion a = 0.5 corresponds to a minimum of Ra if and only if ^ > 0 Va € 
(0.5, 1]. This implies that a = 1 corresponds to a maximum (multi-line FDS) since there is only 
one stationary point in the interval a G [0.5, 1] (see Figure 24). For all a € (0.5, 1], ^ is positive 
iff 

2{X + F)+ G{X + Ff-H>0. 

This implies that 

H-2{X + F) 

^ > ^x+Fy ■ ^^^^ 

The above statements can be sunmied in a test condition to determine the signaling nature 
(multi-line FDS or EQPSD) in a given bin. Using (66) and (67) we can write 

„ ^ 25^ ™ H-2(X + F) 

^ EQPSD (X + F)^ ■ 

Thus, we can write the upstisam capacity c" in a frequency bin k as 

Wlog, [1 + ^^"if^x+F) ] . if« = 0.5, 



(69) 



flog2[l + ^], ifa = l. 



Note: It is globally optimal to employ either multi-line FDS or EQPSD signaling; that is, 
a = 0.5 or 1, only in the case of 2 lines. 
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4»6.9 Flow of the scheme 



1. Perform steps 1-3 of Section 4.5.9. 

2. Compute bins Me2mfds. Mmfds2FDS and Mfds2MFds and employ signaling schemes in bins 
as described in Section 4.6.5. 

3. Transmit and receive data. 

4. Optional: Periodically update noise and crosstalk estimates and transmit spectrum from 
Steps 1-3 of Section 4.5.9, Repeat Step 2 from above. 

" Figure 40 gives a flowchart to obtain the optimal transmit spectrum using EQPSD, FDS, and 
multi-line FDS (MFDS) signaling in the presence of self-interference (self-NEXT and self-FEXT), 
DSIN-NEXT, DSIN-FEXT and AGN. 

4.6.10 Examples and results 

Optimal transmit spectra were used in all examples to compute performance margins and channel 
capacities. 

r 

HDSL2 service: Table 4 lists our simulation results performance margins and channel capacities 
using the EQPSD, FDS and multi-line FDS signaling schemes. 

Notes: 

1. SampUng frequency /, = 1000 kHz, Bin width W = 2 kHz and number of subchannels 
K = 250. Average input power of 20 dBm in each transmission direction, 

2. denotes the upstream capacity of line i using EQPSD and FDS signaling only and 
Cf (MFDS) denotes tiie upstream capacity of line i using EQPSD, FDS and multi-line 
FDS signaling schemes. All the rates are in Mbps. 

3. The column Margin lists the performance margin when the bit rate is fixed at 1.552 
Mbps. In each row in the top half the capacity is fixed at = 1.5520 and in the 
bottom half tiie capacity is fixed at QCMFDS) = 1.5520. 
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Table 4: Uncoded psTfommnce maigins (in dB)md channel edacities (in Mbps) using EQPSD, FDS and 
molti-line FDS for HDSL2 (CSA No. 6). 



XtalkSrc 


-^E2MFDS 


■WmFDS2FDS 


■^FDS2MFDS 




Cf(MFDS) 


Margin 


Diff 


1 HDSL2 


8 


11 


11 


1.5520 


2.3763 


27.682 


9.852 


1HDSL2 


0 


0 


0 


0.8027 


1.5520 


37.534 




2HDSL2 


9 


9 


30 


1.5520 


1.8293 


25.934 


4.543 


2HDSL2 


4 


4 


19 


1.1861 


1.5520 


30.477 




3HDSL2 


8 


8 


112 


1.5520 


1.6067 


24.910 


0.985 


3HDSL2 


7 


7 


100 


1.4792 


1.5520 


25.791 




4HDSL2 


8 


8 


246 


1.5520 


1.5520 


24.186 


0 



Diff = Difference between bottom half and top half of each row of Margin. 



4. The column Diff denotes the gain in performance margins between using EQPSD and 
FDS versus EQPSD, FDS and multi-line FDS signaling, i.e., the difference in margins 
between the bottom half and top half of each row. 

5. Each HDSL2 Une contributes NEXT and FEXT calculated using 2-piece Unger model 
[8], 

6. These runs were done with no different service (DS) interferers. The results would vary 
depending on the particular DS interferer(s) present. 

Conclusions: 

L Significant gains in margin for small number of lines. The gains decrease with increase 
in number of lines. 

2. There is no gain in margin using multi-line FDS for 5 or more lines (4 Crosstalk dis- 
turbers) for these line and interference models. 

"GDSL" service: Table 5 lists our simulation results performance margins and channel capacities 
using the EQPSD, FDS and multi-line FDS signaling schemes in the case of "GDSL". 

Notes: 
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Table 5: Uncoded performance margins (in dB) and channel capacities (in Mbps) using EQPSD, FDS and 
multi-line FDS for "GDSL" (3 kft line). 



XtalkSrc 


■^E2MFDS 


■^MFDS2FDS 


■WpDSSMFDS 




Cf(MFDS) 


Margin 


Diff 


1 GDSL 


505 


1253 


1253 


25.0046 


31.6188 


8.21 


8.49 


IGDSL 


245 


981 


981 


16.5141 


25.0007 


16.70 




2 GDSL 


952 


1214 


1214 


25.0007 


27.3923 


6.13 


2.91 


2 GDSL 


825 


1116 


1116 


22.0076 


25.0030 


9.04 




3 GDSL 


1186 


1212 


1212 


25.0004 


25.6686 


5.05 


0.75 


3 GDSL 


1145 


1186 


1186 


24.2172 


25.0008 


5.80 




4 GDSL 


1222 


1222 


2000 


25.0018 


25.0018 


437 


0 



Diff = Difference between bottoni half and top half of each row of Margin, 



1. Sampling frequency fs = 8000 kHz, Bin width W - 2 kHz and number of subchannels 
K = 2000. Average input power of 20 dBm in each transmission direction, 

2, Cf denotes the upstream capacity of line i using EQPSD and FDS signaling only and 

(MFDS) denotes the upstream capacity of line i using EQPSD, FDS and multi-line 
FDS signalirig schemes. All the rates are in Mbps, 

3. The column Margin lists the performance margin when the bit rate is fixed at 25 Mbps, 
In each row in the top half the capacity is fixed at Cf = 25 and in the bottom half the 
capacity is fixed at Cf (MFDS) = 25, 

4, The column Diff denotes the gain in performance margins between using EQPSD and 
FDS versus EQPSD, FDS and multi-line FDS signaling, i.e., the difference in margins 
between the bottom half and top half of each row, 

5, Each "GDSL" line contributes self-NEXT and self-FEXT calculated using 2-piece 
Unger model [8], In "GDSL" case the self-FEXT level is more dominant than self- 
NEXT. To model this we take only 1% of the self-NEXT power calculated using 2- 
piece Unger model in our simulations. 

6. These runs were done with no different service (DS) interferers. The results would vary 
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Table 6: Uncoded performance margins (in dB) and channel capacities (in Mbps) uang EQPSD, FDS and 
multi-line FDS for '*VDSL2" (3 kft line). 



XtalkSrc 


Me2MFDS 


MmFDS2FBS 


■fl^FDS2MFDS 


Cf 


Cf(MFDS) 


Margin 


Diff 


1 VDSL2 


58 


236 


236 


12.4011 


24.8234 


16.022 


18.913 


1 VDSL2 


8 


50 


50 


2.5552 


12.4001 


34.935 




2 VDSL2 


160 


219 


219 


12.4003 


18.8073 


14.074 


13.476 


2VDSL2 


46 


78 


78 


4.4478 


12.4036 


27.550 




3 VDSL2 


217 


217 


217 


12.4028 


15.6002 


12.985 


7.765 


3 VDSL2 


127 


127 


127 


7.3365 


12.4002 


20.750 




4 VDSL2 


219 


219 


553 


12.4016 


13.7787 


12J50 


3.275 


4 VDSL2 


179 


179 


359 


10.1474 


12.4012 


15.525 




5 VDSL2 


224 


224 


1014 


12.4014 


12.9039 


11.705 


1.005 


5 VDSL2 


211 


211 


878 


11.6945 


12.4014 


12.710 




6VDSL2 


231 


231 


1455 


12.4025 


12.5278 


11.280 


0.212 


6 VDSL2 


229 


229 


1412 


12.2521 


12.4018 


11.492 




7 VDSL2 


240 


240 


1880 


12.4004 


12.4049 


10.945 


0.007 


7 VDSL2 


240 


' 240 


1878 


12.3954 


12.4001 


10.952 





Diff - Difference between bottom half and top half of each row of Margin. 



depending on the particular DS interferer(s) present 
Conclusions: 

1 . Significant gains in margin for small number of lines. The gains decrease with increase 
in number of lines. 

2. There is no gain in margin using multi-line FDS for 5 or more lines (4 Crosstalk dis- 
turbers) for these line and interference models, 

"VDSL2" service: Table 6 lists our simulation results performance margins and channel capaci- 
ties using the EQPSD, FDS and multi-line FDS signaling schemes in the case of "VDSL2'\ 
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Notes: 



L Sampling frequency fs = 8000 kHz, Bin width W = 2 kHz and number of subchannels 
K = 2000. Average input power of 20 dBm in each transmission direction. 

2. Cf denotes the upstream capacity of line i using EQPSD and FDS signaling only and 

(MFDS) denotes the upstream capacity of line i using EQPSD, FDS and multi-line 
FDS signaling schemes. All the rates are in Mbps. 

3. The colunm Margin lists the performance margin when the bit rate is fixed at 12.4 
Mbps. hi each row in the top half the capacity is fixed at Cf = 12.4 and in the bottom 
half the capacity is fixed at Cf (MFDS) = 12.4. 

4. The column Diff denotes the gain in performance margins between using EQPSD and 
FDS versus EQPSD, FDS and multi-line FDS signaling, i.e., the difference in margins 
between the bottom half and top half of each row. 

5. Each VDSL2 line contributes self-NEXT and self-FEXT calculated using 2-piece Unger 
model [8]. In VDSL2 case self-NEXT and self-FEXT both are high but self-NEXT 
dominates self-FEXT. 
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Conclusions: 

1. Significant gains in margin for small number of lines. The gains decrease with increase 
innumber of lines. 

2. There is no gain in margin using multi-line FDS for 9 or more lines (8 crosstalk dis- 
turbers). These runs were done with no different service (DS) interferers. The results 
would vary depending on the particular DS interferer present. 

4.7 Joint signaling for lines differing in channel, noise and interference char- 
acteristics 

We have so far looked at a scenario where all the lines in a binder have the same channel charac- 
teristics and experience similar noise and interference characteristics in both directions of trans- 
mission. These assumptions made the signaling scheme solutions more tractable. We also need 
to look at a scenario between neighboring lines in binder groups where the channel characteristics 
vary (e.g., different length and different gauge lines) and we have different noise and interference 
characteristics between upstream and downstream transmission (e.g., asymmetrical services like 
ADSL and VDSL; different coupling transfer function in different directions). In this Section, we 
derive results for neighboring lines carrying the same service when they differ in channel, noise 
and interference characteristics. Specifically, we develop test conditions to determine the signaling 
nature in a given bin k. 

4J.1 Solution for 2 lines: EQPSD and FDS signaling 

Consider the case of 2 lines with "different channel, noise and interference characteristics. We again 
divide the channel into several equal bandwidth bins (see Figure 16) and continue our design and 
analysis on one frequency bin k assuming the subchannel frequency responses (l)-(3). For ease of 
notation in this Section, for line 1 we set 

Hr^Hi^k, X.^Xij,, F.^Fi^k asin(l)-(3), (70) 
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and let 

Ni = Noifk) + DSMifk) + DS^Uk). (71) 

be the lumped noise PSD in line 1 bin k. Further, let Pmi and Pm2 be the average powers over 
range [0, W] Hz in bin k of line 1 and 2 respectively. Let and sf(/)- denote the PSDs in bin 
k of line 1 upstream and downstream directions and §2 (/) ^iif) denote the PSDs in bin k of 
line 2 upstream and downstream directions (recall the notation introduced in Section 4.1, Item 9). 
The corresponding capacities of the subchannel k are denoted by , cf, and (4- 

We desire a signaling scheme that can have FDS, EQPSD and all combinations in between in 
a frequency bin. Therefore we divide each bin in half and define the upstream and downstream 
transmit spectra as follows (see Figure 41): 



a 



2P^i 
W 



if|/l<f, 



iff <\f\<W, 
0 otherwise, 

(l_a)2^ if|/|<f, 

iff<|/!<W, 
otherwise, 



(72) 



O!- 



W 



(73) 



a 



2Pm2 



if I/I < ^ 



w IJ i ^ 2 ' 

{l-ay-^ iff <\f\<W, 



(74) 



0 



otherwise, 



and 



sfif) = 



(1-0.)%- if|/|<f. 



a 



0 



2R 



ml 



W 



iff <|/|<W, 
otherwise, 



(75) 



where 0.5 < o; < 1. We assume that the upstream and downstream transmit spectra obey power 
complementarity, i,e, line 1 puts less power where line 2 puts more and vice versa. When oj = 0.5, 
^lif) = ^1 (/). ^2 (/) = 4if) V/ e [0, W] (EQPSD signaling); when a = 1, 4(/) and sUf) are 
disjoint (FDS signaling). The capacities of opposite directions are equal for each line: 

c? = cf and c^==c^. 
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The factor ot controls the power distribution in the bin, and W is the bandwidth of the bin. 

Next, we show that the optimal signaling strategy uses only FDS or EQPSD in each subchan- 
nel. We also derive a test condition to detennine the optimal signaling scheme to use. 

The achievable rate for one frequency bin can be written as 



rW 

Ra{sI{S)AU)A{S))= log: 

JQ 



1 + 



iVi + 4(/)Xi + W)Fi 



df. 



Thus, 



cl= m2^^RA{sl{f\si{f),slU)y 



(76) 



(77) 



We will consider the upstream capacity c" expression for our analysis. Further, we will use Ra 
for RA{sl{f), s^{f), sl{f)) in the remainder of this Section. Substituting for the PSDs from (72), 
(73) and (74) into (76) and using (77) we get the following expression for the upstream capacity 



W 

— max 

2 0.5<Q<1 



10g2 



1 + 



W 



w 



+ l0g2 



1 + 



(l-a)2P„igi 
W 



•i'l T -r IV 



w 



(78) 



Let Gi = and G2 = denote the SNRs in the bin due to line 1 and line 2 respectively. 
Then, we can rewrite (78) as 
W 



max — 

0.5<a<l 2 



M2 



1 + 



aGiHi 



1 + {1 - a)G2Xi + aGiFi 



+ log2 



1 + 



(1- 



1 + aGzX- 



1 + (1 - a)G2Fi\ j 



(79) 



Using (77) and differentiating the achievable rate (Ra) expression in (79) with respect to a 
gives us 

dRA 



da 



= {2a - 1) [GliXf -F?) + 2G2{Xi - Fx) - G^H^iG^Fr + 1)] L, 



(80) 



with L > 0 Va e (0, 1]. Setting the derivative to zero gives us the single stationary point a = 0.5. 
Thus, the achievable rate Ra is monotonic in the interval a 6 (0.5, 1] (see Figure 24). If the value 
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Of = 0.5 corresponds to a maximum of Ra, then it is optimal to perform EQPSD signaling in this 
bin. If the value a = 0.5 corresponds to a minimum of Ra, then the maximum is achieved by the 
value a = 1, meaning it is optimal to perform FDS signaling in this bin. No other values of a are 
an optimal option. 

The quantity a = 0.5 corresponds to a maximum o/Ra (EQPSD) if and only if ^ < 0 
Va e (0.5, 1]. For all a G (0.5, 1], ^ is negative if and only if (see (80)) 

GliXf - Ff) + 2G2{Xi - Fx) - GiHi{G2Fi + 1) < 0. 

This implies that 

^ GliX!-F^) + 2G2iX,-F,) ^gj^ 
G^FiHi + Hi 

In a similar fashion or = 0.5 corresponds to a minimum of Ra if and only if ^ > 0 
Vq; G (0.5, 1], This implies that a = 1 corresponds to a maximum (FDS) since there is only 
one stationary point in the interval a E [0.5, 1] (see Figure 24), For all o; e (0.5, 1], ^ is positive 
if and only if (see (80)) 

GliX^i - F^) + 2G2{Xi - Fi) » G.miG^Fi + 1) > 0. 

r 

This implies that 

^ gi(X?-i?) + 2G2(Xi-Fi) ^g2) 
G2F1H1 + Hi 

The above statements can be summed in a test condition to determine the signaling nature 
(FDS or EQPSD) in a given bin. Using (81) and (82) we can write 



2^ ^^>^° GI{X!-F!)+2G2{Xi-Fi) 
' FDS G2F1H1 + H, 

Thus, we can write the upstream capacity of line 1 in bin A; as 



(83) 



c?=< 



(84) 



K\os„\l4 2LniMi 1 ifa-l 

^ 2 ^"52 ^ NiW+2P^2Fi\ ' II a — 1. 
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4J2 Solution for M Unes: EQPSD and FDS signaling 



It is straightforward to generalize the result in the previous Section to M lines where each line i 
has parameters Hi, Gi, Pmu and Fi for z 6 {1, . - . , M}. Further, we assume that the self-NEXT 
and self-FEXT coupling transfer functions between lines 2, • - , M and line 1 are all the same. The 
test condition to determine signaling nature (EQPSD or FDS) in bin k of line 1 for M line case can 
be written as 



NiW 



FDS 



{i:fi2Gi)FiHx + Hi 



We can write the upstream capacity of line 1 in bin A as 



Wlog2 



1 + 



1 + 



iVlW+2(X;]^Pmi)FlJ ' 



, if a = 0.5, 



if a = 1. 



(85) 



(86) 



4.7.3 Solution for 2 Unes: EQPSD and multi-line FDS signaling 



We saw in Section 4.6.8 that in the case of two lines it is optimal to use multi-line FDS instead of 
FDS signaling. In this Section we will derive a test condition to determine the signaling nature in 
a given bin. We use the notation as introduced in Section 4.7.1. 

We desire a signaling scheme that supports multi-line FDS, EQPSD, and all combinations in 
between in a frequency bin. Therefore we divide each bin in half and define the upstream and 
downstream transmit spectra as follows (see Figure 42): 



a 



2P^ 



ifl/l<f, 



slif) = sfif) = { 



w 

iff <|/|<W, 



" 2 

Otherwise, 



(87) 



if!/|<f, 



a 



w 



iff <|/|<W, 
Otherwise, 



(88) 
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where 0.5 < a < 1. We assume that the upstream and downstream transmit spectra obey power 
complementarity, i.e., line 1 puts less power where line 2 puts more and vice versa. In further dis- 
cussion we will use transmit spectra s^(/) and (/)• When a = 0.5, sl[f) = 5^(/), V/ e [0, W] 
(EQPSD signaling); when a = 1, s?(/) and s^{f) are disjoint (FDS signaling). Tlie capacities of 
opposite directions are equal for each line: 

cl^4 andc^=:4 

The factor a controls the power distribution in the bin and W is the bandwidth of the bin. 

Next, we show that the optimal signaling strategy uses only EQPSD or multi-line FDS in each 
subchannel and derive a test condition to determine the signaling scheme to use. 

The achievable rate for one frequency bin can be written as 



RA{sl{f),si{f),st{f))=^ f log; 



1 + 



(89) 



then 



c?-^ma?c^iJ4(5?(/),5^(/),5^(/)). 



(90) 



We will consider the upstream capacity c" expression for our analysis. Further, we will use Ra for 
^a(5i (/), 4{f)y^^ *e remainder of this Section. Substituting for the PSDs from (72) and 
(73) into (89) and using (90) we get the following expression for the upstream capacity 



W 

— max 

2 0.5<a<l 



l0g2 



1 + 



a2PmiHi 
W 



AT t (l"a)2Pm2Xi , (l-a)2P^2Fi 
-i^i "t" "T w 



W 



+ l0g2 



1 + 



{l-a)2PmiHi 

W 



M 4. Q2fin2Xi , a2Pyrt2Fl 

1^1 -r -r ^ 



w 



(91) 



Let Gi = and G2 = ^^'^enote the SNRs in the bin due to line 1 and line 2 respectively. 
Then, we can rewrite (91) as 
W 

0.5<a<l 2 



c? = max 



|log2 [1 + 



aGiHi 



+ (1 - a)G2Xi + (1 - a)G2Fi 



+ log2 



1 + 



(1- 

1 + aGsX: 



X;i+Q!G2FiJj 



(92) 
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Using (90) and differentiating the achievable rate (Ra) expression in (92) with respect to a 
gives us 

^ = (2a - 1) [Gl(Xi + + 2G2(Xi + FO » dH,] L, (93) 

with L > 0 Vof e (0, !]• Setting the derivative to zero gives us the single stationary point a = 0.5- 
Thus, the achievable rate Ra is monotonic in the interval a E (0.5, 1] (see Figure 24). If the value 
a = 0.5 corresponds to a maximum of Ra, then it is optimal to perform EQPSD signaling in this 
bin. If the value a = 0.5 corresponds to a minimum of Ra, then the maximum is achieved by the 
value a = 1, meaning it is optimal to perform multi-line FDS signaling in this bin. No other values 
of a are an optimal option. 

The quantity a = 0.5 corresponds to a maximum of Ra (EQPSD) if and only if ^ < 0 
Va 6 (0.5, 1]. For all a 6 (0.5, 1], ^ is negative if and only if (see (93)) 

Gl{Xi + Fif + 2G2(Xi + Fi) - GiHi < 0. 

This implies that 

^ GUX, + F,f + 2G2{X, + F^) _ ^^^^ 

Hi 

In a similar fashion a = 0,5 corresponds to a minimum of Ra if and only if > 0 Va G 
(0.5, 1], This implies that a = 1 corresponds to a maximum of Ra (multi-line FDS) since there is 
only one stationary point in the interval o; G [0.5, 1] (see Figure 24). For all of G (0.5, 1], is 
positive if and only if (see (93)) 

GliXi + Fxf + 2^2(^1 + Fi) - GiFi > 0. 

This implies that 

Gl(Xx + F0^ + 2G2(Xx + Fi) ^^3^ 

The above statements can be sunmied in a test condition to determine the signaling nature 
(EQPSD or multi-line FDS) in a given bin. Using (94) and (95) we can write 

G, = ^ ^^>^^ (n{Xy + F{f^2G<,{X, + F,) ^^^^ 
multi - lineFDS 
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Thus, we can write the upstream capacity of line 1 in bin A; as 

Vlog, [14- ^,^^^-^^^^^ ifc. = 0.5, 
- I (97) 

4.8 Optimizing under a PSD mask constraint: No self-interference 

In this Section we will impose an additional peak power constraint in frequency, i,e., a limiting 
static PSD mask constraint. This implies that no transmit spectrum can lie above the PSD mask 
constraint. This constraint is in addition to the average power constraint We shall obtain optimal 
transmit spectra for an xDSL line under these constraints, in the absence of self-interference. 

4.5.1 Problem statement 

Maximize the capacity of an xDSL line in the presence of AGN and interference (DSIN-NEXT 
and DSIN-FEXT) from other services under two constraints: 

1 , The xDSL transmit spectra are limited by constraining static PSD masks; (/) for upstream 
and Q^{f) for downstream. 

2. The average xDSL input power in each direction of transmission must be limited to Ptnax 
(Watts). 

Do this by designing the distribution of energy over frequency (the transmit spectrum) of the xDSL 
transmission. 

4.8.2 Solution 

Consider a line (line 1) carrying an xDSL service. Line 1 experiences interference from other 
neighboring services (DSIN-NEXT and DSIN-FEXT) and channel noise No{f) (AGN) but no 
self-NEXT or self-FEXT (see Figure 19). 
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The twisted pair channel can be treated as a Gaussian channel with colored Gaussian noise 
[13]. Recall that DSnU) is the PSD of the combined DSIN-NEXT and DSpif) is the PSD of 
the combined DSIN-FEXT Let 5^(/) and denote the PSDs of line 1 upstream (u) direction 
and downstream (d) direction transmitted signals, respectively. Further, let and C*^ denote the 
upstream and downstream direction capacities of line 1 respectively. Let Hc{f) denote the channel 
transfer function of line L 



The channel capacities (in bps) are given by [14] 

\Hc{f)\' S-if) 



roo 

C = sup / l0g2 



and 



sup Plog 



oo 

2 



\Hcif)fSV) 



1 + 



df (98) 



df. (99) 



Noif)) + DSN{f) + DSF{f)\ 
The supremum is taken over all possible 5"(/) and iS^(/) satisfying the average power constraints 
for the two directions 

2/ 5«(/)rf/<P^ax and 2/ Sy)df<Pra^, (100) 
and the positivity and new peak power constraints 

r 

0<5"(/) <g"(/)V/ and 0 < 5^) < Q^(/) V/, (101) 

Note that these equations are the same as (4>-(6) except for the additional peak power constraint 
in frequency. For discussion purposes, we will focus on the upstream transmission. The same 
analysis can be appUed to the downstream channel. 

We wish to maximize (98) subject to the constraints (100), (101). The constraints (100), (101) 
are differentiable and concave, Rirther, the objective function to be maximized (98) is also con- 
cave (the log function is concave). Any solution to this problem must satisfy the necessary KKT 
(Karush-Kuhn-Tucker) [22] conditions for optimaiity. For a concave objective function and con- 
cave, differentiable constraints, any solution that satisfies the necessary KKT conditions is a unique 
globally optimal solution [22], Thus, we seek any solution that satisfies the KKT conditions, since 
it is automatically the unique optimal solution. 
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The optimal solution to (98), (99), (100), (101) is basically a "peak-constrained water-filling".* 
The optimal transmit spectrum is given by 

5o"pt{/) = | for/€£^„ax, (102) 

0 otherwise, 

with A a Lagrange multiplier. The spectral regions jEpos and are specified by 

^pos = {/ : 0 < < and 

^max - {/ : > Q^/)}. (103) 

We vary the value of A to achieve the optimal transmit spectrum S^^^if) that satisfies the average 
and peak power constraints (100), (101). It can be easily shown that this solution satisfies the KKT 
conditions for optimality. Substituting the optimal PSD SQ^^{f) into (98) yields the capacity C 
under the average and peak power constraints. 

Note that if the maximum allowed average power (Pmax) exceeds the power under the con- 
straining mask then the optimal transmit spectrum is the constraining PSD mask itself. In the 
absence of an average power constraint (but with a peak power constraint) the optimal transniit 
spectrum is again the constraining PSD mask. 



4.83 Examples 

In this Section we consider a line carrying HDSL2 service under the OPTIS [5] constraining PSD 
mask and input power specifications. An average input power (Pmax) of 19.78 dBm and a fixed bit 
rate of 1.552 Mbps was used fordl simulations. 

Figure 43 shows the optimal downstream transmit spectrum for HDSL2 with OPTIS down- 
stream constraining mask in the presence of DSIN-NEXT from 49 HDSL interferers and AGN 
(-140 dBm/Hz). The key features in the case of HDSL interferers are: 

^Peak-constrained water-filling can be likened to filling water in a closed vessel with uneven top and bottom 
surfaces. 
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L Comparing the peak-constrained transmit spectrum in Figure 43 with the unconstrained in 
peak power one in Figure 21 indicates that the peak-constrained optimal solution tries to 
follow the unconstrained in peak power optimal solution. The peak-constrained optimal 
solution has a null in the spectrum around 150 kHz similar to the one in the unconstrained in 
peak power spectrum. The null in the transmit spectra occurs in order to avoid the interfering 
HDSL transmit spectrum, 

2. An OPTIS transmit spectrum, achieved by tracking 1 dBm/Hz below the OPTIS PSD mask 
throughout, does not yield good performance margins (see Table 7). The OPTIS transmit 
spectrum looks different from the peak-constrained optimal spectrum (see Figure 43). The 
null in the peak-constrained optimal spectrum (which is not seen in the OPTIS transmit 
spectrum) indicates that it is suboptimal to distribute power according to the OPTIS transmit 
spectrum. 

Figure 44 shows the optimal upstream transmit spectrum for HDSL2 with OPTIS upstream 
constraining mask in the presence of DSIN-NEXT from 25 Tl interferers and AGN (- 140 dBm/Hz), 
Again, we compare the peak-constrained transmit spectrum in Figure 44 with the unconstrained in 
peak power one in Figure 22. Note that the peak-constrained optimal transmit spectrum puts no 
power in the high-frequency spectrum (to avoid Tl interference) as opposed to an OPTIS transmit 
spectrum. 

4.9 Optimizing under a PSD mask constraint: With self-interference 

The solution outlined in the previous Section applies only in the absence of self -interference. In 
this Section we will find an optimal transmit spectrum in the presence of additional self-NEXT 
and self-FEXT We will impose a peak power constraint in frequency, Ic, a limiting static PSD 
mask constraint, in addition to the average power and symmetric bit-rate constraints. We will 
obtain the optimal transmit spectra for an xDSL line under these constraints in the presence of 
self-interference. 
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Table 7: Uacoded performance margins (in dB) for CSA No. 6; OPTIS vs. Peak-constrained Optimal 
"under OPTIS" 



Crosstalk Sic 


xDSL service 


OPTIS 


Optimal 


Diff 


Dn 


Up 


Dn 


Up 


Dn 


Up 


49 HDSL 


HDSL2 


12.24 


2.7 


13.74 


3.74 


1.54 


1.03 


25 Tl 


HDSL2 


17.5 


19.9 


18.81 


20.43 


1.31 


0,53 


39 self 


HDSL2 


9.0 


2.1 


15.51 


17.58 


6.51 


15.48 


24self+24Tl 


HDSL2 


1.7 


4.3 


4.74 


4.52 


3.04 


0.22 



Bit rate fixed at 1.552 Mbps. 
Average Input power = 19.78 dBm, 

Diff (Dn) = Difference in Downstream margins (Optimal - OPTIS) 
Diff (Up) = Difference in Upstream margins (Optimal - OPTIS) 

4.9.1 Problem statement 

Maximize the capacity of an xDSL line in the presence of AGN, interference (DSIN-NEXT and 
DSIN-FEXT) from other services, and self-interference (self-NEXT and self-FEXT) under three 

r 

constraints: 

1 . The xDSL transmit spectra are limited by constraining static PSD masks; (/) for upstream 

Q^{f) for downstream, 

2. The average xDSL input power in each direction of transmission must be limited to Pmax 
(Watts). 

3. Equal capacity in both directions (upstream and downstream) for xDSL. 

Do this by designing the distribution of energy over frequency (the transnait spectra) of the xDSL 
transmissions. 

Additional assumptions are made in this case as given in Section 4.5.3 or 4.63 depending on 
the signaling scheme used. 



102 



4.9.2 Solution 



Consider a line (line 1) carrying xDSL service. Line 1 experiences interference from other neigh- 
boring services (DSIN-NEXT and DSIN-FEXT), channel noise No{f) (AGN), and self-interference 
(self-NEXT and self-EEXT) (see Figure 3). 

We need to find peak-constrained optimal transmit spectra for upstream and downstream trans- 
mission. We let the constraining PSD mask Q{f) be the maximum of the two upstream and down- 
stream constraining masks (<3"(/) and We then employ the solutions as described in 
Sections 4,5 or 4.6 but limit the peak power to the constraining mask <5(/). Thus, we obtain a 
pe^-constrained transmit spectrum 5opt(/). Using this mask, we optimally group the bins (see 
Section 4.5.10) to obtain optimal upstream and downstream transmit spectra (^f (/) and Si{f)). 



4.9.3 Algorithm for peak-constrained optimization of the transmit spectra 

1. Choose the constraining PSD mask as 

Q{f)=max{Q^{f),Qy)) V/. 

2. Solve for the optimal transmit spectrum S'^^^{f) according to the algorithms in Sections 
4.5.7, 4.5.8, or 4.6 with the following added constraint: 

^ Qif) V/ where 5-(/)>Q(/), 
5opt(/) = \ (104) 
[ S'^if) otherwise, 

where S^{f) is the water-filling solution (refer to [14, 16] if the spectral region employs 
EQPSD or multi-line FDS signaling and to [16] if the spectral region employs FDS signaling) 
(see Sections 4,5 and 4.6). This is the peak-constrained water-filUng solution in the presence 
of self-interference. As argued in the previous Section, this solution satisfies the necessary 
KKT conditions for optimality and therefore is the unique optimal solution. 

3. Denote the spectral region employing FDS signaling as Efds and the spectral region em- 
ploying EQPSD signaling as E'eqpsd- 
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Obtain Si^if) from S^^^{f) by symmetry, i.e., S^,{f) = S^^{f) in EQPSD and multi-line 
FDS regions and ± 5^pt(/) in FDS spectral regions. Merge Sip^if) and S^^{f) to 

form Soptif) as 

SoAf) = '5opt(/) = '5o'pt(/) V/in^EQPSD, 

5opt(/) = 5«pt(/)u5o^p,(/) V/inE-pDs, (105) 

where U represents the union of the two transmit spectra. 
Group the bins to obtain upstream and downstream masks as 

S^if) = Soptif) V/ in EpDs and where Q«(/)>Q''{/), 

Sfif) = Soptif) V/ in £;fds and where Q«(/)<Q''{/) (106) 

in £'fds aiid 

Sni) = sfif) = 5opt(/) V/ in ^EQPSD. (107) 

4. Check if the average power constraint is violated for upstream or downstream transmission, 

5. If the average power constraint is violated for direction o (i.e., the total transmit power in the 
direction o is more than Pmaxf then transfer power from 5i (/) to 5f (/). Transfer power first 
from spectral regions of to Sf{f) with the least Sf{f) ~ Sf (/) difference. Repeat this 
successively in spectral regions with increasing S^f) - S^f) difference until the average 
power in both directions is the same,^ We transfer power from one direction o to the other 
direction o in spectral regions where the difference in power between the two transmission 
directions is the least until the power between the two directions becomes equal This power 
transfer scheme is in a se^nse optimal as it tries to even out the powers between the two 
directions, with the least loss in the total sum of the transmit powers of the two directions. 

''Note that if the total traiivSmit power in direction o is more than Pmax then the transmit power in direction o is less 
thanPmax. 

^This approach of transferring power from direction o to direction o can be likened to "stealing from the rich and 
giving to the poor.'* 
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If the difference S^{f)—Sf{f) is the same (or marginally varying) for a range of frequencies, 
then transfer power from direction o to direction o in those spectral regions that give the 
maximum gain in bit rates for direction a 

43 A Examples and results 

In this Section we consider a line carrying HDSL2 service under the OFITS [5] constraining PSD 
mask and input power specifications. An average input power (Pmax) of 19,78 dBm and a fixed bit 
rate of 1.552 Mbps was used for all simulations. 

- Figure 45 shows the optimal upstream and downstream transmit spectra for HDSL2 with OP- 
US constraining masks in the presence of self-NEXT and self-FEXT from 39 HDSL2 interferers 
and AGN (—140 dBm/Hz). Note that the optimal upstream and downstream transmit spectra are 
separated in frequency (using FDS signaling) in a large spectral region in order to avoid high 
self-NEXT. On the other hand, OPTIS transmit spectra have a large spectral overlap at lower fre- 
quencies (self-NEXT is high here) that significantly reduces its performance margins (see Table 
7). 

Figure 46 shows the optimal upstream and downstream transmit spectra for HDSL2 with OP- 
TIS constraining masks in the presence of self-NEXT and self-FEXT from 24 HDSL2 interferers, 
DSIN-NEXT from 24 Tl interferers, and AGN (-140 dBm/Hz). Again, we see that the upstream 
and downstream optimal spectra are separated in frequency (using FDS signaling) over a large 
spectral region. However, the EQPSD spectral region towards the beginning of the spectrum is 
larger here than in the previous example, since we have more DSIN-NEXT from TL 

Key here is that optimal transmit spectra employ optimal separation in frequency of upstream 
and downstream services in the presence of interference. The "1 dB below OPTIS" transmit spectra 
do not do this, and so have inferior performance. 

Table 7 compares the performance margins of the OPTIS transmit spectra (obtained from the 
OPTIS PSD mask by uniformly subtracting 1 dBm/Hz over the entire frequency range as in [5]) 
with the optimal transmit spectra under the OPTIS PSD mask constraints. Table 7 shows that the 
optimal scheme significantly outperforms OPTIS in the case of self-interference. In cases involving 
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different service interferers (HDSL and Tl) the optimal scheme consistently outperforms OPTIS 
by 1 dB or more. Further, comparing these results with those in Table 1 suggests that the OPTIS 
PSD mask is not a good constraining PSD mask, since the unconstrained in peak power margins in 
Table 1 are significantly higher than the ones in Table 7. Comparing Tables 1 and 7 suggests that 
optimal signaling with no peak power constraint (static PSD mask) gives high peif ormance margin 
gains. 

4.10 Bridged taps 

Bridged taps (BTs) are short segments of twisted pairs that attach to another twisted pair that 
carries data between the subscriber and the CO. BTs are terminated at the other end with some 
characteristic impedance. BTs reflect the signals on the data-carrying line. These reflections de- 
structively interfere with the transmitted signal over certain frequencies. This leads to nulls in the 
channel transfer function and the self-FEXT transfer function at these frequencies (see Figure 48). 
These nulls in the channel transfer function significantly reduce the data transmission rate. Thus, 
bridged taps pose an important problem in achieving high bit rates over xDSL lines.^ 

Bridged taps presence, location, and length vary according to each loop setup. Thus, the effect 

r 

of BTs on the transmission signals is different for each loop. This means that the channel transfer 
function nulls (in frequency) vary for each separate line. We need to adapt the transmit spectrum to 
the channel conditions in order to achieve high bit-rates. We need the optimal power distribution 
that maximizes the bit-rates in the presence of bridged taps and interference. This further enforces 
the need for optimal dynamic transmit spectra and indicates that static transmit spectra are not 
a good idea. In this Section, we present optimal and near-optimal solutions to find the transmit 
spectra in the presence of BTs. • 

4.10.1 Optimal transmit spectra 

Optimal signaling is more computationally expensive to implement in the presence of bridged taps 
[3], as the channel transfer function has nulls and thus loses its monotonicity. In this scenario, 
^Bridged taps can be removed from xDSL lines, but this is an expensive (labor-intensive) procedure. 
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even the self-FEXT transfer function has nulls. In spite of this, the overall optimal solution can be 
obtained by a bin by bin analysis: 

1. Divide the frequency axis into narrow bins or subchannels. Compute channel transfer func- 
tion, various interference transfer functions, and AGN. 

2. Qhoose an initial power distribution of Pmax over all bins. 

3. Given the powers in each bin decide the optimal signaling scheme in each bin. Compute 
capacities for each bin and hence compute channel capacity. 

4. Re-distribute the powers in each bin by water-filling [14], [16], decide the optimal signaling 
scheme in each bin, and re-calculate the channel capacity. Repeat this step until we find 
the maximum possible channel capacity. It can be exceedingly computationally intensive to 

, find the optimal power distribution over all bins. There can be several local maxima for the 
channel capacity curve, and there is no guarantee that a search algorithm will converge to 
the global maximum. 

The optimal power distribution algorithm suggests that EQPSD, FDS, and multi-line FDS bins 
could be randomly distributed throughout the transmission bandwidth. The search for the optimal 
switchover bins from one signaling scheme to the other could be exceedingly expensive (involving 
a multi-dimensional search). 

4.10.2 Suboptimal transmit spectra 

We saw in the previous Section that the optimal transmit spectrum could be very expensive to 
obtain. However, we can always get a good suboptimal solution for line i as follows: 

1. Divide the jfrequency axis into narrow bins or subchannels as in Section 4.1. Compute chan- 
nel transfer function {Hc{f)), the various interference transfer functions (ifiv{/), Hp{f), 
DSnU), and DSF{f)), and AGN (iV^(/)). Obtain subchannel values (Hi^, X,-,^, Fi^k) for 
each bin using (l)-(3) and (13). Let k denote the bin number. 
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2. Use the condition evaluations in (26) and (27) to detennine the signaling scheme (EQPSD 
or FDS) in each bin* For each bin: 

• If (X?jfc - i^jfc ~ Hi^kFi^k < 0) and the right side of (26) < 0, then employ EQPSD 
signaling in that bin (since power in every bin > 0), 

• IfiXf f^-F^^^-Hi^kFi^k > 0) and the right side of (27) <0, then employ FDS signaling 
in that bin (since power in every bin > 0), 

• Employ FDS signaling if both the above conditions are not satisfied, 

3. Perform the optimal power distribution under average power constraint of Pmax using water- 
filling technique [14], [16]. 

4. Use condition evaluations in (46) and (54) to determine bins employing multi-line FDS. Re- 
distribute power optimally using water-filling technique. This step is optional and indicates 
which bins employ multi-line FDS signaling. 

The suboptimal solution determines the signaling strategy in each bin by simple, fast comparisons 
involving transfer functions and SNRs. This is followed by a simple optimal power distribution 
scheme using the water-filling technique. 

Note that the optimal and suboptimal algorithms can be implemented under a peak frequency- 
domain power constraint (static PSD mask). This is achieved by using peak-constrained water- 
filling technique (instead of just water-filling) for optimal power distribution (see Sections 4.8 and 
4.9) in the algorithms given in Sections 4.10.1 and 4.10.2. 

4.10.3 Examples and discussion 

Optimal transmit spectra: Theoretically, the optimal transmit spectrum in die presence of BTs 
can have several switchover bins fi:om one signaling scheme to the other (for e.g., EQPSD to 
FDS and FDS to EQPSD switchover bins). However, we argue that in most of the symmet- 
rical data-rate services (like HDSL2 and "VDSL2'') there is only one switchover bin from 
EQPSD to FDS inspite of bridged taps. 
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As frequency increases, the self-NEXT transfer function rapidly increases but the self-FEXT 
and the channel transfer functions generally decrease even for bridged taps case (see Figures 
17 and 48), Thus, the quantity Xf j^ - Ff^. - Hi^kFijc teiids to be an increasing function 
of frequency or bin number fc, and stays positive once it becomes positive. Similarly, the 
quantity Hi^k — ^{Xi^k — Fi^k) tends to decrease with frequency or bin number k and stays 
negative once it becomes negative. Using the condition evaluations (26) and (27) for all the 
frequency bins indicate that there is only one EQPSD to FDS switchover bin. Our studies 
indicate that is indeed true for a wide range of loops having bridged taps and employing 
HDSL2, "VDSL2" or similar synmietric services. The optimal switchover bin along with 
- the optimal transmit spectrum can be determined using the algorithm in Section 4,5.7. 

Figure 47 illustrates a case of "contiguous" optimal transmit spectra in case of a loop with 
bridged taps (CSA loop 4), We can clearly see that the optimal transmit spectra have only 
one transition region from EQPSD to FDS signaling. The transmit spectra were obtained 
such that we have equal performance margins and equal average powers in both directions 
of transmission. 

Suboptimal transmit spectra: We presented strong arguments in support of only one EQPSD to 
FDS switchover hm in the previous paragraph. However, there can be exceptions when the 
arguments do not hold, and we have multiple EQPSD and FDS regions (see Figure 48). 

Consider a hypothetical case of a short loop (1.4 kft with 3 bridged taps) carrying the 
"GDSL" service. The channel transfer function, self-NEXT, and self-FEXT transfer func- 
tions are illustrated at the top of Figure 48. Note that for "GDSL" service the self-NEXT 
is assumed very low. Since the self-NEXT is low, the non-monotonicity of the self-FEXT 
and the channel transfer function lead to distributed EQPSD and FDS regions across the 
transmission bandwidth as illustrated in the bottom of Figure 48. In such a scenario, the 
optimal power distribution algorithm of Section 4. 10, 1 is exceedingly difficult to implement. 
However we can easily implement the suboptimal solution as given in Section 4.10,2 
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4.11 Optimization: Asymmetrical data-rate channels 



Please refer to Section 5. 

4.12 Extensions 

4.12.1 More general signaling techniques 

The signaling techniques outlined earlier are not the only techniques that can give us improved 
capacity results. One possible scheme is illustrated in Figure 49. In this Figure, U Pi and DOWNi 
refer to line z, upstream and downstream direction PSDs respectively. In this scheme, we use multi- 
line FDS between group of lines (1 and 2) having high self-NEXT and high self-FEXT with other 
group of lines (3 and 4). However, there is EQPSD among group of lines (1 and 2 employ EQPSD 
as do 3 and 4) that have low self-NEXT and low self-FEXT within the group. This scheme can 
be extended for M self-interfering lines (with different self-NEXT and self-FEXT combinations 
between them) using combination of EQPSD, FDS, and multi-line FDS signaling schemes between 
different lines and frequency bins. 

The above scheme c^ be applied in the case of groups of lines with different self-interference 
(self-NEXT and self-FEXT) characteristics between different set of lines. 

4.12.2 More general interferer models 

If the self-NEXT and self-FEXT interferer model cannot be easily characterized by monotonicity 
in regions, (that is, if they vary rapidly and non-monotonously from one subchannel to the other), 
then we must search for the overdl optimal solution on a bin by bin basis. This search is outlined 
in the Section 4. 10 on bridged taps, 

4.12.3 Channel variations 

Some channels (e.g., the geophysical well-logging wireline channel) undergo a significant change 
in channel transfer function i?c(/) as a function of temperature. Temperature variations are a 
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part of nature and hence we need to continuously update our channel transfer functions. Changes 
in channel characteristics can change the channel capacity. We can develop an adaptive optimal 
transmit spectrum to adjust to these as well as any other variations. 

412.4 Broadband niodulation schemes 

We saw in Section 4.5. 10 that we can easily group the bins of the optimal transmit spectrum to make 
it smoother (with fewer discontinuities), so that we could apply different broadband modulation 
schemes. One can apply different broadband modulation schemes (like multi-level PAM, QAM, 
CAP, etc) over laige spectral regions to the optimal transmit spectrum obtained after grouping the 
bins and determine the performance margins. In this case, we need to use a DFE at the receiver 
to compensate for the severe channel attenuation characteristics. All these broadband modulation 
schemes do not suffer from latency as DMT does, but the DFE structure is complex. It is worth- 
while to compare the margins obtained with broadband' modulation schemes with those obtained 
using DMT as well as compare the complexity and implementation issues involved. 

412.5 Linear power constraints in frequency 

We saw in earlier Sections 4.4 — 4. 10, optimal power distribution using water-filling technique un- 
der an average power constraint, and peak-constrained water-filling technique under a peak power 
constraint in frequency or average plus peak power constraint in frequency. In general, we can 
determine the optimal power distribution under any set of general linear power constraints in fre- 
quency. Further, we can employ one of the joint signaling techniques discussed in this document 
under these new constraints using similar analysis, 

4.12.6 CDS signaling under a peak power constraint in frequency 

In case of a limiting static PSD mask, (see Sections 4.8 and 4.9), or odierwise, one may be required 
to limit the peak power in one or all the frequency bins. In this case a power-peaky signaling 
scheme like FDS or multi-line FDS will no longer be optimal as now we have a peak power 
constraint instead of the average power constraint. For this case, CDS or multi-line CDS signaling 
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[20] would be a better orthogonal signaling technique and would give increased capacity benefits 
without compromising spectral compatibility. 

Recall, that in frequency bins where self-NEXT is high and self-FEXT is low, we need orthog- 
onal signaling (FDS, TDS, or CDS) between upstream and downstream transmissions, i.e., 

5?(/)-L4(/), Vi^j. (108) 

Under an average power constraint, FDS signaling is the optimal signahng strategy (see Section 
4.5.12). In FDS signaling s^if) and 5^(/) occupy distinct separate frequency bands that are twice 
as higher than those using EQPSD signaling (see Figure 25). In CDS signaling the transmit spectra 
sf{'f) and s^{f) look similar to EQPSD signaling but the upstream and downstream spectra are 
separated using two orthogonal codes. Under a peak power constraint in frequency CDS signaling 
is preferred. Towards this end, we can group together bins using FDS into one spectral region 
-E'fds- We can implement spread spectrum CDS (SS-CDS) over this spectral region E'fds such 
that 

^(/)±5/(/), Vi^j, V/GJ?FDS. (109) 

Further, recall that in frequency bins where self-FEXT is high, we need to use orthogonal sig- 
naling (multi-line FDS, tDS, or multi-line CDS) between upstream and downstream transmissions 
of aU the M lines, i.e., 

s^{f)±s'j{f), yi^j, o&{u,d}. (110) 

Multi-line CDS separates the M interfering lines using M orthogonal codes and is less power- 
peaky in frequency than multi-line FDS. Under an average power constraint, multi-Une FDS sig- 
naling strategy is preferred. In multi-line FDS each line gets a separate frequency slot within each 
bin for transmission. The PSD_5?(/) in each bin is M times higher (or taller) than the corre- 
sponding PSD using EQPSD signaling (see Figure 18). Clearly, under a peak power constraint an 
alternative orthogonal signaling scheme like multi-line CDS is preferred. We can group together 
bins using multi-line FDS into one spectral region -E'mfds- We can implement SS-CDS over this 
spectral region E'mfds such that 

Stif) ± S]if), Vi # J, o e {u,d}, and V/ € j^mfds. (HI) 
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Note that implementatioa of SS-CDS cannot give perfectly orthogonal codes; instead we have 
only codes with very low cross-correlation. However, use of CDS or multi-line CDS signaling 
yields similar capacity (in the limit as cross-correlation between codes -> 0) as FDS or multi-line 
FDS schemes. 

4.12.7 Multi-user detector at central office 

We have seen that self-interference is a major limiter in achieving higher channel capacity. We 
can extend the work in previous Sections and construct a multi-user detector [21] at the central 
office that uses the self-interference for joint detection of each user (or line). In this sense the self- 
interference is not treated as only noise but can be used as information to achieve further significant 
gains in capacity of twisted pair lines. 

5 Optimized Signaling techniques for Asymmetrical data-rate 
channels 

Section 4 developed optimal signaling techniques for symmetric data-rate channels. In this Section, 
we will generalize the results in the previous Sections to general asymmetric data-rate channels. 
Asymmetrical data-rate channels have different upstream and downstream transmission rates, for 
e.g., ADSL and VDSL services. These channels also employ different average powers in the 
two transmission directions. We find joint signaling strategies and optimal power distribution for 
these channels using similar approaches as described in previous Sections (see Sections 4.5, and 
4.6). In this case we assume the knowledge of the ratio of average powers between upstream and 
downstream directions. ^ 

We will employ similar assumptions as in Section 4.1 with additional assumptions noted in 
particular cases. Further, we again employ joint signaling techniques comprised of one or more 
of the signaling schemes discussed in Section 43 with scaled-PSD (SPSD) signaling replacing 
EQPSD. SPSD is similar to EQPSD in that both upstream and downstream occupy the same band- 
width but we relax the constraint that both have to have the same PSD, instead one is a scaled 
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version of another in each narrow bin. 



5.1 Asymmetric data-rate optimization: Interference from other services 
(DSIN-NEXT and DSIN-FEXT) - Solution: SPSD Signaling 

In this scenario, each xDSL line experiences no self-interference (Figure 19 with neither self- 
NEXT nor self-FEXT). There is only DSIN-NEXT and DSIN-FEXT from other neighboring ser- 
vices such as Tl, ADSL, HDSL, etc., in addition to AGN. The solution is well known (classical 
water-filling [14]), but the idea is to optimally distribute power in a dynamic fashion. 

S.IA Problem statement 

Maximize the upstream and the downstream capacity of an xDSL line in the presence of AGN and 
interference (DSIN-NEXT and DSIN-FEXT) from other services under the constraint: 

1. The average xDSL input power in the upstream and downstream dhrections of transmission 
must be limited to Pxjp (Watts) and (Watts). 

r 

Do this by optimally designing the distribution of energy over frequency (the transmit spectrum) 
of the xDSL transmission. 

5.1.2 Additional assumption 

We add the following assumption to the ones in Section 4,1 for this case: 

12. Both directions (upstream^d downstream) of transmission experience the same channel 
noise (AGN) and different service interference (DSIN-NEXT and DSIN-FEXT). 

5.1.3 Solution 

Consider a line (line 1) carrying xDSL service. Line 1 experiences interference in both directions 
from other neighboring services (DSIN-NEXT and DSIN-FEXT) and channel noise No{f) (AGN) 
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but no self-NEXT or self-FEXT (see Figure 19), 

The DSIN-NEXT and DSIN-FEXT interference can be modeled as colored Gaussian noise for 
calculating capacity [13], Recall that DSnH) is the PSD of the combined DSESf-NEXT and let 
DSfU) is the PSD of the combined DSIN-FEXT. Let S^{f) and S^{f) denote the PSDs of line 
1 upstream (u) direction and downstream (d) direction transmitted signals respectively. Further, 
let and denote the upstream and downstream direction capacities of line 1 respectively. Let 
Hcif) denote the channel transfer function of line 1, The twisted pair channel is treated as a 
Gaussian channel with colored Gaussian noise. In this case the channel capacity (in bps) is given 
by [14] 



roo 

C" = sup / l0g2 



\Hc{f)\' S-if) 



and 



1 + 



No{f) + DSN{f)+DSFif) 
\Hc{f)fS'{f) 



df (112) 



No{f)) + DSM)+DSF{f)\ 



df. (113) 



C^= sup /°°log. 
The supremum is taken over all possible S^{f) and S'^{f) satisfying 

5"(/)>0V/, 5V)>0V/, 
and the average power constraints for the two directions 

2 r 5«(/)d/ < Pup, and 2 H S''{f)df < Fdn- (114) 

^0 JO 

Note that in general both directions experience different AGN, DSIN-NEXT and DSIN-FEXT 
interference. The upstream and downstream capacities can be independently calculated by sepa- 
rately maximizing (112) and (113). For illustration purposes we optimize the upstream transmdt 
PSD 5~(/). 

The optimal power distribution in this case is obtained by the classical "water-filling" tech- 
nique [14]. The optimal S'^if) is given by 

5«p.(/)= IW '^'^^^ (115) 

[ 0 otherwise, 

with A a Lagrange multipHer and E the spectral region where S^{f) > 0. We vary the value of 
A such tiiat 5jpt(/) satisfies with equality the average power constraint in (114). The equality is 
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satisfied for a single value of A giving us a unique optimal PSD S^^^if), Plugging the optimal PSD 
•^opt(/) (112) yields the capacity under the average power constraint. This procedure yields 
a unique optimal upstream transmit spectrum 5qp^(/) [14], Similarly we can find the optimal 
downstream transmit spectrum S^^^{f). The transmit spectra in general. are different (5^(/) ^ 
S^{f)) for the opposite directions of transnussion. 

Figure 20 gives a flowchart to obtain the optimal transmit spectrum in each direction using 
only SPSD signaling in the presence of DSIN-NEXT, DSBSf-FEXT and AGN, It uses the classic 
water-filling solution to obtain the transmit spectrum. The novelty is in applying this to xDSL 
scenario to achieve a dynamic transmit spectrum (different for each interference type) in the two 
dkections of transmission, 

5,1.4 Examples 

In this Section, we present some examples for the ADSL service [9]. We used an upstream average 
input power (over positive frequency) (0.5 * Pup) of 13 dBm and a downstream input power (over 
positive frequency) (0.5*Pdn) of 20.4 dBm, An uncoded performance margin of 6 dB was fixed for 
both upstream and downstream. We measured the achievable rates for the two directions. Figure 
50 shows the optimal up'stream and downstream transmit spectrum for ADSL in the presence of 
DSIN-NEXT from 49 HDSL interferers and AGN (-140 dBm/Hz). Note the deep null in the 
upstream transmit spectrum from approximately 110 to 220 kHz. This results fi:om "water-filling'" 
— the peak of the first main lobe of HDSL lies in the vicinity of 80 to 255 kHz. 

Figure 51 shows the optimal upstream and downstream transmit spectra for ADSL in the pres- 
ence of DSIN-NEXT from 25 Tl interferers and AGN (-140 dBm/Hz), 

The optimal transmit spectra for the two OHDSL and Tl interferers) cases are signifi- 
cantly different, evidence of the fact that the optimal transmit spectra will change depending 
on the nature of the interference* 

Table 8 lists achievable bit rates (in Mbps) for ADSL service in presence of DSIN-NEXT from 
HDSL and Tl interferers. Note that a desired ratio (r) between the upstream and downstream rates 
can be achieved by appropriately changing the input powers Pyp and Pdn- 
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Table 8: Achievable bit rates (in Mbps) using SPSD signaling 



Crosstalk source 


xDSL service 


Bit rates 


Ratio of rates 








49 HDSL 


ADSL 


4.39 


6.61 


11.00 


1.51 


25 Tl 


ADSL 


2.25 


2.95 


5.20 


1.31 



Uncoded performance margin fixed at 6 dB. 

Summary: We have seen that the optimal transmit spectra vary significantly with the interfer- 
ence combination. The water-filling solution yields unique transmit spectra for each interference 
combination [14], The transmit spectra adapt to minimize the effects of the interference combina- 
tion. The optimal transmit spectra for upstream and downstream direction are scaled replicas of 
each other in narrow bins, i.e. they employ SPSD signaling. 

5.2 Asymmetric data-rate optimization: Interference from other services 
(DSIN-NEXT and DSIN-FEXT) plus self-interference (self-NEXT and 
low self-FEXT) - Solution: SPSD and FDS Signaling 

In this scenario each asymmetric data-rate xDSL line experiences self-interference (high self- 
NEXT and low self-FEXT) in addition to AGN and DSIN-NEXT and DSIN-FEXT from other 
services (see Figure 3) in a generic xDSL service. This is the case of interest for the ADSL and 
VDSL services. 

We assume self-NEXT dominates self-FEXT and self-FEXT is not very high. So, we will 
design a system that rejects self-NEXT using orthogonal signaling as described in Section 4.5.1, 

5.2*1 Problem statement 

Maximize the sum of the upstream and the downstream capacity of an xDSL line in the presence 
of AGN, interference (DSIN-NEXT and DSIN-FEXT) from other services, and self-NEXT and 
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self-FEXT under the constraints: 

L The average xDSL input power in the upstream and downstream directions of transmission 
must be limited to Pup (Watts) and Pdn (Watts). 

Do this by designing the distribution of eneigy over frequency (the transmit spectrum) of the up- 
stream and downstream xDSL transmissions. 

We can also add an additional constraint of having the upstream and downstream powers in 
the SPSD region to have the same ratio as the input powers (see Section 5.2.6), 

5.2.2 Additional assumptions 

We add the assumptions mentioned in Section 4.5.3 to the ones in Section 4,1 for this case. 

5.2.3 Signaling scheme 

Since the level of self-NEXT will vary with frequency (recall Figure 17), it is clear that in high 
self-NEXT regions of the spectrum, orthogonal signaling (FDS, for example) might be of use in 
order to reject self-NEXT. However, in low self-NEXT regions, the loss of transmission bandwidth 
of FDS may outweigh any gain in capacity due to self-NEXT rejection. Therefore, we would like 
our signaling scheme to be general enough to encompass both FDS signaling, SPSD signaling, and 
the spectrum of choices in between. Our approach is similar to that in Section 4.5.4. 

Key to our scheme is that the upstream and downstream transmissions use different transmit 
spectra. All upstream (to CO) transmitters 27* transmit with the spectrum S^{f) and all down- 
stream (from CO) transmitters I^transmit with the spectrum S^{f), Implicit in our scheme is the 
fact that in this case, self-NEXT dominates selfFEXTand self-FEXT is small. If not, it would not 
be wise to constrain all T^' to the same transmit PSD. 

Our goal is to maximize the sum of the upstream capacity (C^) and the downstream 
capacity (C*^) given an average upstream and downstream total power constraint of Pup and 
Pdn respectively. 
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Consider the case of two lines with the same service. Line 1 upstream capacity is and line 2 
downstream capacity is C^. Under the Gaussian channel assumption, we can write these capacities 
(in bps) as 



roo 

sup / l0g2 

and 

roo 

sup / logs 

S«(/),S<'(/)-/0 



1^ \Hcif)f S^if) 



df, 
(116) 



df. 
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The supremum is taken over all possible 5"(/) and S^{f) satisfying 

5"(/)>0V/, 5'^(/)>0V/, 
and the average power constraints for the two directions 

' roo roo 

2 / 5'«(/)d/ < Pup, and 2 / S^{f)df < Pdn- (118) 
Jo Jo 

We can solve for the capacities and using "water-filling" if we impose the restriction 
of SPSD. However, this gives low capacities. Therefore, we employ FDS (5^(/) orthogonal to 
5'^(/)) in spectral regions where self-NEXT is large enough to hmit our capacity and SPSD in the 
remaining spectrum. This gives much improved performance. 

We assume the subchannel frequency responses (l)-(3) and the notation introduced in (12) 
and (13). Let s^{f) denote the PSD in bin k of line 1 upstream direction and s'^{f) denote the PSD 
in bin k of line 2 downstream direction (recall the notation introduced in Section 4.1, Item 9). The 
corresponding capacities of the subchannel k are denoted by c" and c^. 

We desire a signaling scheme that includes FDS, SPSD and all combinations in between in 
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each frequency bin. We divide each bin of width W in two parts^°: ^ and and define the 
upstream and downstream transmit spectra as follows (see Figure 52): 



if|/l<& 



0 otherwise 



(119) 



and 



g(l-a)(ft+X)Pm -f ifi ^ W 



h*W 



if^^ 



(120) 



.-M < I/I < 

0 otherwise. 
Here and 5 * F^^ are the average upstream and downstream powers over frequency range 
[OjVT] in bin A: and 0.5 < o: < 1 (Note that when a = 0.5, s"(/) and s^(/) occupy the entire 
bin V/ 6 [0, (SPSD signaling); when a = 1, 5^(/) and 5^(/) are disjoint (FDS signaUng). 
These two extreme transmit spectra along with other possible spectra (for different values of a) are 
illustrated in Figure 52. The PSDs 5"(/) and s^{f) are "symmetrical" or power complementary 
to each other. This tries to ensure the the upstream and downstream capacities scale according to 
some desired ratio. 

Next, we show that given this setup, the optimal signaling strategy uses only FDS or SPSD in 
each subchannel under some additional constraints. 



5.2.4 Solution: One frequency bin 



Recall the definition of the achievable mte firom (16) and that of c" and c"^ from (17). Substituting 
for the PSDs from (119) and (120) into (16) and using (17) we get the following expressions for 
the upstream and downstream capacities 

W 



max 




1 + 



a(h~i-l)Prr.H 

w 



AT t (l-a)gih^-l)Pn.X , a(h+l)PmF 
^ W W 

(l-ft)(fe+llPmff 

hW 



+ 



^ hW "T" hW 



!) 



(121) 



^^'The bin split-up does not necessarily have to be g% to the left side of the bin and (1 - g)% to the right side of 
the bin as shown in Figure 52, In general any 5% - (1 - ^)% power-complementary (using a and (1 - a)) split-up 
between opposite direction bins will work. 
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W 

max < - — - 
0.5<a<l \h-\-l 

hW 



1 + 



[l-a)9(h+l)P,nn 

w 



N + 



h + 1 



10g2 



1 + 



w 

ag(/i+l)fn.g 
hW 



+ 



(l-Q)g(h+l)PrnF 
W 



+ 



hW 



(122) 



Let G = ^^^jf"* denote the SNR in the bin. Then, we can write the sum of the capacities 
(ctot = c« + c^) using (121) and (122) as 

W f, L aGH 



Ctot 



h + 

/ilog2 



- max < logo 1 + 7 

1 0.5<a<l ( L 1 



loga 
/iloga 



1 + 
1 + 



+ 5(1 - a)GX + aGF 
(1 - a)GH 



+ 



h + gaGX + {l-a)GF 
g{l - a)GH 



+ 



1 + 



1 + aGX + 5(1 - a)GF 
gaGH 



+ 



/i + (1 - ot)GX + gaGF 
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(123) 



Note: To get a truly global maximum for Ctot we need to jointly maximize (123) over a, K and 
g. This is a hard problem to solve analytically. But we could numerically solve it for the range of 
bins by using constrained maximization routines over several variables if computational resources 

r 

are available. We note this but take the alternative approach of finding some sort of analytical 
relationship between the variables to decide the signaling scheme in each bin. 

We found that the problem can be solved if we assume that the power scaling factor between 
upstream and downstream (5), is equal to the scaling factor between the bin widths {h), i.e., g = h. 
This is a reasonable approach since in FDS regions this would give us the same scaling ofg = h 
between the upstream and downstream rates. Substituting for g in (123) we get 

aGH 



/il0g2 



logs 

/ll0g2 



1 + 
1 + 



+ h{l-a)GX + aGF_ 
(1 - a)GH 



+ 



h + haGX + (1 - a)GF 
h{l - a)GH 



+ 



1 + 



l + aGX + h{l-a)GF 
haGH 



+ 



h + {l-a)GX + haGF 
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We define the quantity after max in the above equation as i^Atot, i-e., the total achievable rate 
between the upstream and downstream directions. We search for the maxima of i^Atot in the above 
problem by treating it first as an unconstrainted optimization problem. We follow similar steps 
as described in Section 4.5.5, i.e., we differentiate i?Atot with respect to a and equate it to zero. 
Still, we are unable to solve the problem analytically. But, numerical studies show that we have a 
maximum of i^Atot at either a = 0.5 or a = 1. Thus, we study the behavior of the derivative at 
a = 1, i.e., ^^q"^ [^_^» We determine whether this derivative is positive {a = 0.5 is the maximum 
point) or negative (a = 1 is the maximum point). 

Interestingly the signaling scheme-deciding conditions that we get are similar to the ones in 
(26) and (27) and can be stated as 



If X^-F^-HF <Q then 
SPSD 



^ {h + l)Pm H-2{X-F) 
FDS X^-F-^-HF- 



NW 



(125) 



G 



Jf X^-F^-HF>0 then 

{h + l)Pm ^ H-2iX-F) 
NW pgg X^-F^-HF' 



(126) 



Thus, we can write the upstream and downstream capacities c" and c** in a frequency bin k as 



' W'l0g2[l + 



c" = < 



NW+gPmX+PraF) 



ifQ; = 0.5, 



(127) 



w 

h+1 



lo& [l + iT^Pr^] , if«=l. 
' ^log2 [l + ^.h^JXp,^) ] , if a = 0.5, 



(128) 



ifa = 1. 



5^.5 Solution: All frequency bins 

We can extend this joint signaling technique to all frequency bins similar to that described in 
Section 4,5.6. The algorithm for doing so is similar to the one described in Section 4.5.7 with the 
following exceptions: 
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1. We maximize the joint capacity of the lines, i^e,, + C^. 

2. Note that the ratio between upstream and downstream powers can vary with each bin. In the 
SPSD signaling region we can additionally impose the restriction that each bin should have 
the same ratio between upstream and downstream powers equal to 5spsd where 5spsd = 
fftot ~ where ^tot is the ratio between the upstream and downstream input powers. 

A fast, suboptimal solution can be obtained by using Me as the SPSD to FDS swith-over bin 
as described in Section 4.5.8. Further, we can obtain contiguous spectra by using the grouping 
ideas discussed in Section 45.10. 

5.2.6 Examples and results 

In this Section, we present some examples and results for the ADSL service. AGN of -140 
dBm/Hz was added to the interference combination in all simulations. We present results for joint 
capacity maximization under the following two scenarios: 

Fixed power ratio in SPSD bins: We constrain the SPSD bins to have the same downstream to 
upstream power ratio (^spsd) as the ratio between the total downstream and upstream powers 
{gspsD == Qtot = Table 9 lists achievable rates (in Mbps) obtained from our simula- 
tions for CSA loop number 6, We used an upstream average input power (over positive 
frequency) (0.5 * P^p) of 13 dBm and a downstream input power (over positive frequency) 
(0.5 * Pjy^) of 20.4 dBm. A fixed 6 dB of uncoded performance margin was used for all 
simulations along with a BER of 10"^. The last column lists the ratio of rates (r) between 
the downstream and upstream rates. We see significantly higher joint capacities than 6 Mbps 
in presence of most interference combinations. 

Table 10 shows the difference between the optimal solution (using the optimal Me2f) and the 
fast approximate suboptimal solution (using Me2f = Me) for a variety of interfering lines. 
The column Diff (in Mbps) shows the difference in achievable rates between the optimal and 
the suboptimal scheme. We make a similar note as in Section 4.5.11 that the suboptimal 
scheme gives near-optimal results that requires very little computational effort. 

123 



Table 9: Achievable bit rates (in Mbps) using SPSD and FDS signaling for CSA No. 6 



Crosstalk source 


xDSL service 


Bit rates 


Ratio of rates 
' c 








Iself 


ADSL 


1.99 


8.40 


10.39 


4.22 


10 self 


ADSL 


1.65 


7.34 


8.99 


4.46 


19 self 


ADSL 


1.54 


6.98 


8.52 


4.53 


29 self 


ADSL 


1.48 


6.72 


8.20 


4.54 


39 self 


ADSL 


1.43 


6.55 


7.98 


4.58 


10 self+ 10 HDSL 


ADSL 


1.41 


5.65 


7.06 


4.00 


10 self+ 10 Tl 


ADSL 


0.74 


2.37 


3.11 


3.19 



Uncoded performance margin fixed at 6 dB. 



Table 10: Achievable bit rates (in Mbps) for CSA No. 6; Optimal vs. Suboptimal, 



Crosstalk source 


xDSL 
service 


Optimal 
scheme (Mbps) 


Me2F 


Fast, suboptimal 
scheme (Mbps) 


Me 


Diff 


1 self 


ADSL 


10.39 


16 


10.39 


16 


0 


10 self 


ADSL 


8.99 


12 


8.99 


12 


0 


19 self 


ADSL 


8.52 


11 


8.52 


11 


0 


29 self 


ADSL 


8.20 


11 


8.20 


11 


0 


39 self 


ADSir 


7.98 


10 


7.98 


10 


0 


lOself+lOHDSL 


ADSL 


7.06 


30 


6.98 


12 


0.08 


10 self+lOTl 


ADSL 


3.11 


12 


3.11 


12 


0 



Uncoded performance margin fixed at 6 dB. 
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Table 11: Achievable bit rates (in Mbps) using SPSD and FDS signaliag for CSA No, 6: Variable power 
ratio in SPSD bins 



Crosstalk source 


xDSL service 


Bit rates 


Ratio of rates 

r = ^ 

' c« 


Diff 


C" 






Iself 


ADSL 


5.01 


5.72 


10.74 


1.42 


0.35 


10 self 


ADSL 


4.36 


5.67 


10.02 


0.77 


1.03 


19 self 


ADSL 


5.83 


3.95 


9.78 


0.68 


1.26 


29 self 


ADSL 


5.84 


3.75 


9.60 


0.64 


1.40 


39 self 


ADSL 


5.88 


3.57 


9.45 


0.61 


1.47 


10 self+ 10 Tl 


ADSL 


1.45 


1.76 


3.20 


1.22 


0.09 



Uncoded performance margin fixed at 6 dB, 

Variable power ratio in SPSD bins: We relax the constraint of fixed ratio of powers in SPSD 
bins in our optimization process. In this case, see better performance but the elegance of 
uniform complementarity of powers in each bin amongst SPSD or FDS bins is lost. 

Table 1 1 lists achievable rates (in Mbps) obtained from our simulations for CSA loop number 
6. We used similar simulation conditions as in the case for fixed power ratio case. The last 
colunm Usts the difference (Diff) between the variable power ratio and fixed power ratio 
cases. We note that variable power ratio case always performs better than the fixed ratio 
case. In the variable power ratio we mostly observe that the upstream and downstream 
power spectra employ FDS and occupy two distinct contiguous bandwidths. 

All the transmit spectra are spectrally compatible with existing services by design. This is true 
again, due to water-filling techniques that optimize power-distribution without causing significant 
interference in existing neighboring services. 
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52J Extensions 



Based on the work described in earlier Sections we can easily write the following extensions for 
the asymmetrical data-rate signaling scheme: 

1. In case of significant self-FEXT we can apply signaling techniques incorporating multi-line 
FDS as described in Section 4,6, 

2. The ideas on joint signaling in presence of differing channel, noise and interference char- 
acteristics as described in Section 4.7 can also be easily extended to asymmetrical data-rate 
lines. 

3. The optimization can also be easily done under a PSD mask constraint as described in Sec- 
tions 4.8 and 4.9. 

4. The results on bridged taps (see Section 4,10) and the extensions in Section 4.12 hold for 
asymmetrical data-rate lines. 

6 Summary of Contributions 

The key differences from the prior art are: 

L hicreased capacity for xDSL lines using optimal and suboptimal transmit spectra involving 
joint signaling schemes. 

2. "Symmetrical " ( or power complementary) upstream/downstream optimal transmit spectrum 
for a xDSL line in presencejof self-NEXT, self-FEXT, AGN, and otiier interfering lines like 
Tl, HDSL, and ADSL using EQPSD and FDS signaling. 

3. Fast near-optimal solution for the transmit spectrum which is computationally very attractive 
and very easy to implement for xDSL lines. 

4. Spectral optimization gives good spectral compatibility with other services (FDS better than 
CDS for spectral compatibility under an average power constraint). 
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5. Dynamic transmit spectram that adjusts automatically according to the interference type, 

6. Multi-line FDS signaling technique to combat self-FEXT. 

7. Increased capacity for HDSL2, "GDSL", and "VDSL2" lines using multi-Une FDS signaling 
when appropriate, 

8. Increased capacity in generic xDSL lines when neighboring lines have different channel, 
noise and interference characteristics, 

9. Concept of static estimation of interference values by reading look-up table of the topology 
of the cables (which self-interfering lines are where) at powerup. The self-interference values 
can be estimated in this manner. Dynamic measurement of interference values is done by 
"listening" to the interference during powerup, (Subtract the estimated self-interference from 
this measured interference to get the different service interference.) 

10, We can also interpret our results as capacity estimates given a fixed margin in the presence 
of fixed interferers. 

11, Extension of above points to asymmetrical data-rate channels like ADSL and VDSL. 

r 

Final notes; 

1. We have framed our work within the context of the HDSL2, "GDSL", and "VDSL2" trans- 
mission formats. However, our results are more general, and apply to all channels that exhibit 
crosst^ interference from neighboring channels. We summarize a few channels where this 
technique could be potentially applied: 

(a) Twisted pair lines (standard telephone lines) 

(b) Untwisted pairs of copper lines 

(c) Unpaired cables 

(d) Coaxial cables 

(e) Power Unes 
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(f) Geophysical well-logging telemetry cables 

(g) ^^eless channels. 

2. If a static mask is desired (e.g., for ease of implementation), we propose that a thorough study 
be made of the optimal solutions in different interference and noise scenarios as proposed in 
this document and then a best static compromising PSD mask be chosen. 
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Glossary 

ADSL: Asymmetrical digital subscriber line 
AGN: Additive Gaussian noise 
BER: Bit error rate (or probability) 
BT: Bridged tap 

CAP: Carrierless amplitude/pulse modulation 

CDMA: Code-division multiple access 

CDS: Code-division signaling 

CO: Central office 

CSA: Carrier serving area 

DFE: Decision feedback equalization 

DMT: Discrete multitone technology 

DSL: Digital subscriber line 

EQPSD: Equal pOAver spectral density signaling 

FDS: Frequency division signaling 

FEXT: Far-end crosstalk 

"GDSL": General digital subscriber line 

HDSL: High bit-rate digital subscriber line 

BDOSLl: High bit-rate digital subscriber line 2 

ISDN: Integrated services digital network 

ESI: Intersymbol interference 
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MFDS: Multi-line Frequency division signaling 

NEXT: Near-end crosstalk 

PAM: Pulse amplitude modulation 

POTS; Plain old telephone services 

PSD: Power spectral density 

QAM: Quadrature amplitude modulation 

SNR: Signal to noise ratio 

SPSD: Scaled power spectral density signaling 

SS-CDS: Spread spectrum code-division signaling 

Tl: Transmission 1 standard 

TDS: Time division signaling 

VDSL: Very high bit-rate DSL 

"VDSL2": Very high bit-rate DSL 2 

xDSL: Any generic DSL service 
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Notation 

±: Orthogonal 
U: Union 

A: Kind of service, such as ADSL, HDSL, HDSL2, VDSL, etc, 
jB: Channel transmission bandwidth 
Ci Channel capacity or line capacity 
D: -Difference between two capacities 
E: Spectral region 

F: Magnitude squared Far-end crosstalk (self-FEXT) transfer function in a single bin 

G: Signal to noise ratio (SNR) in a single bin 

H: Magnitude squared channel transfer function in a single bin 

J : Kind of signaling scheme, such as EQPSD, FDS, multi-line FDS, etc, 

r 

K: Total number of bins within channel transmission bandwidth 

L: Function of line parameters (G, F, X, if) in a single bin; it is always a positive quantity 
Ml Number of interfering lines carrying the same service 

N: Total additive Gaussian noise (AGN) power plus total different service interference 
P: Power 

Q: Constraining PSD mask 
R: Receiver 

S: Power spectral density (PSD) 
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T: Transmitter 

U : Positive quantity equal toY + Z + N 
V: Positive quantity equal to r + Z + iV + 5 
Wi Bandwidth of a bin or a subchannel 

X: Magnitude squared Near-end crosstalk (self-NEXT) transfer function in a single bin 

Y: Part of crosstalk power that couples into another service line 

Z: Part of crosstalk power that couples into another service line 

a: An amplitude level of a transmit spectrum 

b: An amplitude level of a transmit spectrum 

c; Capacity of a bin or a subchannel 

rf: Downstream direction 

/: Frequency 

g: Ratio between downstream and upstream powers in a single bin (g = ~) 

h: Ratio between downstream and upstream bandwidths in a single bin employing FDS 

i: Line number 

j: Line number 

ki Bin index 

n: Fraction to choose power distribution, 0 < n < 1 
o: Direction index o e {tx, d} 

r: Ratio between downstream and upstream rates (r = ^) 
ui Upstream direction 
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: Capacity of line i in transmission direction o 
C^: Capacity of a line in direction o 
Cii Capacity of a line i 

E-pDs: Spectral region employing FDS signaling 

EuFDs- spectral region employing multi-line FDS signaling 

Fi^. Magnitude squared self-FEXT transfer function on line i and bin k 

Fit Magnitude squared self-FEXT transfer function of line i in a single bin 

Gji Ratio of signal power in line i to noise power in line 1 in a single bin 

Hiy. Magnitude squared channel transfer function of line i and bin k 

Hii Magnitude squared channel transfer function of line i in a single bin 

Noif): Channel noise 

Nil AGN plus different service interference on line i 

Pjnii Power in positive frequency range ([0, W]) of a single bin of line i 

Pm: Power in positive frequency range ([0, W]) of a single bin 

^^max^ Total average power over the entire frequency range ([-B, B]) of the channel 

Q^ifY- Constraining PSD mask in direction o 

Ra^ Achievable rate in a single bin or subchannel 

R^: Receiver on line i in direction o 

Sf{f): PSD of line z in direction o 

S'^if): PSD of a line in direction o 

TP: Transmitter on line i in direction o 
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Xiy. Magnitude squared self-NEXT transfer function on line i and bin k 
Xii Magnitude squared self-NEXT transfer function on line i in a single bin 
^i, J* Capacity of a single bin of line i using signaling scheme J. 
cf: Capacity of a single bin of line i in direction o 
c"": Capacity of a single bin in direction o 
Si (/) : PSD in a single bin of line i in direction o 
PSD in a single bin in direction o 
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WHAT IS CLAIMED IS: 

1. A communications system that comprises: 

a subscriber modem; 
a central office modem; and 
5 a communications channel coupled between the subscriber modem and the central 

office modem and configured to transport uplink signals fi'om the subscriber modem to 
the central office modem, and fiirther configured to transport downlink signals fi'om the 
central office modem to the subscriber modem, 

wherein the power spectral density of the transmitted uplink signals is 
10 proportional to the power spectral density of the transmitted downlink signals. 

2. A communications system that comprises: 

a subscriber modem; 
15 a central office modem; and 

a communications channel coupled between the subscriber modem and the central 

office modem and configured to transport uplink signals fi'om the subscriber modem to 

the central office modem, and fiuther configured to transport downlink signals fi'om the 

central office modem to the subscriber modem, 
20 wherein at frequencies below a selected frequency Me2f? the power spectral 

density of the transmitted uplink signals is proportional to the power spectral density of 

the transmitted downlink signals by a positive scale factor, and 

wherein at frequencies above Meif, the power spectral density of the uplink 

signals are limited to one or more uplink frequency bands and the downlink signals are 
25 limited to one or more downlink frequency bands that are disjoint from the uplink 

frequency bands, and wherein the total bandwidth of the uplink frequency bands is 

proportional to the total bandwidth of the downlink frequency bands by the same positive 

scale factor. 

30 
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3. A communications system that comprises: 
a subscriber modem; 
a central office modem; and 

a communications channel coupled between the subscriber modem and the central 
office modem and configured to transport uplink signals fi'om the subscriber modem to 
the central office modem, and fiirther configured to transport downlink signals fi:om the 
central office modem to the subscriber modem, 

wherein when the connection is initialized, firequency bands are allocated to the 
uplink and downlink power signals so that the total uplink and downlink capacity is 
maximized over the channel for predetermined uplink and downUnk average signal 
powers, 

wherein the predetermined uplink and downlink average powers are unequal. 
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ABSTRACT 



An improved system and method for communicating information on a 
transmission medium is disclosed. The disclosed system and method provide improved 
communications on the transmission medium in the presence of interference, which may 
take the form of uncorrelated interference from different services on neighboring 
transmission mediums. In one embodiment, the communications system includes a 
subscriber modem, a central office modem; and a communications channel coupled 
between the subscriber modem and the central office modem. The commimications 
channel transports uplink signals from the subscriber modem to the central office modem, 
and downlink signals from the central office modem to the subscriber modem. At 
frequencies below a selected frequency Me2f, the power spectral density of the 
transmitted upUnk signals is proportional to the power spectral density of the transmitted 
downlink signals by a positive scale factor, and at frequencies above Me2f? the power 
spectral density of the uplink signals is limited to one or more uplink frequency bands and 
the power spectral density of the downUnk signals is limited to one or more downlink 
frequency bands that are disjoint from the uplink frequency bands. Preferably, the total 
bandwidth of the uplink frequency bands above Me2f is proportional to the total 
bandwidth of the downlink frequency bands above Meif by the same positive scale factor. 
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APPENDIX 



Source code describing one preferred embodiment of the invention. 



function [DsRbmax, Sfopt, fopt, UsRbmax] ^ psdestim (LoopNo^ df, symrate, 
start f, Dsfs^ Usfs, maxbits, margin, PdBmin, agn, . . . 
ADSL, HDSLnext, ISDNnext, Tlnext, self, UporDn) ; 



10 



15 



20 



25 



30 



35 



The function 'psdestim' estimates the capacity (datarates) of an xDSL 
loop given the average input power, performance margin, 
the channel noise and the crosstalk sources. 



-5 

9; 



Input : 

LoopNo 
df 

symrate 
startf 
Dsfs : 
Usfs : 
maxbits 
margin 
PdBmin : 



CSA Loop Number 

Bandwidth (kHz) of each subchannel 
DMT symbol rate 
: Start frequency (kHz) of the channel 
Stop frequency (kHz) of the downstream channel 
Stop frequency (kHz) of the upstream channel 

Max. number of bits per carrier (subchannel) 
Performance or noise margin (dB) 



Total average Input Power (dBm) 
agn : Additive gaussian noise (dBm/Hz) 

ADSL : Number of interfering (NEXT+FEXT) ADSL lines 
HDSLnext: Number of interfering (NEXT) HDSL lines 
ISDNnext: Number of interfering (NEXT) ISDN lines 
Tlnext : Number of interfering (NEXT) Tl lines 

self : Number of self-interfering (NEXT+FEXT) lines 
UporDn : Upstream =1, Downstream = 0 
Output : 

DsRbmax : Downstream maximum optimal bit rate 
Sfopt : Transmit spectrum of optimal scheme 
fopt : Frequency vector used by optimal scheme 
UsRbmax : Upstream maximum optimal bit rate 



% Example: 

% psdestim(206,2,2,l, 1000, 1000, le5, 18 . 39, 20, -14 0, 0, 0, 0, 0, 39, 2) ; 

% 

40 % Rohit Gaikwad 

% Rice University 
% 2/2/98 



45 



50 



% Set Fontsize and linewidth 

fsize = 14; 

Iwidth = 1.5; 

largefonts (fsize, Iwidth) ; 



55 



% Nomenclature 
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DS: Different service 



10 



15 



% DS IN-NEXT 

% DSIN-FEXT 

% self-NEXT 

% self-FEXT 

% Mo 

% Ml 

% Mac 

% Mbc 



DS NEXT into line under consideration 

DS FEXT into line under consideration 

NEXT from neighboring lines carrying same service 

FEXT from neighboring lines carrying same service 

M_E2F 

M_E 

M_E2MFDS 
M FDS2MFDS 



20 



% Compute magnitude squared channel transfer function {|Hc(f)|^2) 

% Ref: ADSL standard Tl. 413-95, Annex E 

T = 1/ (symrate*le3) ; % DMT symbol period 

DsW = Dsfs/2; % downstream bandwidth 

N = 1+f loor ( (DsW-startf ) /df ) ; % number of subchannels 



25 



30 



35 



40 



[Hfsqr, freq, len] = cablemodel (LoopNo, N, startf, start f+ (N-1 ) *df) ; 

% compute the channel transfer function 
endf = max(freq); % Last tone frequency 

P = 0.5*le-3*10^ (PdBmin/10) ; % Input power (Watts) on pos freq. range 
% 

% 

% 

% 

% 

% 

% Do DSIN-NEXT and DSIN-FEXT interference computations 
% Ref: ADSL standard Tl. 413-95, Annex B 

% 

% 

% Initialization 

Phdslnext = 0; Pdslnext = 0; Ptlnext = 0; Padsl = 0; 



45 



50 



55 



% DSIN-NEXT+DSIN-FEXT due to ADSL 
if (ADSL 0) 

[finkHzF, Padslfext] = adslf extmain (ADSL, LoopNo) ; 

[finkHzN, Padsl, Padslnext] = adslUsnextmain (ADSL) ; 

[fadslF, Padslfext] resamp ( finkHzF, Padslfext , df, startf, endf ) ; 

[fadslN, Padslnext] = 
resamp (finkHzN, Padslnext, df , f inkHzN (1) , max (f inkHzN) ) ; 

indlow = length ( find (fadslF < fadslN(l))); 

indup = length (find (fadslF > max ( fadslN) )) ; 

Padslnext = [ zeros ( 1, indlow-1) Padslnext zeros ( 1, indup) ] ; 



Padsl 



= Padslnext + Padslfext; 



end 



60 % DSIN-NEXT due to HDSL 

% 2B1Q signaling 
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if (HDSLnext 49) 

load /home/rohitg/Research/Nortel/HDSL_NEXT/hdslnext49 .mat f inkHz 
Phdslnext 

elseif (HDSLnext == 20) 
5 load /home/rohitg/Research/Nortel/HDSL_NEXT/hdslnext20 .mat f inkHz 

Phdslnext 

elseif (HDSLnext 10) 

load /home/rohitg/Research/Nortel/HDSL_NEXT/hdslnextlO .mat f inkHz 
Phdslnext 
10 elseif (HDSLnext == 1) 

load /home/rohitg/Research/Nortel/HDSL_NEXT/hdslnextl .mat f inkHz 
Phdslnext 
end 

15 if (HDSLnext 0) 

[fhdsl, Phdslnext] = resamp ( f inkHz, Phdslnext , df, startf, endf ) ; 
end 



o — — — — — — 

g, 
o 

% DSIN-NEXT due to ISDN 
% 2B1Q signaling of basic access DSLs 
25 if (ISDNnext ==4 9) 

load /home/rohitg/Research/Nortel/ISDN_NEXT/isdnnext4 9.mat f inkHz 
Pisdnnext 

elseif (ISDNnext == 24) 

load /home/rohitg/Research/Nortel/ISDN_NEXT/isdnnext24 .mat f inkHz 
30 Pisdnnext 

elseif (ISDNnext == 10) 

load /home/rohitg/Research/Nortel/ISDN_NEXT/isdnnextlO .mat f inkHz 
Pisdnnext 
end 

35 

if (ISDNnext -= 0) 

[fdsl, Pdslnext] = resamp ( f inkHz, Pisdnnext , df, startf, endf); 
end 

% 

40 % 

% 

% 

% DSIN-NEXT due to Tl 
45 % Same binder Tl 

if (Tlnext == 49) 

load /home/rohitg/Research/Nortel/Tl_NEXT/tlnext4 9.mat f inkHz Ptlnext 
elseif (Tlnext == 25) 

load /home/rohitg/Research/Nortel/Tl_NEXT/tlnext25.mat f inkHz Ptlnext 
50 elseif (Tlnext == 24) 

load /home/rohitg/Research/Nortel/Tl_NEXT/tlnext24 .mat f inkHz Ptlnext 
elseif (Tlnext == 10) 

load /home/rohitg/Research/Nortel/Tl_NEXT/tlnextlO .mat f inkHz Ptlnext 
end 

55 

if (Tlnext -= 0) 

[ftl, Ptlnext] = resamp (f inkHz, Ptlnext, df,startf, endf ) ; 
end 



60 



Conley, Rose & Tayon, P.C. 



A-3 



Atty. Docket No. 1789-01910 



% self-NEXT + self-FEXT 
if (self 0) 
5 % keyboard; 

% Ref: 2-piece Unger model from TlEl . 4 /95-127 
Ntr = [49 10 1] ; 
chilr - [4,8e-9 7.1e-10 6.6e-ll]; 
chihr = [2.4e-13 9,5e-14 8.9e-15]; 
10 explr = [0.4 0.5 0.6] ; 

exphr - [1.4 1.4 1.5] ; 

chil = interpl (Ntr, chilr, self, ^linear'); 
chih = interpl (Ntr, chihr, self, 'linear'); 
15 expl = interpl (Ntr, explr, self, 'linear'); 

exph = interpl (Ntr, exphr, self, 'linear'); 

Hnextsqr=(freq<=20) *chil. * { (freq^leS) .^expl) +. . . 
( f req>2 0 ) *chih . * ( ( f req^leS ) . ^exph) ; 
20 % Reduce Hnextsqr to 1% of its value in case of '^GDSL'' 

% Hnextsqr = 0 . 01 . *Hnextsqr ; 

Hfextsqr = ... 

( ( (3.083e-20) ./ (10^0.6) ) ^ (self^ (0-6) ) *len* ( (leB^freq) .^2) .^Hfsqr) ; 
Hxfsqr = Hnextsqr + Hfextsqr; 

25 else 

% Zero self-interference if no self interfering lines 
Hxfsqr = zeros ( size (Hfsqr) ) ; 

Hnextsqr= zeros (size (Hfsqr) ) ; 
Hfextsqr= zeros ( size (Hfsqr )) ; 
30 Mo = 0; % EQPSD is optimal, no FDS required 

end 



35 



% Add-up DS-interference and channel noise (AGN) 
40 interf - Ptlnext + Pdslnext + Phdslnext + Padsl; % Total DS- 

interference 

No = 2*le-3*10^ (agn/10) ; % AWGN in W/Hz 

noisePSD= ( {No/2 ) +interf ) ' ; % in W/Hz 

noisePSD= noisePSD. *ones ( size ( f req) ) ; % convert into a vector 

45 

df = {freq(2)-freq(l) ) *le3; % convert df into Hz 

SNRgap = 9.8; % dB % factor which allows one to compute 

% achievable bit rates instead 
% of capacity 

50 gamma = 10^ ( (margin + SNRgap) /lO); % Total performance margin (units) 

noisePSD = gamma*noisePSD; % weighted Noise PSD 

Hxfsqr = gamma^Hxf sqr ; % weighted total self- 

interference 

Hnextsqr = gamma* Hnextsqr; 
55 Hfextsqr = gamma ^Hfextsqr; 

H = Hfsqr; 

X = Hnextsqr; 

F = Hfextsqr; 

60 % 
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10 



15 



20 



25 



30 



% Initialization 






pht5tyi r 1 • 








L^eqpsQ — 


LJ f 




% Lagrange multiplier EQPSD soln, with self-interf. 


uwgn — [ J , 




Q. 


Water-filling with no self-interference 


cx - [ J ; 




■6 


self-interference limit when using EQPSD 






o 
-6 


FDS signalling in all subchannels 


Pwgn = 


L J ; 






Bits = [ J ; 








Bitsx = [ ] ; 








Bitsras = L J / 








Bitseqpsd = 


[]; 






DsRbmaxnoise 






[]; % Bit rate doing EQPSD water-filling with no 






% 


self-interf erence 


DsRbmaxx = 


[]; 


% 


Bit rate for self-interference limit 






% 


self-interference dominates AGN and DS-interf erence 


DsRbmaxeqpsd 






[]; % Bit rate doing EQPSD with self-interference 


DsRbmaxfds = 


[]; 


Q, 
O 


Bit rate doing FDS in all subchannels 


DsRbmaxopt = 


[]; 




Bit rate doing EQPSD/FDS optimally 


SNR = []; 








Mzeros = 


[]; 






Bracket = 


[]; 






fMo = []; 








w = noisePSD. 


*df; 


W = 2* (endf-startf ) ; 


% 



35 



%9.3.g,g,9.g,g.3,g.3,9,g,9,^5.g.g.aQ,aaaQ.aQ.Q,aQ.Q.Q,Q.o,a o o o oo o o o o o o o o o o 

oo oo OOOOOO oo O oo 0 oo O oo O'OOOOOOOOOO OOO OOOOOOO OOOO'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'S'S'S'S'S'S'S'S'S 



% Lagrange multiplier solution for the transmit spectrum using EQPSD 
% signaling in presence of self-NEXT+self-FEXT 

% 

40 % Ref: 'Achievable Information Rates on Digital Subscriber Loops: 

% Limiting information rates with crosstalk noise' — J.T, Aslanis and 
% J.M. Cioffi in IEEE Trans, on Commun., vol. 40, no. 2, Feb. 1992 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%^ 



45 



50 



55 



if (Hxfsqr) 

% Using bisection method to calculate the correct value of lamda 
lamdal = (8.06)*lel4; 
lamdau le-4; 

a = Hxf sqr. (Hxfsqr + Hfsqr) ; 

b = w.*(2^Hxfsqr + Hfsqr); 

c = w.*(w - (W.^Hfsqr) ./ (lamdal*2*N) ) ; 

Pestl = (-b + sqrt((b.^2) - 4*a . *c) ) . / (2*a) ; 

Pestl (find(Pestl<0) ) = 0; 

c = w.*(w - (W. *Hfsqr) . / (lamdau*2*N )); 

Pestu = (-b + sqrt((b.^2) - 4*a . *c) ) . / {2*a) ; 
Pestu{find(Pestu<0) ) = 0; 



60 



Pest = Pestle- 
while ( (abs (sum (Pest) - P) > {le3^eps))& 
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(abs (lamdal-lamdau) > leO*eps) ) 

lamda = (lamdal+lamdau) /2 ; %Search at the midpoint 
c =3 w.*(w - (W.^Hfsqr) ./ (lamda*2*N )); 

Pest = (-b + sqrt(b.^2 - 4 . ^a . *c) ) . / (2 . *a) ; 
5 Pest (find (Pest<0) ) = 0; 

% Set the new boundaries for the bisection to work 
if (sum(Pest) > P) 

lamdau = lamda; 
10 elseif (sum(Pest) <= P) 

lamda 1 = lamda; 
end % end if 
end % end while 

% Bisection method converges to yield correct value of lamda 



15 



30 



55 



9'9-9-S-9-&&-9-e-&-9-9-9-9-&-2'3-2-9-9-9-S-9-&S-9-S-9-9-9.5.9.9.&&9.9.S-Q.Q.Q.a o o o o o o o o o o o o o o o o o 



% Computation of capacity and achievable bit rate 

Ceqpsd = [Ceqpsd sum (df*log2 ( 1+ { Pest . *Hfsqr) ./ (w-f- (Hxf sqr . *Pest ) ) ) ) ] ; 
20 Bitseqpsd = [Bitseqpsd min(maxbits, (log2 ( 1+ (Pest . *Hf sqr ) . / . , . 

(w + (Hxfsqr. *Pest) ))))]; 
DsRbmaxeqpsd = [DsRbmaxeqpsd ( 1/T) *sum (Bitseqpsd) ] ; 

else 

25 Pest = zeros (size (freq) ) ; 

DsRbmaxeqpsd = 0; 
end % fi Hxfsqr 

9-9-9-9-9-9-9-&'&-9-9'&-S-&'9-9-5-9-9-9-9-9'9'9-9-&9-B-S-&-5-S-2-9-9.5.9.9.9.9.9.9.9-9-9.9.g-9-g.9.&.9,o o o o o o o o o o o o o o o o o o o o o o o 
o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o'5'5'S'5'6'6'o'6'6'6'6'6'5'5'3'6'6'6'6'5'6'6'6'6"6'5'5'6'6'5'6'6'6'6 

%% 

S-&-9-&-&'B-9-9-5-9r9-9-9-2-Sr9-Sr9-9-S-9-5-&'&'S-S-0-9'9-9-9-9-9-S-9'&-^9'9'9-9'9'9-e-9-5-S-9-Q-9.9.9-Q-e.^a.S.o o o o o o o o o o o o o o o o o o 
%% 

% Transmit spectrum calculation using FDS in all the bins 

35 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
%% 

if (Hfextsqr) 

noisefds= noisePSD; 
[Rbmax,Sfds] = 

40 fdsxtaik (P, Hfsqr, noisefds, Hfextsqr, N,maxbits, symrate*le3, df ) ; 

else 

noisefdsup = interpl (1 :N, noisePSD, 1 : (N-1) / (2*N-1) :N) ; 
Hfsqrup = interpl ( 1 : N, Hf sqr , 1 : (N-1 ) / (2*N-1) : N) ; 

45 % Using bisection method to calculate the correct value of lamdafds 

lamdafdsL = le-14; 
lamdafdsU = le4; 

SfdsL = lamdafdsL - ( (noisefdsup) . /Hfsqrup) ; 
SfdsL(find(SfdsL<0) ) = 0; 
50 SfdsL = down (SfdsL, 2, 0) ; 

SfdsU = lamdafdsU - ( (noisefdsup) , /Hfsqrup) ; 
SfdsU(find(SfdsU<0) ) - 0; 
SfdsU = down (SfdsU, 2, 0) ; 



Sfds = SfdsL; 



while ( (abs ( {df/2) *sum(Sfds) - P) > (lel*eps))& ... 
(abs (lamdafdsL-lamdafdsU) > leO*eps) ) 
60 lamdafds = (lamdaf dsL-t-lamdaf dsU) /2 ; %Search at the midpoint 

Sfds = lamdafds - ( (noisefdsup) . /Hfsqrup) ; 
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Sfds (find(Sfds<0) ) = 0; 

Sfds = down (Sfds, 2, 0) ; 

% Set the new boundaries for the bisection to work 
5 if ( {df/2) *suni(Sfds) > P) 

lamdafdsU = lamdafds; 
elseif ( {df/2) *sum(Sfds) <- P) 

lamdafdsL = lamdafds; 
end % end if 
10 end % end while 

% Bisection method converges to yield correct value of lamdawgn 

if (mod (length (Hf sqrup) , 2 ) ) 

Sfds = up (Sfds, 2,0) ; 

15 else 

Sfds - up(Sfds, 2, 1) ; 

end 

end % fi Hfextsqr 

20 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

g, Q, 
o o 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
%% 

25 % Transmit spectrum calculations for SUBOPTIMAL scheme with self- 

interference 

% Do EQPSD in bins [1,M1] and FDS in bins [M1,K] where Ml is obtained 
% by computing the inequalities as listed below 

ooooooooooooooooooooooooooooooooooooooooooooooooooo'5'6'5'8'3'S'6'6'5"6"6'6'3'6'6'5'6'6'6'6'6'6'6'6 

30 %% 

if (Hxfsqr) 

NumF = (H - 2* (X-F) ) ; 

DenF = ((X.^2) - (F.^2) - (H.*F)); 
35 Fxt = NumF. /DenF; 

Fxtdec = 0 , Stones (size (freq) ) ; 

Fxtindeqpsd = find((NumF> 0) & (DenF <0) ) ; 

Fxtdec (Fxtindeqpsd) = zeros (size (Fxtindeqpsd) ) ; 

Fxtindfds = find( (NumF <0) & (DenF > 0)); 
40 Fxtdec (Fxtindfds) ones (size ( Fxtindfds) ) ; 

if (isempty (Fxtindfds ) ) 

istar = length (Fxtdec) ; 

else 

45 istar = min (Fxtindfds) ; 

end 

Ml = max (Fxtindeqpsd) ; 
if (isempty (Ml) ) 
Ml = 0; 

50 end 

optin - [0 le-12*P] ; 

% Water-filling or peak-constrained water-filling to optimize 
% power distribution given the signaling strategy in each bin 
55 [PI optOut]=fmin('rbmax',0,P,optIn,Ml,P,freq, . . . 

Hf sqr, Hnextsqr, Hfextsqr, noisePSD,maxbits, T) ; 
% Compute the achievable bit rate and the optimal transmit spectrum 
[DsRbmaxeqpsdfds, Sfeqpsdfds]= rbmaxpsd ( PI, Ml, P, freq, , . . 

Hf sqr, Hnextsqr, Hfextsqr, noisePSD, maxbits , T ) ; 
60 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
%% 
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%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
%l 

% Transmit spectrum calculations using OPTIMAL scheme of 
5 % EQPSD in bins [l,Mo] and 

% FDS in bin [Mo, K] 

3-3-&-5-9-9'9'9-9-9'9-9-9-9-&-9-&-&'5-&-9-9-9'9-9-9-9-9-9-&-9-3-5.5-9'9-S.Q-Q-9.o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o 
o oo 0 ooooo oo o oo ooooooooo o oooooooooooooo oo'6'6'6'S'6'6'6'6^'6'6'6'6'6'6'6'6'6'6'6'6'6'6'6"5'6'6'5'6'6'6'6'5'5'6 

% Using Golden Section Search in One dimension 

10 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
%% 

Cap = [] ; Sfopt = [] ; 

R - 0.61803399; % Golden Ratios 

C = 1-R; 

15 

MoL = Ml; % make use of the new knowledge that 

% all bins [1,M1] do EQPSD 
MoH istar; % All bins [istar,K] do FDS 

else 

20 istar - 0; 

end 

if (Hxfsqr) 

MoM = round { (MoL+MoH) /2) ; 

25 optin = [0 le-12*P] ; 

% We use the four point method where the min is bracketed by 
% the three abscissas MoL< MoM < MoH and the capacity 

% CMoM > CMoL and CMoM > CMoH, i.e. the oridnate of MoM is more than the 
30 other 

% two ordinates. 

% Ref: 'Numerical Recipes in C -Cambridge University Press, 2nd ed. pp. 
397-402 

35 xO - MoL; 

x3 = MoH; 

if (abs (MoH-MoM) > abs (MoM-MoL) ) 
xl ~ MoM; 

x2 = round (MoM+C* (MoH-MoM) ) ; 

40 else 

x2 = MoM; 

xl - round (MoM-C* (MoH-MoM) ) ; 

end 

45 % Water-filling or peak-constrained water-filling to optimize 

% power distribution given the signaling strategy in each bin 



50 



[Pll optOutl] =fmin ( 'rbmax' , 0, P, opt In, xl, P, freq, , . . 

Hf sqr, Hnextsqr, Hfextsqr, noisePSD, maxbits , T) ; 



[P12 optOut2]=fmin('rbmax',0,P,optIn,x2,P,freq, . . . 

Hfsqr, Hnextsqr, Hfextsqr, noisePSD, maxbits, T) ; 

[xO xl x2 x3] ; 
while (abs (x3-x0) > 1) 
55 if ( optOut2(8) < optOutl (8) ) 

[xO,xl] = shft2 (x0,xl,x2) ; 
x2 = floor (R*xl+C*x3) ; 

[PI optOutN2]=fmin('rbmax',0,P,optIn,x2,P,freq, . . . 

Hfsqr, Hnextsqr, Hfextsqr, noisePSD, maxbits , T) ; 
60 [optOutl, optOut 2] = shft2 (optOutl, optOut2,optOutN2 ) ; 

else 
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[x3,x2] == shft2 (x3,x2,xl) ; 
xl = floor (R*x2+C'^x0) ; 

[PI optOutNl]=fmin('rbmax\0,P,optIn,xl,P,freq, . , . 

Hf sqr, Hnextsqr, Hf extsqr, noisePSD^maxbits, T) ; 
5 [optOut2,optOutl] = shft2 (optOut2,optOutl,optOutNl) ; 

end 

[xO xl x2 x3] ; 
end 

10 % Output the best of the two current values between xl and x2 

if (optOutl{8) < optOut2{8)) 
Mo = xl; 

opt Out = optOutl; 

else 

15 Mo - x2; 

optOut = optOut2; 

end 

Bracket = [Bracket; xO xl x2 x3] ; 
20 Mzeros = [Mzeros Mo] ; 

if (Mo) 

fMo - [fMo freq(Mo) ] ; 

end; 

DsRbmaxopt = [DsRbmaxopt -optOut(8)]; 



25 



30 



50 



else 

DsRbinaxopt = 0; 
DsRbmaxeqpsdfds = 0; 
end % fi (Hxfsqr) 



35 % Transmit spectrum calculations in presence of no self-NEXT or self-FEXT: 

% EQPSD signaling 

% Ref: 'The Multitione channel' — Irving Kalet, IEEE trans on comm. 
% vol. 37, no. 2, Feb. 1989 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
40 %% 

% Using bisection method to calculate the correct value of lamdawgn 
lamdaL = le-14; 
lamdaU = le4; 
45 SfL = lamdaL - (noisePSD. /Hf sqr ) ; 

SfL(find(SfL<0) ) - 0; 



SfU = lamdaU - (noisePSD . /Hf sqr ) ; 
SfU(find(SfU<0) ) = 0; 

Sf = SfL; 



while ( abs (df *sum(Sf ) - P) > (le2*eps) ) 

lamdawgn = (lamdaL+lamdaU) /2; %Search at the midpoint 
55 Sf = lamdawgn - (noisePSD . /Hf sqr) ; 

Sf (find(Sf<0) ) = 0; 
fsubA = find{Sf > 0) ; 

% Set the new boundaries for the bisection to work 
60 if (df*sum(Sf ) > P) 

lamdaU = lamdawgn; 
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elseif {df*sum(Sf) <= P) 

lamdaL = lamdawgn; 
end % end if 
end % end while 

5 % Bisection method converges to yield correct value of lamdawgn 




% Computation of capacity and achievable bit rate 

10 

PdBm = [PdBm lO^loglO (2*P/le-3) ] ; 

SNR = [SNR 10*loglO{P./(noisePSD.*DsW) ) ] ; 

Cwgn = [Cwgn df '^sum (log2 (1 + (Sf . *Hf sqr) . / (noisePSD) ) ) ] ; 

Pwgn = [Pwgn df *sum (Sf ) ' ] ; 
15 Bits = [Bits min(maxbits, (log2(l+(Sf .^Hfsqr) ./noisePSD) ) ) ] ; 

% Bits per transmission 

% in each subchannel or bits /symbol 

DsRbmaxnoise = [DsRbmaxnoise (1/T) *sum(Bits) ] ; 

% Total bit rate in bps 
20 % T is the DMT symbol period 

g,9,g,5,g,5,g,5.g,^^g,^g,5,5,^g,5,g,3,3,g,g,g,g,C[,5.g,g,^g,g.g.2,g,g,g,g,g,g,^^^Q,o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o 
ooooooooooooooooooooooooooooooooooooooooooooo'5'6'6'6'6'6'6'6'S'5'S'6'6'6'6'6'3'5'6'6'o"6"6'5'6'6'5'S'6'6 

%% 

% Peak PSD constraint in presence of AGN: EQPSD solution 
25 % 5/14/98 

[DsRbmaxnoiseP,SfP,SfOPTIS] = peakpsdwf (noisePSD, Hf sqr, P, df ,maxbits, . . . 
T, freq, UporDn) ; 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

Q, Q, 

30 

% Capacity crosstalk limit 
if (Hxfsqr) 

Cx = sum( (W/ (2*N) ) *log2 (1 + (Hf sqr . /Hxfsqr) )) ; 
Bitsx = [Bitsx min (maxbits^ (log2 (1 + Hf sqr , /Hxfsqr) ))] ; 
35 DsRbmaxx= [DsRbmaxx ( 1/T ) *sum (Bitsx) ] ; 

end 

2-9-&-S-9-5-9-9'9'9-9-9'9-9-9-&-&-5-&-$-9-9-S-9'9-9-&-9-&'S-9-9.S.9.9.B.5.9-9.9-9.9.9.9-0 o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o 
ooooooooooooooooooooooooooooooooooooooo'6'6'5"6'6'5'6'6'6'6'6'6'5'6'5'6'6'6'6"6'6'6'5'6'S'6'6'6'6'6'5'6'6'6'6'6 

%% 

40 % Capacity of scheme using FDS through all the bins 

if (Hfextsqr) 
Cfds = [Cfds 

0.5*df*sum{log2(l+(Sfds* .*Hfsqr) ./{Sfds' . *Hf extsqr+noisef ds ) ) ) ] ; 
Bitsfds = [Bitsfds min (maxbits, 0 . 5* (log2 ( 1+ (Sfds * . ^Hf sqr) . / . . . 

45 (Sfds' .^Hfextsqr + noisefds) ) ) ) ] ; 

% Bits per symbol 
DsRbmaxfds = [DsRbmaxfds (1/ (T) ) *sum{Bitsf ds) ] ; 
else 

Cfds = [Cfds {df/2) *sum(log2 (1 + (Sf ds . *Hf sqrup) . / {noisef dsup) ) ) ] ; 
50 Bitsfds = [Bitsfds min(maxbits, 0.5*(log2(l + 

(Sfds. ^Hfsqrup) , / (noisefdsup) ) ) ) ] ; 

% Bits per symbol 

DsRbmaxfds = [DsRbmaxfds ( 1/T) *sum (Bitsfds) ] ; % bit rate in bps 

Sfds = down (Sfds, 2, 0) ; 
55 end 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%!%%%% 
%% 

60 %% 

% Upsample all the transfer functions 
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Hfsqrup - interpl ( 1 : N, Hf sqr, 1 : (N-1 ) / (2*N-1 ) : N) ; 

Hxfsqrup = interpl { 1 : N, Hxf sqr , 1 : (N-1) / (2*N-1) :N) ; 

Hnextsqrup = interpl (1 :N, Hnextsqr, 1: (N-1) / (2*N-1) :N) 

Hfextsqrup = interpl ( 1 : N, Hfextsqr, 1 : (N-1) / (2^N-1) :N) 

5 noisePSDup = interpl { 1 : N, noisePSD, 1 : (N-1 )/ {2*N-1) : N) 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
%% 

% Marshalling output arguments 
10 if (self) 

UsRbmax = 0; 

[DsRbinax, Sfopt]= rbmaxpsd ( PI, Mo, P, f req, . . . 

Hf sqr, Hnextsqr , Hf extsqr, noisePSD, maxbits, T) ; 

fopt = freq; 
15 fbit = DsRbmax; % Bit rate in bps 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
%% 

% Upsample the transmit spectra 

Sfupl = [Sfopt (l:Mo) ; Sf opt ( 1 : Mo) ] ; % EQPSD region 
20 Sfupl = Sfupl ( : ) ; 

Sfup2 - [Sfopt (Mo+1: end) ; zeros ( 1, N-Mo) ] ; % FDS region 

Sfup2 = Sfup2 ( : ) ; 

"2 25 Sfoptup = [Sfupl; Sfup2] ; % Upsampled correct transmit spectra 

foptup = interpl (1:N, fopt, 1: (N-1)/ (2*N-1) :N) ; 

dfl = df/2; 
M fupl = foptup (1: length (Sfupl) ) ; 

O fup2 = foptup{l+length{Sfupl) rend) ; 

r^f: SfoptupDn = [Sfupl;cshift (Sfup2' , 1) ' ] ; % Upsampled correct downstream 

% transmit spectra 

- %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
Q 35 %% 

m % Multi-line FDS signaling (Scheme C) 

fU %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

LI 

\Z G = [2^Sfopt (l:Mo) ' ./noisePSD(l:Mo) ; Sfopt (Mo+1 : N) /noisePSD (Mo+1 : N) ] ; 

U 40 exprBC = {(1+G.*H) - ( (1+ ( (self+1) /2) ^G. *H) , ^ (2/ (self+1) )))./.. . 

□ (G.* ( (l+( (self+1) /2) *G.*H) , ^ ( 2/ ( self+1 ) ) -1) ) ; 

indSchmBC = f ind ( F<exprBC) ; % FDS (Scheme B) in these bins 

if (-isempty (indSchmBC) ) 

Mbc = max (max (indSchmBC) , Mo) ; 

45 else 

Mbc = Mo; 

end 

% Lets keep the same G 
50 exprAC = ( (2+G. * (X+H) )-( (1+ ( (self+1) /2) *G. *H) . (1/ (self+1) ). ^ (2+G. *X) 

))./... 

(G.*((l+( (self+1) /2)^G.^H) . ^ (1/ (self +1) ) -1 ) ); 
indSchmAC = f ind (F<exprAC) ; % EQPSD (Scheme A) in these bins 

if { ^is empty (indSchmAC) } 
55 Mac = max (indSchmAC) ; 

else 

Mac - 0; 

end 

if (Mac > Mo) 
60 Mac = Mo; 

end 
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[DsRbmaxCSfoptC] = rbmaxpsdC {PI, Mo, P, freq, . . . 

Hfsqr,Hnextsqr,Hfextsqr,noisePSD,maxbits, T,Mac,Mbc, self) ; 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%^^%^^%^^^^^^^^^^^^ 

% Grouping of bins: Contiguous spectra (several ways to optimally group) 



% Optimal grouping: equal performance margins in both transmission 
10 % directions 

Sfdsupdn = [Sf opt (Mo+1 rend) ; Sf opt (Mo+1 : end) ] ; 
Sf dsupdn = Sfdsupdn ( : ) ; 

15 Bitsopteqpsddn= [log2 ( 1+ (Sf opt ( 1 :Mo) • . *Hf sqr ( 1 :Mo) . . 

(Sfopt (l:Mo) » .*Hxfsqr(l:Mo} + noisePSD (1 :Mo) )))]; 
Bitsopteqpsddn = [Bitsopteqpsddn ^' Bitsopteqpsddn ' ] ; 
Bitsopteqpsddn = Bitsopteqpsddn (:) ; 

Bitsoptfds = [log2 (1+ (Sfdsupdn* . *Hf sqrup {2*Mo+l : 2*N) )./... 
20 (Sfdsupdn' .*Hfextsqrup(2*Mo+l:2*N) + noisePSDup (2*Mo+l : 2*N) ) ) ] ; 

Bitsopt = [Bitsopteqpsddn* Bitsoptfds]; 

Bitsoptcum = cumsum (Bitsoptfds ) ; 
eqpsddn = sum (Bitsopteqpsddn) ; 
fds = sum (Bitsoptfds ). /2; 
25 fds+eqpsddn; 

if (^isempty (Bitsoptcum) ) 

Mg = 2*Mo + max (find (Bitsoptcum < (Bitsoptcum (end) /2 ))) ; 
else 
Mg = 2*Mo; 
30 end 

sum (Bitsopt (2*Mo+l:Mg) ) ; 
sum (Bitsopt (Mg+1:2*N) ) ; 

SfoptUpDn = [Sfupl; Sfdsupdn]; 
35 SfoptsmthUp = [Sfupl' Sf dsupdn ( 1 : Mg-2*Mo) ' eps*ones ( 1, 2*N-Mg) ] » ; 

SfoptsmthDn = [Sfupl' eps^ones ( 1, Mg-2^Mo) Sfdsupdn (Mg-2^Mo+l : end) ']' ; 



% Optimal grouping: equal performance margins and equal average powers in 
40 % both directions of transmission 

% Note that Mg = 2*Mo + Mc + somemore 

Mc = 114; % for 39 self-NEXT + 39 self-FEXT 
% Mc =60; % for 24 Tl NEXT + 24 HDSL2 NEXT 
45 % Mc = 92 ; % CSA 4: 39 self-NEXT + 39 self-FEXT 

% Mc = 112 ; % for 39 self-NEXT + 39 self-FEXT + 1 HDSL 
% Mc = 113 ; % for 4 9 self-NEXT + 49 self-FEXT -f 1 HDSL 

Sfdscum = cumsum ( Sfdsupdn (Mc+1 : end) ) ; 

50 Mg = 2*Mo + Mc + max ( find (Sfdscum < (Sfdscum (end) +sum (Sfdsupdn ( 1 : Mc) )) /2 

) ) ; 

SfoptUpDn = [Sfupl; Sfdsupdn]; 

SfoptsmthUp = [Sfupl' eps*ones (l,Mc) Sfdsupdn (Mc+1 :Mg-2*Mo) ' ... 
eps^ones (1,2*N-Mg) ] ' ; 
55 SfoptsmthDn = [Sfupl' Sf dsupdn ( 1 : Mc) ' eps^ones (1, Mg-2*Mo-Mc) ... 

Sfdsupdn (Mg-2^Mo+l : end) ' ] ' ; 
PDn = lO^loglO (le3*2*sum(SfoptsmthDn) * (df/2) ) 
PUp = lO^loglO (le3*2^sum(SfoptsmthUp) * (df/2) ) 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%^^^^^^^ 
60 %% 
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^^^^^^^^%:^^^^^5'^9'9'9'^9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9-9'9'9f9'9'9'Q'9'9'9'^Q'^^Q-'3~<i.o o o o o o o o o o o o o o o o o o o o o o 
o o o o o oo o o o o o o o o o-6-6-5-6-6'6'6'6'6'6-6-6-6-6-6'6'5-6-6-6-6'6'6'6'6'5-5-6'6-6^-S-6-6'6'6-6'6-6%^^^^^%%%%%%%%%%%%%% 

% Broadband modulation schemes for the contiguous spectra (specifically 
% multi-level PAM using Decision Feedback Equalizers (DFE) at the receiver 

5 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

9- 9- 
o o 

Sfdsupdnciom = cumsum (Sf dsupdn) ; 
if (--isempty (Sfdsupdncum) ) 

dfeupdnbrk = 2*Mo +max ( f ind (Sf dsupdncum < (Sfdsupdncum (end) /2) ) ) ; 
10 else 

dfeupdnbrk = 2*Mo; 
end 

dfeupdnbrk = Mg; 

15 dfeSfoptUp - [SfoptUpDn (1: dfeupdnbrk) ; zeros (2^N-dfeupdnbrk, 1 )] ; 

dfeSfoptDn - [SfoptUpDn (1 : 2*Mo) ; zeros (dfeupdnbrk-2*Mo, 1 ) ; ... 

SfoptUpDn (dfeupdnbrk+1: 2 *N) ] ; 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

g, Q, 
'6 ■& 

20 % DFE margin estimation of the contiguous spectra 

baud = (l*fbit) / {3*le3) ; 

% baud = foptup (max(find(Sfoptsmthup>0) ) ) / (df/2e3) 

reqdsnr = 27.7; 

% Check DFE margin of smooth PSD using marginest program 
25 BWUp = 2*le3*foptup (dfeupdnbrk) ; % Upstream bandwidth 

BWDn = 2*le3* ( 1+foptup (end) -foptup (dfeupdnbrk) +foptup (2*Mo) ) ; 
% BWDn = 2*le3*foptup(end) ; 

bu = fbit/BWUp; % bits/symbol for Upstream 

bd = fbit/BWDn; % bits/symbol for Downstream 

30 baudUp = fbit/ (df l*bu) ; % symbols/s 

baudDn = fbit/ (df l^bd) ; 
reqdsnrUp = pamreqdsnr [2^bu, le-7); 
reqdsnrDn = pamreqdsnr ( 2^bd, le-7); 

noisePSDdfeUp = noisePSDup+ (Hf extsqrup . * (df eSf optUp M ) + ... 
35 (Hnextsqrup.^ (dfeSfoptDn' ) ) ; 

dfeM_smthUp = dfemarginest (SfoptsmthUp, Hfsqrup, noisePSDdfeUp ' ./gamma, . . . 

baudUp, reqdsnrUp) ; 
noisePSDdfeDn = noisePSDup+ (Hf extsqrup, ^df eSf optDn *) + ... 

(Hnextsqrup. *dfeSfoptUp' ) ; 
40 dfeM_smthDn = dfemarginest (SfoptsmthDn, Hfsqrup, noisePSDdfeDn' . /gamma, .. . 

baudDn, reqdsnrDn) ; 



45 



%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 



% Check DMT margin of contiguous spectra using marginest program 
noisePSDUp = noisePSDup+ (Hf extsqrup. ^Sf optsmthUp' ) + ... 
50 Hnextsqrup . *Sf optsmthDn ' ; 

chkM_smthUp = dmtmarginest (Sf optsmthUp ' , Hfsqrup, noisePSDUp . /gamma, . . . 

SNRgap, maxbits, (2*DsRbmax) / (symrate^le3 ) , 0 ) 
noisePSDDn = noisePSDup+ (Hf extsqrup . *Sf optsmthDn ' ) + . . . 

Hnextsqrup. *Sf opt smthUp' ; 
55 chkM_smthDn = dmtmarginest (SfoptsmthDn ', Hfsqrup, noisePSDDn, /gamma, ... 

SNRgap, maxbits, (2*DsRbmax) / (symrate*le3) , 0) 

%%%%%%^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ o O O O O O O O O O O O O O O O O O O O O O O O O O O O O 0 o o o o 
%% 

60 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
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% No self-interference 
else 

DsRbmax = DsRbmaxnoise; 
Sfopt = Sf; 
5 fopt - freq; 

fbit = DsRbmax; % Bit rate in bps 

foptup = interpl (1:N, fopt, 1: (N-1) / (2^N-1) :N) ; 

dfl = df/2; 

10 Sfoptup = interpl (1:N, Sfopt, 1: (N-1) / (2^N-1) :N) ; 

Sfopt up == Sfoptup ( : ) ; 

chkM_disc=dmtmarginest (Sfopt \ Hfsqr,noisePSD. /gamiua, SNRgap,maxbits, . . . 

DsRbmax/ (symrate^leS) , 0) ; 
bsym = fbit/ (2*le3* (fopt (length (find(Sfopt>eps) )))) ; 
15 baud = fbit/(df*bsym) ; 

reqdsnr = pamreqdsnr (2''bsym, le-7); 

dfeM_disc = dfemarginest (Sfopt , Hfsqr, noisePSD /gamma, . 
baud, reqdsnr) ; 

end 

20 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
%% 

9,9,9,9,9,9,5.9,9,5,2,^9,5,9,9,5,5.9,0,2.0.^5.0.9,9,0,0.9,0,0,0,0,5,0,0^ O O O O 0 o o o o o o o o o o o o o 

ooooooooooooooooooooooooooooooooooooooooooooooooooooo'o'o'o'0"o'o'o'o'o'o'5'6'6'6'6'6'6'6'S'5'6'6 

25 % Check DMT margin of discrete PSD using marginest program 

noisePSDdmtup == noisePSDup+ (Hfextsqrup . ^Sfoptup ^ ) + ... 

[Hnextsqrup(l:2*Mo) . *Sf optup ( 1 : 2^Mo) ' zeros (1, 2*N-2*Mo) ] ; 
%noisePSDdmtup = interpl ( 1 :N, noisePSDdmt, 1 : (N-1) / (2*N-1) :N) ; 
chkM_disc=dmtmarginest (Sfoptup' , Hfsqrup, noisePSDdmtup . /gamma, SNRgap, maxbits 



30 



50 



(2*DsRbmax) / (symrate*le3) ,0) ; 



% baud = {l*l,5520e3) /3; 

% reqdsnr = 27.7; 
35 if (self) 

bsym = fbit/ (2*le3* (foptup (2*Mo) + ( 1+foptup (end) -foptup (2*Mo) ) /2 )) ; 

else 

bsym = fbit/ (2*le3* (foptup(length(find(Sfoptup>eps) ) ) ) ) ; 

end 

40 baud = fbit/ (df l*bsym) ; 

reqdsnr = pamreqdsnr (2^bsym, le-7); 

dfeM_disc = dfemarginest (Sfoptup, Hfsqrup, noisePSDdmtup ' . /gamma, . . . 
baud, reqdsnr) ; 

45 baud - 2*baud; 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
%% 



i-9-9-S'&-9-&-9-&'9-9-5'2-9-9-9-9-&-9'9-5-9-9-9'9-9-S-5-9'9-9-9-9.9.9.5.9,9.9.9.o 9-2-0 5,0 0 o o o o o o o o 0 o o o o o o o 00 o o o 00 o o o 00 
Doooooooooooooooooooooooooooooooooooo'5-6'6'6'6'5'6'6-6'6'6'S'5'6'6'6'6'5'6^'S'6'6'6'6'6'5'5'6'6'6'S'5'6^'6"6'6 



% Print the various margins 

fprintf ( ' \nCheck margins for various modulation techniques \n' ) ; 
fprintf ( 'DiscDMT 

DiscDFE\tContDMTup\tContDFEup\tContDMTdn\tContDFEdn\n' ) ; 
55 fprintf ('%2.2f %2,2f \t %2.2f \t %2.2f \t %2.2f \t 

%2.2f\n' , . . . 

chkM_disc,dfeM_disc,chkM_smthUp,dfeM_smthUp,chkM smthDn,dfeM smthDn) ; 
fprintf ( '\n' ) ; ~ ~ 

9-^9-9-9-9-&-&'S-&'9-5-5-9-5-5-3-9-&.5-5-5-9-@.5.2.9-9.9-9,5.5.2.9,9.9,9,9,5,9,g,9,9,5,5,5,g,o o o o o o o o o o o o o o o o o o o o o o o o o o o o 
o 0 o o o o o o o o o o o 00 o o o o 00 o 0 o o o^ 0 o o o o o'5'5'5'6'6'6'6'5'5'5'5'6'6'6'5'6'6'6'6'6'5'6'S'6'6'5'5'6'6'5'5'5'5'6'6'6'5'6'6'6'5'6 

60 %% 
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% Plot the different PSDs 
5 if(self) 

plot (fupl,10*loglO(le3*Sfupl) ) ; hold on; 
Sfdsupdn(l) = Sfupl{end); 

plot{fup2,10*loglO(le3*Sfdsupdn), 'r~M; hold off; 
legend ('EQPSD signaling', ' FDS signaling' , 0) ; 
10 else 

plot (foptup, lO^loglO (leB^Sfoptup) ) ; 

end 

hold off; grid on; largef onts ( f size, Iwidth) ; 
axis ( [0 500 -54 -34] ) ; 
15 xlabel (' Frequency (kHz)'); 

ylabel ( 'Amplitude (dBm/Hz) ' ) ; 
if (HDSLnext) 

% title ( 'Optimal HDSL2 transmit spectrum with 49 HDSL DSIN-NEXT'); 
elseif (Tlnext) 

20 % title ( 'Optimal HDSL2 transmit spectrum with 25 Tl DSIN-NEXT'); 

elseif (self ) 

title ( 'Optimal HDSL2 transmit spectrum with 39 self-NEXT + 39 self-FEXT'); 
end 

orient tall; 

25 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
%% 

% Print the various Powers 
Qtmp = le-3'^(10.^( (SfOPTIS-1) ./lO) ) ; 

Ptmp = df *sum (Qtmp) ; 
30 PtmpdBm = 10*loglO ( le3*Ptmp) ; 

fprintf('\n Powers (Positive freq. spectrum) for various tranmit 
spectra\n ' ) ; 

f print f ( 'OPTISXt ContigUpXt ContigDn\t (in dBm) \n' ) ; 
35 f print f('%2.2f\t %2.2f\t %2.2f \t \n', PtmpdBm,... 

10*loglO (le3* (df /2) •^sum(SfoptsmthUp) ) , 
lO'^loglO (le3* (df /2 ) *sum (Sf optsmthDn) ) ) ; 
fprintf ('\n' ) ; 

9.^5-9-Q.9.9.9-5.9-9.Q-9.9-B-Q.5-5-9-S-9-9-9-5-9-9-5-9-9-9-5-9-9-9-9.9-5-5.9.5'9.5.Q-5.5.5'5.Q.Q.Q.9.S.Q-Q.^Q.<3-0 o o o o o o o o o o o o o o o o o 
'O'O'&'S'&'O'&'&'O'O'O'O'S'O'O'O'O'O'O'D'S'O'O'O 

40 %% 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
%% 

45 % SPECTRAL COMPATIBILITY 

% Define optional parameters 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
%% 

srcno = self; 

50 cdmaOrfds - 2; % 0 = water-filling, no self-NEXT or self-FEXT 

% 1 = EDS in all bins 
% 2 = EQPSD in 1 to Mo bins and FDS from Mo+1 bin onward 
% 3 = EQPSD in 1 to Mo bins and CDS from Mo+1 bin onward 

&9-9-&&-&S-9-&'9-S-&^&9-9-&'&^9-5-9-9'2-5-9'9-9-9-9-9-&&-9-S-9-9-9-9-9-5.S.9.9.9-9^ 9-0 o o o o o o o o o o o o o o o o o o o o o o 

ooooooooooooooooo'ooooooooooooooooooooooooooooooooo'6o'6'6'6'6'6'6'6'6'6'8"6"6'6'5'6'6'6'6'6'5'6'6'6 



55 



freqH2 = freq; 

dfH2 = freq (2) - freq(l) ; 



if (Hnextsqr) 

60 switch cdmaOrfds % 1 if next exists otherwise 0 

case 0, 
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SfH2 = Sf; 

case 1, 

SfH2 = Sfds; 

case 2f 

5 SfH2 = Sfoptup; 

freqH2 ^ foptup'; 

dfH2 - freqH2(2) - freqH2{l); 

case 3, 

[CH2, SfH2]= rbmaxpsd(Pl,Mo,P,freqH2, . . . 
10 Hf sqr, Hnextsqr, Hfextsqr, noisePSD^maxbits, T) ; 

SfH2(Mo+l:N) = (Sf H2 (Mo+1 : N) . /2 ) ; 
SfH2 (Mo+2:2:N) = Sf H2 (Mo+1 : 2 : N-1 ) ; 
otherwise, 

SfH2 = Sf; 

15 end % switch 

else 

SfH2 - Sf ; 
end % fi 

20 startf = 1; 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
% 

% Spectral compatibility with HDSL 

0^ 5.9'9-&9-5-9'9-2r9-&'5-9'9-9'9-5-9'9-9-9-9-9-&9-9-&-9-9-&&9'&9-9-9-S-&9-9-S-9-9-9-9-9-9-9-Q-5-OOQ,oooooooooooooooooooooo 
^ J ooooooo 0*0000 ooo ooo ooo oooooo o oooo'oo ooooooo 0 ooooooo'oo'5'o'6'o'o'6'o"o'o'o'5'6'S'6'6'6'S'5'5'3'5'6'6'6 

% 

[freqH, Phdsltx, Phdslnext] = hdslnextmain ( 4 9) ; % 49 HDSL DSIN-NEXT power 
baud = 392/dfH2; 
wgn = le-3*10^ (-140/10) ; 
30 reqdsnrhdsl = 21.5; 

f reqH2 = f reqH2 ' ; 

O O 0 o o o o o o o o o o o o o o 

oo'ooooooooooooooooooooooooooo'ooooooooooooooooooooooo'oo'o'o'o'6'6'6'S"6'6'6'6'6'6'6'6'6'6'5'5'6'6 
% 

35 % HDSL2 into HDSL 

% Computing NEXT of HDSL2 lines 
% Ref: 2-piece Unger model from TlEl . 4 /95-127 
Ntr - [49 10 1] ; 
chilr = [4.8e-9 7.1e-10 6.6e-ll]; 
40 chihr = [2.4e-13 9.5e-14 8.9e-15]; 

explr - [0.4 0.5 0.6]; 
exphr = [1.4 1.4 1.5] ; 

chil = interpl(Ntr, chilr, srcno, 'linear'); 
45 chih = interpl(Ntr, chihr, srcno, 'linear'); 

expl = interpl(Ntr, explr, srcno, 'linear'); 
exph = interpl(Ntr, exphr, srcno, 'linear'); 

nextsqr=(freqH2<=20) ^chil.* ( (freqH2*le3) .^expl)+. . . 
50 (freqH2>20) *chih.* ( (freqH2*le3) .'^exph) ; 

[fH2, HfsqrH2] = resamp ( f opt , Hfsqr, dfH2, startf , max ( freqH2 )) ; 
f extsqr = ... 

( ( (3.083e- 

20) ./ (10-^0.6) ) * (srcno^ (0-6) ) *len* ( (le3*freqH2) .^2) ,*HfsqrH2) ; 
55 Phdsl2x = SfoptsmthUp' .^nextsqr + SfoptsmthDn' . *f extsqr; 

% Check for contiguous spectra 
Phdsl2x = Phdsl2x.* (1. / (l+(6500. / (freqH2*le3) ) ) ) ; 

% LPF cutoff frequency =6.5 kHz 
[HfsqrH,freqtH]=cablemodel(LoopNo,l+(max(freqH) ) /df H2, f reqH ( 1) , . . . 
60 max (f reqH)); 

[freqtH, PhdsltxH] = resamp (f reqH, Phdsltx, dfH2, startf, max ( f reqtH) ) ; 
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Phdsl2xH = [Phcisl2x zeros ( 1, length { freqtH) -length ( freqH2 ))] ; 

hdsl2__in_hdsl =df emarginest (PhdsltxH, Hf sqrH, Phdsl2xH+wgn, baud, reqdsnrhdsl) 

B-&-9-B-9-9-6'9-&-9-9-9-S'9'9-9'9-9-9-9-9-S-9-9'&&-9^ o o o o o o o o o o o o o o o o o o o o o o o 

o 0 o o o o o o o o o o o o o o o o o o o o o o o'o'6"5'0'3'6'6'6'5'5'5'5"S-s-6'3'S'6tI'6'6'6'6'S^'6-S-6'6'6'6'6'S"6'6'S-5-6"3'6-^ 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
% 

% Spectral compatibility with echo-canceled ADSL (EC ADSL) 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
% 

% HDSL2 into EC ADSL 
dfadsl = 4.3125; 

adslstartf = 25; 

adslendf = 1104; 

[freqH2A, Phdsl2xA] = resamp (f reqH2 , Phdsl2x, dfH2, adslstartf, 
max (f reqH2 ) ) ; 

[freqadsl, Padsltxhp, Padslnext] = adslnextmain ( 1) ; 

% Note that we use adslnextmain, this does not include the sine rolloff 
% as suggested by modifications in TlEl . 4 /95-127 

[freqA, PadsltxhpA] = resamp (f reqadsl, Padsltxhp, dfH2 , adslstartf , adslendf ) ; 

N = l+round( (adslendf-startf ) ./dfH2) ; 

[HfsqrA, f] = cablemodel (LoopNo, N, startf, adslendf ) ; 

[freqA, HfsqrA] = resamp(f, HfsqrA, dfH2 , adslstartf , adslendf); 

Phdsl2xA = [Phdsl2xA zeros ( 1, length (freqA) -length (freqH2A) )] ; 
hdsl2_in_ecadsl=dmtmarginest (PadsltxhpA, HfsqrA, Phdsl2xA+wgn, . . . 
9,8, 15, 6.784e3/dfH2,0) 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
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